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Abstract 

The NuMI facility at Fermilab will provide an extremely intense beam of neutrinos for the MINOS 
neutrino-oscillation experiment. The spacious and fully-outfitted MINOS near detector hall will be 
the ideal venue for a high-statistics, high-resolution and F-nucleon/nucleus scattering experiment. 
The experiment described here will measure neutrino cross-sections and probe nuclear effects es- 
sential to present and future neutrino-oscillation experiments. Moreover, with the high NuMI beam 
intensity, the experiment will either initially address or significantly improve our knowledge of a 
wide variety of neutrino physics topics of interest and importance to the elementary-particle and 
nuclear-physics communities. 
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1 Executive Summary 



The imminent completion of the NuMI beamUne, which will be the highest intensity neutrino beamline 
in the worid for many years after its completion, offers the particle and nuclear physics community a 
new opportunity. By constructing a fully active neutrino detector to run for the first time in a high rate 
neutrino beam, the MINERz^A experiment, a collaboration between the high energy physics community 
already working at Fermilab and groups of new users from the medium energy nuclear physics com- 
munity, proposes to exploit this opportunity to access a broad and rich program in neutrino scattering 
physics. 

MINERz^A will be able to complete a physics program of high rate studies of exclusive final states 
in neutrino scattering, as described in Chapters I^HEl of elucidation of the connection between pQCD 
and QCD in non-perturbative regime, as described in Chapter [TOl and of studies of the axial current in 
the elastic (Chapter |6ll, DIS (Chapter flOb and off-forward (Chapter fTTTi regimes, as well as inside the 
nucleus (ChapterfT^. MINERz^A then seeks the application of its data to aid present and future neutrino 
oscillation experiments (Chapter fT3l. where understanding the details of neutrino cross-sections and 
final states is essential for separating backgrounds to oscillation from signal. 

MINERz^A can address all these topics, and can bring a new physics focus to the Fermilab program 
with a simple, low-risk detector of modest cost, as detailed in Chapters (141 and fT6HT7l The performance 
of this detector is expected to be excellent for resolving individual final states as well as measuring 
kinematics in inclusive reactions as documented in Chapter (TSl 

As we submit this proposal to Fermilab, we are also preparing to request funding from sources 
outside Fermilab to pursue this interdisciplinary experiment at the intersection of particle and nuclear 
physics. We request that the lab and its Physics Advisory Committee support this physics and our efforts 
to seek outside funding by granting stage one approval to MINERz/A in time for this approval to enter 
into funding deliberations this spring and summer. 
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2 Overview of the MINERz/A Experiment 



Upcoming neutrino oscillation experiments in the United States, Europe and Japan are driving the con- 
struction of new, very intense neutrino beamlines required to achieve reasonable event rates at detectors 
located hundreds of kilometers away. These new beamlines will allow us to initiate a vigorous research 
program at a detector, located close to the production target, where event rates are much higher than 
at the previous generation of neutrino beam facilities. In addition, it is neutrino oscillation experi- 
ments, with their low-energy neutrinos and massive nuclear targets, which highlight the need for much 
improved knowledge of low-energy neutrino-Nucleus interactions, the overall goal of this experiment. 

At Fermilab, the new neutrino facility NuMI, designed for the MINOS neutrino oscillation exper- 
iment, will be based on the Main Injector (MI) accelerator. The neutrino beams from the MI yield 
several orders of magnitude more events per kg of detector per year of exposure than the higher-energy 
Tevatron neutrino beam. This highlights the major improvement of this next generation of neutrino 
experiments. One can now perform statistically-significant experiments with much lighter targets than 
the massive iron, marble and other high-A detector materials used in the past. That these facilities are 
designed to study neutrino oscillations points out the second advantage of these neutrino experiments: 
An excellent knowledge of the neutrino beam will be required to reduce the beam-associated system- 
atic uncertainties of the oscillation result. This knowledge of the neutrino spectrum will also reduce the 
beam systematics in the measurement of neutrino-scattering phenomena. 

To take advantage of these major improvements in experimental neutrino physics possible with the 
NuMI beam and facility, a collaboration of elementary particle and nuclear physics groups and insti- 
tutions named "MINERi/A" (Main INjector ExpeRiment: v-K) has been formed. This collaboration 
represents the combined efforts of two earlier groups that submitted Expressions of Interest (EOI) 10 
to the Fermilab PAC in December, 2002. The goal of the MINERz^A experiment is to perform a high- 
statistics neutrino-nucleus scattering experiment using a fine-grained detector located on-axis, upstream 
of the MINOS near detector. 

2.1 The Fermilab NuMI Facility 

The Fermilab NuMI on-site facility is made up of the beamline components, the underground facilities 
to contain these components and a large, on-site experimental detector hall to contain the MINOS near 
detector, located just over 1 km downstream of the target and ~ 100 meters underground. 

2.1.1 The NuMI near experimental hall 

This experimental hall is being constructed and completely outfitted for the MINOS near detector. The 
hall is 45 m long, 9.5 m wide and 9.6 m high. There is a space upstream of the MINOS near detec- 
tor amounting to, roughly, a cylindrical volume 26 m long and 3 m in radius for additional detector(s) 
which, were it desired, could use the MINOS near detector as an external muon-identifier and spec- 
trometer. 

2.1.2 The NuMI neutrino beam 

The neutrino energy distribution of the NuMI beam can be chosen by changing the distance of the target 
and second horn from the first horn, as in a zoom lens. The energy of the beamline can also be varied, 
essentially continuously, by simply changing the target's distance from the first horn and leaving the 
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second horn in a fixed position. There is a loss of event rate with this procedure compared to also 
moving the second horn, and the most efficient energy tunes will always require moving the second 
horn. However, moving the target and second horn involves considerably more time and expense then 
simply moving the target. It is now expected that the Main Injector will deliver 2.5 x lO^'^ POT/year at 
the start of MINOS running, and will ramp up to to higher proton intensities if the required funds can be 
obtained. The charged-current neutrino event rates per ton (of detector) per year at startup of MINOS 
would then range from just under 200 K to over 1200 K depending on the position of the target. 

For the MINOS experiment the beamline will be operating mainly in its lowest possible neutrino 
energy configuration to probe the desired low values of Am^. However, to minimize systematics, there 
will also be running in higher-energy configurations that will significantly increase the event rates and 
kinematic reach of MINERz/A. 

The f e content of the low-energy beam is estimated at just over 1 % of the flux. An important 
function of MINERz/A will be to provide a far more accurate measurement of the f e flux and energy 
spectrum within the NuMI beam than is possible with the much coarser MINOS near detector. This 
important figure-of-merit is needed for the design of next-generation neutrino-oscillation experiments 
using the NuMI beam, as well as Vf, studies in the MINOS experiment. 

2.2 Neutrino Scattering Physics 

A neutrino scattering experiment in the NuMI near experimental hall offers a unique opportunity to 
study a broad spectrum of physics topics with measurement precision heretofore unachievable. Several 

of these topics have not yet been studied in any systematic, dedicated way. For other topics, the few 
results that do exist are compromised by large statistical and systematic errors. Topics particularly open 
to rapid progress upon exposure of MINERz/A in the NuMI beam include: 

• Precision measurement of the quasi-elastic neutrino-nucleus cross-section, including its and 

dependence, and study of the nucleon axial form factors. 

• Determination of single- and double-pion production cross-sections in the resonance produc- 
tion region for both neutral-current and charged-current interactions, including a study of isospin 
ampUtudes, measurement of pion angular distributions, isolation of dominant form factors, and 
measurement of the effective axial- vector mass. 

• Clarification of the W (= mass of the hadronic system) transition region wherein resonance pro- 
duction merges with deep-inelastic scattering, including tests of phenomenological characteriza- 
tions of this transition such as quark/hadron duality. 

• Precision measurement of coherent single-pion production cross-sections, with particular atten- 
tion to target A dependence. Coherent -k^ production, especially via neutral-currents, is a signifi- 
cant background for next-generation neutrino oscillation experiments seeking to observe v^j^ Ug 
oscillation. 

• Examination of nuclear effects in neutrino-induced interactions including energy loss and final- 
state modifications in heavy nuclei. With sufficient V running, a study of quark flavor-dependent 
nuclear effects can also be performed. 

• Clarification of the role of nuclear effects as they influence the determination of sin^ 9^ via 
measurement of the ratio of neuttal-current to charged-current cross-sections off different nuclei. 
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• With sufficient v running, much-improved measurement of the parton distribution functions will 
be possible using a measurement of all six v and V structure functions. 

• Examination of the leading exponential contributions of perturbative QCD. 

• Precision measurement of exclusive strange-particle production channels near threshold, thereby 
improving knowledge of backgrounds in nucleon-decay searches, determination of Vus, and en- 
abling searches for strangeness-changing neutral-currents and candidate pentaquark resonances. 
Measurement of hyperon-production cross-sections, including hyperon polarization, is feasible 
with exposure of MINERzvA to V beams. 

• Improved determination of the effective charm-quark mass (mc) near threshold, and new mea- 
surements of Vf-d, s{x) and, independently, s{x). 

• Studies of nuclear physics for which neutrino reactions provide information complementary to 
JLab studies in the same kinematic range. 

In addition to being significant fields of study in their own right, improved knowledge of many 
of these topics is essential to minimizing systematic uncertainties in neutrino-oscillation experi- 
ments. 

2.2.1 Low-energy neutrino cross-sections 

This is a topic of considerable importance to both present and proposed future (off-axis) neutrino os- 
cillation experiments. Available measurements of both total and exclusive cross-sections from early 
experiments at ANL, BNL, CERN and FNAL all have considerable uncertainties due to low statistics 
and large systematic errors, including poor knowledge of the incoming neutrino flux|El. A working 
group|3 1 to assemble all available data on u and V cross-sections and to determine quantitative require- 
ments for new experiments has been established by members of this collaboration. MINERi/A will be 
able to measure these cross-sections with negligible statistical errors and with the well-controlled beam 
systematic errors needed for the MINOS experiment 

2.2.2 Quasi-elastic scattering 

Charged-current quasi-elastic reactions play a crucial role in both non-accelerator and accelerator neu- 
trino oscillation studies, and cross-section uncertainties - often expressed as uncertainty in the value of 
the axial-vector mass - are a significant component in error budgets of these experiments. There have 
been recent advances in the measurement of the vector component of elastic scattering from SLAC and 
Jefferson Lab. Measurement of the neutrino quasi-elastic channel is the most direct way to improve 
our knowledge of the axial-vector component to this channel. MINERi/A's ability to carefully measure 
da I dQ"^ to high allows investigation of the non-dipole component of the axial- vector form factor to 
an unprecedented accuracy. Combining these MINERz/A measurements with present and future Jeffer- 
son Lab data will permit precision extraction of all form factors needed to improve and test models of 
the nucleon. In addition, due to the well-constrained kinematics of this channel, a careful study of the 
muon and proton momentum vectors allows an important probe of nuclear effects. 
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2.2.3 Resonances and transition to deep-inelastic scattering 

Existing data on neutrino resonance-production is insufficient for the task of specifying the complex 
overlapping A and N* resonance amplitudes and related form-factors which characterize the 1-5 GeV 

regime. Neutrino Monte-Carlo programs trying to simulate this kinematic region have used early 
theoretical predictions by Rein and Sehgal|4| or results from electro-production experiments. Recently 
Lee and Sato|5| have developed a new model for weak production of the A resonance. Paschos and 
collaborators! 6 I have also contributed to this effort. It is noteworthy that the theoretical and exper- 
imental picture of the resonance region is far more obscure than the quasi-elastic and deep-inelastic 
scattering (DIS) regions which border it and that much of the relevant MINOS event sample falls inside 
this poorly-understood resonance region. 

Recent work at Jefferson Lab|7| shows strong support for quark/hadron duality, which relates the 
average resonance production cross-section to the DIS F2 structure function. How to incorporate this 
new paradigm into neutrino Monte Carlos is currently being studied. An analysis by Bodek and YangjHJ 
offers a very promising procedure for fitting F2 in the low-Q^/high-x region. Extrapolating their results 
through the resonance region yields values of F2 consistent with duality arguments and the Jefferson Lab 
results mentioned above. The resonance and transition region will be carefully examined by MINERz^A. 



2.2.4 Coherent pion production 

Both charged- and neutral-current coherent production of pions result in a single forward-going pion 
with little energy transfer to the target nucleus. In the neutral-current case, the single forward-going vr'^ 
can mimic an electron and be misinterpreted as a Uf. event. Existing cross-section measurements for this 
reaction are only accurate to ~ 35% and are only available for a limited number of target nuclei. 



2.2.5 Studying nuclear effects with neutrinos 

The study of nuclear effects with neutrinos can be broadly divided into two areas. The first area involves 
the kinematics of the initial interaction (spectral function of the struck nucleon within the nucleus and 
Pauli-excluded interactions) and the evolution of the hadronic cascade as it proceeds through the nu- 
cleus. This aspect has direct and important application to the MINOS neutrino oscillation experiment 
since it can drastically distort the initial neutrino energy by mixing final states to such an extent that the 
visible energy observed in the detector is much different than the initial energy. 

The second area involves modification of the structure functions, Fj and, consequently, the cross- 
section of v-A scattering compared to z^-nucleon scattering. Nuclear effects in DIS have been stud- 
ied extensively using muon and electron beams, but only superficially for neutrinos (in low-statistics 
bubble-chamber experiments). High-statistics neutrino experiments have, to date, only been possible 
using heavy nuclear targets such as iron-dominated target-calorimeters. For these experiments, results 
from e/;U-A analyses have been applied to the results. However, there are strong indications that the 
nuclear corrections for dfi-A and u-A are different. Among these differences is growing evidence for 
quark-flavor dependent nuclear effects. A neutrino-scattering program at NuMI would permit a sys- 
tematic, precision study of these effects, by using a variety of heavy nuclear targets and both u and V 
beams. 
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2.2.6 Strangeness and charm production 

MINERz/A will allow precise measurement of cross-sections for exclusive-channel strangeness associated- 
production (AS = 0) and Cabbibo-suppressed (AS' = 1) reactions. Detailed studies of the hadronic 
systems will be carried out, including dependence, resonant structure, and polarization states for 
produced lambda hyperons. A detailed study of coupling strengths and form-factors characterizing the 
AS weak hadronic current is envisaged, which will hopefully reawaken efforts at detailed modelling of 
these reactions ||9l. MINERi^A observations of strangeness production near threshold will have ramifi- 
cations in other areas of particle physics, as for example with estimation of atmospheric-neutrino AS* 
backgrounds for nucleon decay searches with megaton-year exposure. Searches for new resonant states 
and new physics will of course be possible: we envisage a dedicated search for strangeness-changing 
neutral-current reactions and investigation of unusual baryon resonances such as the recently reported 
candidate pentaquark state (in K^n and K^p systems). Clean measurement of Vus should be feasible; 
it may be possible to address long-standing discrepancies between theory and experiment concerning 
hyperon beta-decay by exploring the related inverse reactions obtained via AS = 1 single-hyperon pro- 
duction by antineutrinos lllOI . The production of hyperons by neutrinos and antineutrinos would provide 
new information in the form of hyperon polarization which would reduce ambiguities which currently 
compromise the analysis of hyperon beta-decay processes. 

Although the neutrino energy spectrum is relatively low for a high-statistics charm study, it does 
cover the important threshold region where production rates are highly dependent on the mass of the 
charm quark. Depending on the value of nic, the expected number of charm events could change by as 
much as 50% in this sensitive region. 

2.2.7 Extracting parton distribution functions 

Neutrinos have long been a particularly sensitive probe of nucleon structure. One obvious reason is the 
neutrino's ability to directly resolve the flavor of the nucleon's constituents: u interacts with d, s, u and 
c while the V interacts with u, c, d and s. This unique ability of the neutrino to "taste" only particular 
flavors of quarks enhances any study of parton distribution functions. Study of the partonic structure 
of the nucleon, using the neutrino's weak probe, would complement the on-going study of this subject 
with electromagnetic probes at Jefferson Lab as well as earlier studies at SLAC, CERN and FNAL. 
With the high statistics foreseen for MINERz^A, as well as the special attention to minimizing neutrino 
beam systematics, it should be possible for the first time to determine the separate structure functions 
2xFf^(x,Q2), 2xFf^(x,g2), F^^{x,Q^), F^^{x,Q^), xF!(^{x,Q^) and xF^^{x,Q^). This in 
turn would allow much-improved knowledge of the individual sea-quark distributions. 

2.3 The MINERz/A Detector 

To perform the full spectrum of physics outlined in this proposal, the MINERi^A target/detector must 
be able to: 

• Identify muons and measure their momentum with high precision, 

• identify individual hadrons and vr*^ and measure their momentum, 

• measure the energy of both hadronic and electromagnetic showers with reasonable precision. 
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• minimize confusion of neutral-current and charged-current event classifications, and 

• accommodate otiier nuclear targets. 

These goals can be met by a relatively compact and active target/detector consisting of a central sec- 
tion of essentially solid scintillator bars (Figure O. This central detector is surrounded on all sides by 
an electromagnetic calorimeter, a hadronic calorimeter and a magnetized muon-identifier/spectrometer. 
The detector has the approximate overall shape of a hexagon (to permit three stereo views) with a cross- 
section of 3.55 m minor and 4.10 m major axis. The length is up to 5.9 m depending on how close 
MINERzvA can be placed to the MINOS near detector. The active plastic scintillator volume is 6.1 
tons allowing variable sized fiducial volumes depending on the channel being studied. At the upstream 
end of the detector are nuclear targets consisting of 1 ton of Fe and Pb. Significant granularity and 
vertex-reconstruction accuracy can be achieved by the use of triangular-shaped extruded plastic scin- 
tillator(CH) bars with 3.3 cm base, 1.7 cm height and length up to 4.0 m, with an optical fiber placed 
in a groove at the base of the bar for readout. A second triangular shape with base 1.65 cm and height 
1.7 cm (1/2 of the larger triangles) will be used in the barrel and downstream calorimeter detectors. 
Recent work at the Fermilab Scintillator R&D Facility has shown that using light division across trian- 
gularly shaped scintillator strips of this size can yield coordinate resolutions of a few millimeters. The 
orientation of the scintillator strips are alternated so that efficient pattern recognition and tracking can 
be performed. 

Following the downstream end of the central detector are electromagnetic and hadronic calorime- 
ters. MINERz^A should be placed as close as possible to the upstream face of the MINOS near detector 

Muon Ranger (MR) Outer Detector (OD) Veto 








Figure 1 : A schematic side view of the MINERz^A detector with sub-detectors labeled. The neutrino 
beam enters from the right. 
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in order to use that detector's magnetic field and steel as an external muon-identifier and spectrometer 
for the forward-going muons, and as a calorimeter for any hadronic energy exiting MlNERz/A itself. 
Moving the MlNERz^A detector further upstream from the MINOS detector will decrease the accep- 
tance for muons in the MINOS detector. If necessary a "muon ranger", consisting of 1 .2 m of segmented 
and magnetised iron, will be added to help identify and measure the momentum of low-energy muons. 
For high-energy muons, the MINOS near detector will provide much better momentum resolution than 
the muon ranger 

With this design, even at the lowest beam-energy setting, MlNERz/A will collect more than 580 K 
events per 2.5 x lO^'^ POT in a 3 ton active target fiducial volume and just under 200 K events in each 
of the nuclear targets. 

The statistics from a several-year MINOS run will suffice to study all the physics topics listed 
above, although some measurements would be limited in kinematic reach by the beam energies used for 
MINOS. In addition, all studies involving V channels would be somewhat limited with the currently- 
planned MINOS exposure, which includes relatively Uttle V running. 
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3 Low-Energy Neutrino Scattering Overview 



3.1 Form Factors and Structure Functions 

Several formalisms are used to discuss electron-nucleon and neutrino-nucleon scattering, and the cor- 
responding reactions on nuclear targets. 

Inclusive lepton scattering can be described in the language of structure functions or in terms of 
form factors for the production of resonant final states. The two descriptions are equivalent and there 
are expressions relating form-factors to structure functions. In electron scattering, the vector form 
factors can be related to the two structure functions Wi and W2 (which are different for neutrons and 
protons), or equivalently F2 and R. 

In neutrino scattering, there are three structure functions Wi, W2 and W3 (or F2, R and xFs), 
different for neutrons and protons, and containing both vector and axial-vector components. There are 
also two other structure functions (important only at very low energies) whose contributions depend on 
the final-state lepton mass; these can be related to the dominant structure functions within the framework 
of theoretical models. 

3.2 Electron versus Neutrino Scattering 

From the conservation of the vector current (CVC), the vector structure functions (or form-factors) mea- 
sured in electron scattering can be related to their counterparts in neutrino scattering for specific isospin 
final states. For elastic scattering from spin-i quarks or nucleons, these relationships between vector 
form factors are simple. For production of higher spin resonances, the relations are more complicated 
and involve Clebsch-Gordon coefficients. 

In contrast, the axial structure functions in neutrino scattering cannot be related to those from elec- 
tron scattering, except in certain limiting cases (for example, within the quark/parton model at high 
energies with V=A). At low Q^, the axial and vector form factors are different, e.g. because of the 
different interactions with the pion cloud around the nucleon. 

Another difference arises from nuclear effects in inclusive neutrino vs. electron scattering. Nuclear 
effects on the axial and vector components of the cross-section can differ due to shadowing, and can 
also affect valence and sea quarks differently. 

3.3 Sum Rules and Constraints 

Several theoretical constraints and sum rules can be tested in electron and neutrino reactions (or applied 
in the analysis of data). Some of the sum rules and constraints are vaUd at all values of Q^, and some 
are valid only in certain limits. 

The Adler sum rules apply separately to the axial and vector parts of Wi, W2, and W3 and are valid 
for all values of (since they are based on current algebra considerations). At high Q"^, these sum 
rules are equivalent to the statement that the number of u valence quarks in the proton minus the number 
of d valence quarks is equal to 1 . 

Other sum rules, such as the momentum sum rule (sum of the momentum carried by quarks and 
gluons is 1) and the Gross/Llewelyn-Smith sum rule (number of valence quarks is equal to 3), have 
QCD corrections and break down at very low Q^. 

As 0, the vector structure functions are further constrained by the measured photoproduction 

cross-section. Conversely, as — > cjo it is expected that the structure functions are described by QCD 



13 



and satisfy QCD sum rules. 
3.4 Final States 

Quasi-elastic' reactions, resonance production, and deep-inelastic scattering are all important compo- 
nents of neutrino scattering at low energies. 

To describe specific final states, one can use the language of structure functions, combined with 
fragmentation functions, at high values of Q^. At low values of Q^, many experiments describe the 
cross-sections for specific exclusive final states. Both of these pictures need to be modified when the 
scattering takes place on a complex nucleus. 

Figure |2l shows the total neutrino and anti-neutrino cross-sections (per nucleon for an isoscalar 
target) versus energy (at low energies) compared to the sum of quasi-elastic, resonant, and inelastic 
contributions. These two figures also show the various contributions to the neutrino and anti-neutrino 
total cross-sections that will be investigated in this experiment. 



'We should clarify that the neutrino community uses the term 'quasi-elastic' to describe a charged-current process in 
which a neutrino interacts with a nucleon to produce a charged lepton in the final state. The nucleon can be a free nucleon or 
a nucleon bound in the nucleus. The term 'quasi-elastic' refers to the fact that the initial state neutrino changes into a different 
lepton, and there is a single recoil nucleon in the final state (which changes its charge state). In contrast, the electron scattering 
community refers to electron-nucleon scattering with a single recoil nucleon as 'elastic' scattering. The term 'quasi-elastic' 
scattering is used by the electron scattering community to describe elastic electron-nucleon scattering from bound nucleons in 
a nucleus. Here the term 'quasi-elastic' refers to the fact that the bound nucleon is quasi-free. Both nomenclatures are used in 
the literature. 
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Figure 2: Total neutrino (top) and anti-neutrino (bottom) cross-sections divided by energy versus energy 
compared to the sum of quasi-elastic, resonant, and inelastic contributions from the NUANCE model. 
The sum is constructed to be continuous in W (= mass of the hadronic system) as follows. For W > 
2 GeV the Bodek-Yang model is used. The Rein-Sehgal model is used for M/^ < 2 GeV. In addition, 
a fraction of the Bodek-Yang cross-section is added to the Rein-Sehgal cross-section between W = 
1.7 GeV and W = 2 GeV. The fraction increases linearly with W from to 0.38 between W = 1.7 
and = 2 GeV. 
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4 Existing Neutrino Scattering Data 



Neutrino experiments dating back to the 1960's have played an important role in particle physics, in- 
cluding discovery of neutral currents, electroweak measurements, determination of the flavor compo- 
sition of the nucleon, measurements of the weak hadronic current, and QCD studies based on scaling 
violations in structure and fragmentation functions. 

In the 1-10 GeV energy range of interest to the current and future generation of neutrino oscillation 
studies, relevant data comes from bubble-chamber experiments that ran from the 1960's through the 
1980's. Gargamelle, the 12-foot bubble chamber at the Argonne ZGS, the 7-foot bubble chamber at 
the AGS at Brookhaven, the Big European Bubble Chamber (BEBC) at CERN, the Serpukhov bubble 
chamber SKAT, and the FNAL 15-foot bubble chamber studied neutrino and anti-neutrino interactions 
off free nucleons and heavy liquid targets. Spark-chamber and emulsion experiments from this era 
played a less prominent role but did make crucial measurements in a number of areas. 

Despite limited statistics, the excellent imaging capabilities of bubble chambers made a wide range 
of physics topics accessible. It is primarily this data that is used to tune our Monte Carlo simulations and 
provides the basis for our present understanding of low-energy neutrino cross-sections. While adequate 
for validating the models at some level, most of the bubble-chamber data-sets are limited in size and 
do not cover the full range of neutrino energy, nuclear targets and neutrino species {vlv) required for a 
complete understanding of neutrino interactions. Some of the main topics of interest for experiments of 
this era are described below. For each topic, an approximate count of the number of SPIRES publication 
references is included. 

4.1 Quasi-elastic Scattering 

(8 pubs) Studies of quasi-elastic charged-current (CC) interactions were among the first results from 
bubble-chamber neutrino exposures, and are the primary tool for studying the axial component of the 
weak nucleon current. While data were taken on both light (H2/D2) and heavy (Neon/propane/freon) 
targets, no attempts were made to extract measurements related to the nuclear system. Rather, the nu- 
clear system was treated as a complication requiring corrections. In many instances even this correction 
was not done, and the published data are for interactions on nucleons in a particular nucleus. This helps 
account for the large spread in data points between different experiments in Figure |9l 

4.2 Other Exclusive Charged-current Channels 

(19 pubs) Total cross-section measurements and studies of differential distributions were made for both 
light and heavy targets in each of the three charged-current single-pion channels. In nearly all cases cuts 
were placed on the hadronic invariant mass {e.g. < 1.4 GeV/c^) to limit the analysis to the resonant 
region. The results are shown in Figure |4l Fewer experiments published cross-sections for two- and 
three-pion channels. 

4.3 Neutral-current Measurements 

(22 pubs) Neutral-current (NC) measurements fall into three categories: elastic measurements in dedi- 
cated experiments, single-pion exclusive final-state measurements, or NC/CC ratio measurements in the 
deep-inelastic scattering (DIS) regime. NC/CC ratio measurements were made at high energies and ap- 
plied cuts on the energy transfer v to isolate the DIS regime. Single-pion studies of the NC electroweak 
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couplings and the isospin characteristics of the hadronic current in the resonance region suffered from 
lack of statistics. Tabled summarizes the published data. These processes are of particular interest, as 
they constitute one of the primary backgrounds to appearance in oscillation experiments. 

4.4 Hadronic Final States 

(32 pubs) A number of publications were devoted to inclusive measurements of the hadronic system 
produced in neutrino interactions. Multiplicity measurements, transverse momentum distributions, in- 
clusive particle production, fragmentation functions, and evaluation of the universality of hadron dy- 
namics were studied. In this area, hadronic mass cuts {e.g. W > 2 GeV) were applied to limit the 
analysis to the DIS region. 

4.5 Strange and Charmed Particle Production 

(27 pubs) Because of their clear signatures in photographic quality bubble chambers, exclusive and 
inclusive measurements of strange and charm particle production were popular topics. A survey of 
these results is given in Section 1931 

4.6 Total Cross-sections 

(19 pubs) Total charged-current cross-section measurements were a staple of bubble-chamber experi- 
ments. Their data is shown in Figure |3l The large errors are due to a combination of low statistics and 
poor knowledge of the v beam. 

4.7 Structure Functions 

(18 pubs) Numerous experiments, particulary those at higher energies (and of course all the large calori- 
metric neutrino detectors like CDHS, CCFR, NuTeV, etc that followed) measured structure functions. 
Neutrino experiments are complementary to studies with electron and muon beams as they allow ex- 
traction of the valence quark distributions through measurement of xF^ as well as independent analysis 
of the strange quark content via di-muon production. These experiments made possible precision elec- 
troweak and QCD measurements with the NC/CC ratio and scaling violation in the structure functions. 

4.8 Summary 

Viewed from a historical perspective, the results from these experiments clearly reflected the topics of 
interest (and the theoretical tools available) at the time they were performed, and some general trends 
are clear. These experiments focused on two regimes. First, low scattering: the non-perturbative 
regime where the scattering takes place from a single nucleon. By measuring total and differential cross- 
sections for exclusive channels (like quasi-elastic and A production), these experiments studied in detail 
the weak hadronic current of the nucleon. Parton-model studies form a second, complementary class of 
experiments, studying scaling phemonena like total cross-sections, structure functions, scaling-variable 
distributions, and inclusive final-state dynamics, and applying kinematic cuts to remove resonant and 
quasi-elastic contributions. 

This dichotomy reflects the fact that decent models only existed for the extreme perturbative and 
non-perturbative limiting cases. The resonant/DIS transition region, where perturbative QCD breaks 
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Table 1 : Neutral-current measurements 



down, was avoided because a clear theoretical framework for it was not available. With the current 
generation of duahty studies at JLab and elsewhere, this complex but fundamental region is just now 
being fruitfully probed. 

Another area of difficulty was treatment of nuclear effects. While heavy targets gave bubble cham- 
bers increased target mass, the confounding effects of the nuclear environment on the target kinematics 
and observed final states were a topic which was largely ignored. Very few nuclear physics studies 
were ever carried out with neutrinos, and these had only the most naive models available for compar- 
ison. These studies focused on nuclear rescattering of produced pions, shadowing and EMC effects, 
formation-zone studies, and inclusive production of slow particles. Neither the small samples nor the 
models available allowed neutrinos to probe the nuclear environment in detail. 

These "holes" in existing neutrino data and related phenomenology are now becoming increasingly 
evident. The MINOS experiment, for instance, will see a wide-band beam of 1-50 GeV neutrinos. 
Since a significant fraction of the interactions in MINOS are in the "transition" region, and nearly all 
take place on Iron nuclei, the areas of study neglected during the bubble-chamber era begin to loom 
large. MINOS, and the neutrino-oscillation experiments that will follow it, will be forced to confront 
them to achieve maximum sensitivity. Part |n] of this proposal explains in detail how MINERi^A will 
not only address fundamental topics in nuclear and neutrino physics, which are compelling in their own 
right, but also substantially improve the quality of results from future oscillation experiments. 
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Figure 3: The NEUGEN prediction for the charged-current cross-section {a jEy) from an isoscalar 
target compared with data from a number of experiments. Quasi-elastic and resonance contributions 
are also shown. Data are from: CCFRR Ol, CDHSW QS, GGM-SPS d, BEBC (ID, ITEP (TtI, 
SKAT iJii, CRS EH, ANL |20|, BNL Ell, GGM-PS El, ANL-QEL El, BNL-QEL 
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Charged Current Single Pion Production 




Figure 4: Cross-sections for chai-ged-current single-pion production. Plot A: Uf^+p /i +p+Tr^, Plot 
B: v^ + n fi^ + n + 7r+, Plot C: Ufj^ + n —>■ fi~ +p + Tr^. Solid lines are the NEUGEN predictions for 
W<1.4 GeV (plot A) and W<2.0 GeV (plots B and C). The dashed curve is the NEUGEN prediction 
with no invariant mass cut, for comparison with the BNL data. Data are from: ANL I35II20I . BNL 1361 . 
FNAL[37|,BEBCk38| 
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5 The NuMI Beam and MINERi/A Event Sample 



The NuMI neutrino beam is produced from vr- and K-decay in a 675 m decay pipe beginning 50 m 
downstream of a double horn focusing system. At the end of the decay pipe a 10 m long hadron 
absorber stops the undecayed secondaries and non-interacting primary protons. Just downstream of the 
absorber, 240 m Dolomite is used to range out muons before the v beam enters the near detector hall. 
Figure 121 shows the beam component and near detector hall layout. 




Figure 5: Layout of NuMI beamline components and near detector hall. 



5.1 Energy Options 



The neutrino energy distribution of the NuMI beam can be chosen by changing the distance of the target 
and second horn with respect to the first horn, as in a zoom lens. These three configurations result in 
three beam energy tunes for the low (LE), medium (ME), and high (HE) energy ranges respectively. 
However, to switch from one the beam mode into an alternate configuration will require down time to 
reconfigure the target hall and a loss of beam time. An alternative to this which allows the peak energy 
to be varied is to change the distance of target from the first horn and leave the second horn fixed in the 
LE position. This can be accomplished remotely with maximum transit of -2.5 m motion of the target 
upstream of the first horn from its nominal low energy position. The configurations corresponding to 
target -1.0 m from nominal results in a "semi-medium" energy beam tune (sME) and target -2.5 m 
from nominal will produce "semi-high" energy beam (sHE). These semi-beam configurations are less 
efficient and result in lower event rates than the ME and HE beams. A considerably more efficient 
sHE beam is possible with three-day downtime to allow the target to be moved back to its nominal 
HE position of -4.0 m. This more efficient sHE(-4.0) beam would yield over 50% more events than 
the sHE(-2.5) beam. For the MINOS experiment the beamline will be operating mainly in its lowest 
possible neutrino energy configuration to be able to reach desired low values of Am^. However, to 
minimize systematics, there will also be running in the sME and sHE configurations described above. 
The neutrino energy distributions for the LE, sME, and sHE running modes are shown in Figure |6l 
Figure shows the event energy distributions for the ME and HE beam configurations for comparison. 
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Figure 6: The neutrino charged-current event energy distribution for the three configurations of the 
NuMI beam corresponding to low-energy (LE), medium-energy (sME) and high-energy (sHE). 
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Figure 7: The neutrino charged-current event energy distribution for the high-rate medium and high- 
energy beam configurations (ME and HE) which involve movement of the second horn as well as target 
position. 
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5.2 MINERz/A Event Rates 



Table |2l shows the charged-current event rates per 10^° protons on target (PoT) per ton for the three 
beam configurations discussed above. In addition, the same configurations but with horn-current re- 
versed provide anti-neutrino beams. Event rates for charged-current events using anti-neutrino beam 
configurations (LErev, MErev, and HErev) are also shown along with their background components. 
Running in these modes would be highly desirable for MINERi/A physics. 



CC Events/10^'^ 


PoT/ton 


Beam 




CCVe 


LE 


78 K 


1.1 K 


sME 


158 K 


1.8 K 


sHE 


257 K 


2K 




CCV^ 


CCu^ 


LErev 


26 K 


34 K 


MErev 


56 K 


10 K 


HErev 


75 K 


13 K 



Table 2: MINERi/A event rates for different beam configurations. 
5.3 Baseline MINOS Run Plan 

Table |3l shows a scenario for predicted PoT over a conservative hypothetical four-year MINOS run. 
From this table the total integrated charged-current event samples for a four-year MINERzvA run would 
be 940 K charged-current events per ton and 275 K 17^ charged-current events per ton. 







Scenario for PoT per year (x 10 


20^ 




year 


total PoT LE 


sME 


sHE 


LErev 


MErev 


HErev 


2006 


3.0 


3.0 


0.0 


0.0 


0.0 


0.0 


0.0 


2007 


4.0 


3.0 


0.7 


0.3 


0.0 


0.0 


0.0 


2008 


4.0 


0.0 


0.0 


0.0 


2.5 


1.0 


0.5 


2009 


4.0 


1.0 


0.5 


0.5 


0.5 


0.5 


1.0 


Total 


15.0 


7.0 


1.2 


0.8 


3.0 


1.5 


1.5 



Table 3: Hypothetical proton luminosity scenario for a four-year run. 



5.4 MINERz/A Data Samples 

The event rates for physics processes of interest to MINERz/A for the four-year scenario discussed in 
the previous section are summarized in the Table |4] 

The distribution of the number of interactions expected for different XBj and values are shown 
for the quasi-elastic, resonant and deep-inelastic channels in Tables |5l|6l and respectively. The spread 
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of the quasi-elastic events in x is due to the smearing from the Fermi motion of the target nucleon. For 
clarity the xbj and distributions of the total and deep-inelastic event samples are shown in Figure[Sl 
These tables are based on the four-year scenario outlined in Table |3] 

The number of interactions expected during the full four-year exposure of the detector in the NuMI 
beam eclipses the number of events recorded in the bubble-chamber experiments described in Section|3] 
by several orders of magnitude. The implications of this unprecedented event sample for physics are 
described in later sections. It would not be an exaggeration to observe that this large sample of neutrino 
interactions will reduce many of the systematic errors currently limiting the sensitivity of neutrino 
oscillation experiments and allow detailed study of kinematic regions that are presently rather poorly 
understood. 

Were MINERz/A the prime user of NuMI, the beamline would be run in the high-energy config- 
uration with energies in the 5-25 GeV range. This configuration offers the ability to study neutrino 
interactions across an appreciable fraction of the xbj range at reasonable Q^. In HE beam mode ex- 
pected event rates would be 580 K charged-current events per 10^'^ PoT per ton, over twice as many 
as the sHE(-2.5) beam. 

5.5 Accuracy of Predicted Neutrino Flux 

As mentioned earlier, one of the significant advantages of MINERi/A over previous wide-band neutrino 
experiments is the expected accuracy with which the neutrino absolute and energy dependent flux is 
known. Since the NuMI beamline has been designed for the MINOS neutrino oscillation experiment, 
particular attention has been paid to control and knowledge of the beam of neutrinos being used in the 
experiment. 

The biggest uncertainty in the predicted energy spectrum of the neutrinos comes directly from the 
uncertainty of hadron prodution spectrum of the vr and K parents of the neutrinos. To help reduce this 
uncertainty, there is an approved Fermilab experiment E-907 1 011119311 which has as it's main goal the 
measurement of hadron production spectra off various nuclear targets. One of the measurements that 
will be made by E-907 is an exposure of the NuMI target to the 120 GeV Main Injector proton beam. 
By using the NuMI target material and shape, E-907 will be able to provide the spectra coming off 
the target including all of the secondary and tertiary interactions which can significantly modify the 
produced spectra. It is expected that with the input from E-907, the absolute and energy dependent 
shape of neutrinos per POT will be known to ^ 3%. 

For the absolute flux of neutrinos there is a second uncertainty which must be considered and that 
is the accuracy with which we know the number of protons on target. With the planned NuMI primary 
proton beamline instrumentation! 12 1, the number of protons on target will be known to between (1 - 
3)%, the range being determined by the calibration techniques used to control drift of the primary beam 
toroid devices. 

To summarize, the energy shape of the NuMI beam should be known to 3% while the absolute flux 
should be known to between (3 - 5)%. 



26 



Event Rates per ton 



Process 


CC 


NC 


Elastic 


103 K 


42 K 


Resonance 


196 K 


70 K 


Transition 


210 K 


65 K 


DIS 


420 K 


125 K 


Coherent 


8.3 K 


4.2 K 


Total 


940 K 


288 K 



Table 4: Total event rates for different reaction types, per ton, for the four-year scenario outlined in 
Table m 




10 15 20 
Q' (GeV/c) 



10 15 20 
Q' (GeV/c)' 



Figure 8: Kinematic distributions {xbj and Q^) expected for deep-inelastic (top left and bottom left) 
and all event types (top right and bottom right). 
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Quasi-elastic events per ton vs. {i^Bj,Q'^) for four-year scenario 



Q2 (GeV/c)2 



XBj 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-2.5 


2.5-3.0 


3.0-3.5 


3.5-4.0 


4.0+ 


Total 


0.0-0.1 


1068 


























1068 


0.1-0.2 


2372 


























2372 


0.2-0.3 


3366 


12 























3378 


0.3-0.4 


5552 


199 


49 




















5800 


0.4-0.5 


7812 


361 


49 




















8222 


0.5-0.6 


9974 


1217 


173 


37 


39 











12 


11452 


0.6-0.7 


11377 


3717 


956 


161 


198 


37 





12 


12 


16470 


0.7-0.8 


9663 


4682 


2061 


906 


427 


173 


136 


12 


62 


18122 


0.8-0.9 


8645 


4744 


2235 


1296 


608 


322 


161 


99 


248 


18358 


0.9-1.0 


6868 


4956 


2198 


1246 


509 


360 


223 


86 


248 


16694 


Total 


66697 


19888 


7721 


3646 


1781 


892 


520 


209 


582 


101936 



Table 5: Quasi-elastic interactions expected per ton for the four-year scenario of Table|3] 



Resonant events per ton vs. {xBj,Q'^) for four-year scenario 



Q2 (GeV/c)2 



XBj 


0.0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-2.5 


2.5-3.0 


3.0-3.5 


3.5-4.0 


4.0+ 


Total 


0.0-0.1 


48169 


49 























48219 


0.1-0.2 


46132 


7763 


173 




















54069 


0.2-0.3 


27649 


15104 


2881 


310 


12 














45958 


0.3-0.4 


16135 


15613 


6508 


1689 


298 


62 


12 


12 





40331 


0.4-0.5 


5974 


13576 


6359 


2943 


1416 


521 


86 


86 


24 


30990 


0.5-0.6 


1018 


8968 


5924 


3279 


1738 


1018 


496 


211 


236 


22892 


0.6-0.7 


74 


4012 


4571 


2956 


1577 


993 


434 


459 


571 


15650 


0.7-0.8 


12 


1217 


2583 


1788 


1217 


919 


558 


496 


770 


9564 


0.8-0.9 





260 


844 


745 


757 


732 


472 


347 


633 


4794 


0.9-1.0 





111 


347 


285 


397 


310 


124 


136 


534 


2248 


Total 


145163 


66673 


30190 


13995 


7412 


4555 


2182 


1747 


2768 


274715 



Table 6: Resonant interactions expected per ton for the four-year scenario of Table |3] 
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Deep-inelastic events per ton vs. ()^Bj,Q^) for four-year scenario 



g2 (GeV/c)2 



XBj 


0-2 


2-4 


4-6 


6-8 


8-10 


10-12 


12-14 


14-16 


16-20 


20+ 


Total 


0.0-0.1 


100276 


1987 


198 


24 




















102485 


0.1-0.2 


123988 


13688 


2670 


832 


310 


86 


12 





12 





141598 


0.2-0.3 


79632 


24954 


5738 


1676 


956 


360 


223 


111 


37 


86 


113773 


0.3-0.4 


39598 


23028 


8011 


2633 


1279 


521 


211 


186 


136 


211 


75814 


0.4-0.5 


15091 


15104 


5614 


2571 


1291 


658 


322 


248 


173 


322 


41394 


0.5-0.6 


4670 


7154 


3316 


1726 


894 


645 


161 


111 


99 


223 


18999 


0.6-0.7 


1366 


2620 


1664 


1043 


397 


236 


186 


86 


24 


136 


7740 


0.7-0.8 


472 


670 


509 


273 


223 


111 


149 


12 


37 


74 


2530 


0.8-0.9 


99 


173 


149 


24 


99 


37 


12 


12 





12 


617 


0.9-1.0 


74 


37 


37 





24 











12 





184 


Total 


365276 


89423 


27906 


10802 


5473 


2654 


1276 


766 


530 


1064 


505134 



Table 7: Deep-inelastic interactions expected per ton for the four-year scenario of Table |3] 
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Part II 

Physics Motivation and Goals 
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6 Quasi-Elastic Scattering 



6.1 Quasi-elastic Cross-sections 

Quasi-elastic scattering makes up the largest single component of the total v-N interaction rate in the 
threshold regime Ey < 2 GeV. Precision measurement of the cross-section for this reaction, includ- 
ing its energy dependence and variation with target nuclei, is essential to current and future neutrino- 
oscillation experiments. Figures |9l and QOl summarize current knowledge of neutrino and anti-neutrino 
quasi-elastic cross-sections. Among the results shown, there are typically 10-20% normalization uncer- 
tainties from knowledge of the fluxes. These plots show that existing measurements have large errors 
throughout the Ey range accessible to MINERi/A (Figure |9l upper plot), and especially in the threshold 
regime which is crucial to future oscillation experiments (Figurej^J lower plot). Figure shows these 
large uncertainties extend to anti-neutrino measurements as well. 

MINERz^A will measure these quasi-elastic cross-sections with samples exceeding earlier (mostly) 
bubble-chamber experiments by two orders of magnitude. MINERz^A will also perform the first preci- 
sion measurement of nucleon form-factors for > \ (GeV/c)^ using neutrinos. 

Consistent and up-to-date treatment of the vector and axial-vector form-factors which characterize 
the nucleon weak current is essential to a realistic cross-section calculation. MINERi/A collaborators 
have been active in this area for some time|40|. Recent parameterizations and fits published by Budd, 
Bodek and Arrington are hereafter referred to as "BBA-2003" results. The curves in Figuresl9landfT0lare 
based on BBA-2003 form-factors, with the axial form-factor mass parameter set to Ma = 1-00 GeV/c^. 
The solid curves are calculated without nuclear corrections, while the dashed curves include a Fermi gas 
model. The dotted curves are calculations for Carbon nuclei and include Fermi motion, Pauli blocking, 
and the effect of nuclear binding on the nucleon form-factors as modeled by Tsushima et a/.|59|. 

Figure ^2 shows predictions for the cross section assuming the BBA form factors, with the axial 
mass = 1.00. The number of events assumes a 3 ton fiducial volume, out of the 6 tons of completely 
active target. We assume the efficiency for detections is 100% and perfect resolution. We would take 
all events for which there is a recoil proton traversing at least one plane or for which the recoil proton is 
absorbed in the nucleus. Therefore, efficiency and acceptance would be 100% since we have a hermetic 
detector with a side muon/hadron absrober. The error in the quasi-elastic cross section is then the 
statistical error and the error from the subtraction of resonance events which, through nuclear effects 
and close-in scatters, simulate the quasi-elastic signature.. This background uncertainty is not shown. 
Not shown also is the 4% normalization flux error. 

Nuclear effects reduce the calculated cross-sections by >10%; this sensitivity to the details of nu- 
clear physics shows that an understanding of final-state nuclear effects is essential to interpretation of 
quasi-elastic neutrino data. As as fine-grained tracking calorimeter, MINERz^A is designed to facilitate 
systematic comparison of quasi-elastic scattering (and other exclusive channels) on a variety of nuclear 
targets, providing a vastly improved empirical foundation for theoretical models of these important 
effects. 

6.2 Form-factors in Quasi-elastic Scattering 

MINERz^A's large quasi-elastic samples will probe the response of the weak nucleon current with 
unprecedented accuracy. The underlying V-A structure of this current include vector and axial-vector 
form-factors whose response is approximately described by dipole forms. The essential formalism 
is given by|41J 
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Figure 9: Compilation of neutrino quasi-elastic cross-section data. The data have large errors and 
are only marginally consistent throughout the Ei^ range accessible to MINERi^A (upper plot), and 
particularly in the threshold region (lower plot). Representative calculations are shown using BBA- 
2003 form factors with Myi=1.00 GeV. The solid curve is without nuclear corrections, the dashed curve 
includes a Fermi gas model l56l . and the dotted curve includes Pauli blocking and nuclear binding. 
The data shown are from FNAL 1983 [47|, ANL 1977 El, BNL 1981 Gl, ANL 1973 |@9|, SKAT 
1990 ©1, GGM 1979 |51|, LSND 2002 i52|, Serpukov 1985 f53l, and GGM 1977 El. 
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Figure 10: Compilation of anti-neutrino quasi-elastic cross-section data. As in Figure |9l the data have 
large errors, and considerable scatter among the different experiments. Theoretical predictions without 
(solid curve) and including nuclear corrections (dashed, dotted curves) are shown for comparison. The 
data shown are from SKAT 1990 1^, GGM 1979 L55J, Serpukov 1985 L53J, and GGM 1977 154J. 




< PiP2)\Jt\n{pi) >-- 



U{p2) 



where q = — k^, = (/ip — 1) — /in> and M = {nip + m„)/2. Here, /ip and /i„ are the proton and 
neutron magnetic moments. It is assumed that second-class currents are absent, hence the scalar (Fy) 
and tensor (F^) form-factors do not appear. 
The form-factors Fy{q'^) and iFy{q^) are given by: 



F'viq' 



4M 



AM' 



According to the conserved vector current (CVC) hypothesis, G^(g^) and G\^{q'^) are directly related 
to form-factors determined from electron scattering G^{q'^), G^{q'^), G\j{q^), and G\i{q^): 



gW) = gW)-gW 

The axial {Fa) and pseudoscalar {Fp) form-factors are 

9A 



gUq') 



FA{q' 



1 



Mi 



2 ' 



Fp{q' 



Gl,{q')-GUQ')- 

2M^FA{q^) 
Ml-q^ ■ 



In the differential cross-section, Fp{q'^) is multiplied by {mi/My, consequently its contribution to 
muon neutrino interactions is very small, except below 0.2 GeV. In general, the axial form-factor F^(g^) 
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Figure 1 1 : Cross section for MINERz/A assuming a 3 ton fiducial volume, 4 year run, perfect resolution, 
100% detection efficiency, BBA-2003 form factors with Ma = 1-00 GeV, and the Fermi gas model. 

can only be extracted from quasi-elastic neutrino scattering; at low q^, however, its behavior can also 
be inferred from pion electroproduction data. 

Until recently, it has been universally assumed that the form-factors' dependence is described by 
the dipole approximation. For example, the vector form factors are normally expressed: 



1 



M^ = 0.71GeV^ 



As discussed below, the dipole parameterization is far from perfect, and MINERi^A will be able to 
measure deviations from this form. 



6.2.1 Vector form-factor discrepancy at high Q"^ 

Electron scattering experiments at SLAC and Jefferson Lab (JLab) have measured the proton and neu- 
tron electromagnetic (vector) form-factors with high precision. The vector form-factors can be deter- 
mined from electron scattering cross-sections using the standard Rosenbluth separation technique|43 1, 
which is sensitive to radiative corrections, or from polarization measurements using the newer polariza- 
tion transfer technique|45 |. Polarization measurements do not directly measure form-factors, but rather 
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the ratio GeIGm- These form-factors can be related to their counterparts in quasi-elastic neutrino scat- 
tering by the CVC hypothesis. Naturally, more accurate form-factors translate directly to improved 
calculations of neutrino quasi-elastic cross-sections. 

Recently, discrepancies in electron scattering measurements of some vector form-factors have ap- 
peared; study of quasi-elastic reactions in MINERi/A may help reveal the origin these discrepacies. 
Figure ^1 shows the BBA-2003 fits to /XpG^/G^^^. There appears to be a difference between the two 
different methods of measuring this ratio. The fit including only cross-section data {i.e. Rosenbluth 
separation) is roughly fiat in and is consistent with form-factor scaling. This is expected if the elec- 
tric charge and magnetization distributions in the proton are the same. However, the newer polarization 
transfer technique yields a much lower ratio at high Q^, and indicates a difference between the elec- 
tric charge and magnetization distributions. The polarization transfer technique is believed to be more 
reliable and less sensitive to radiative effects from two-photon corrections. 
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Figure 12: Ratio of to G\[ as extracted by Rosenbluth separation measurements (diamonds) and 
polarization measurements(crosses). The data are in clear disagreement at high Q^. 



If the electric charge and magnetization distributions of the proton are indeed different, a test of the 
axial form-factor's high-Q^ shape can provide important new input to help resolve differences between 
electron scattering measurements. As discussed below, MINERi/A will be able to accurately measure 
the high-Q^ behavior of Fa- 



6.2.2 Form-factor deviations from dipole behavior 

Electron scattering shows that dipole amplitudes provide only a first-order description of form-factor 
behavior at high Q^. Figure shows the deviation of G\j from dipolar dependence. In general, 
these deviations are different for each of the form factors. 

Figure^Jshows the ratio of the BR\-2003 neutrino and anti-neutrino quasi-elastic cross-sections to 
the prediction using dipole vector form-factors (with G^ = and Ma kept fixed). This plot shows that 
the importance accurately parameterizing the form-factors. In MINERz^A, it will be possible to study 
the dependence of the form-factors beyond the simple dipole dipole approximation which has been 
assumed by all previous neutrino experiments. 
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Figure 14: Ratio of the neutrino and anti-neutrino quasi-elastic cross-sections calculated with BBA- 
2003 form-factors to the simple dipole approximation with = 0. 
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6.3 Axial Form-factor of the Nucleon 



Electron scattering experiments continue to provide increasingly precise measurements of the nucleon 
vector form-factors. Neutrino scattering, however, remains the only practical route to comparable pre- 
cision for the axial form-factors, in particular Fa{Q'^)- The fall-off of the form-factor strength with in- 
creasing is traditionally parameterized using an effective axial- vector mass Ma- Its value is known 
to be ^ 1.00 GeV/c^ to an accuracy of perhaps 5%. This value agrees with the theoretically-corrected 
value from pion electroproduction[42 |, 1.014 it 0.016 GeV/c^. Uncertainty in the value of Ma con- 
tributes directly to uncertainty in the total quasi-elastic cross-section. 

The fractional contributions of Fa, G^^j,G^j,G^^, and to the Q"^ distribution for quasi-elastic 
neutrino and anti-neutrino scattering with the NuMI beam are shown in Figure The contributions 
are determined by comparing the BBA-2003 cross-sections with and without each of the form-factors 
included. MINERi/A will be the first systematic study of Fa, which accounts for roughly half of the 
quasi-elastic cross-section, over the entire range of shown in the figure. 
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Figure 15: Fractional contributions of G^^,Gjfj,G^, G^ and Fa to the distributions for quasi-elastic 
neutrino (top) and anti-neutrino (bottom) samples with the NuMI beam. Because of interference terms, 
the sum of the fractions does not necessarily add up to 100%. 



6.3.1 Vector form-factors and Ma 

Earlier neutrino measurements, mostly bubble-chamber experiments on Deuterium, extracted Ma using 
the best inputs and models available at the time. Changing these assumptions changes the extracted 
value of Ma- Hence, precision measurement of Ma requires starting with the best possible vector 
form-factors, coupling constants, and other parameters. 

Figure^Jshows the Q"^ distribution from the Baker et al-f24\ neutrino experiment compared to the 
dipole form-factor approximation with = and Ma = 1.100 GeV/c^. Also shown are BBA-2003 
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Figure 16: Comparison of distributions using two different sets of form-factors. The data are from 
Baker et al.\2A\. Tiie dotted curve uses dipole form-factors with = and Ma = 1-10 GeV/c^. 
The dashed curve uses more recent BBA-2003 form-factors and Ma = 1-05 GeV/c^. It is essential 
to use the best possible information on vector form-factors from electron scattering experiments when 
extracting the axial form-factor from neutrino data. 



predictions with Ma = 1.050 GeV/c^. Use of more accurate electromagnetic form-factors requires a 
different Ma value to describe the same distribution. Thus, with the same value of qa, adopting the 
dipole approximation (and = 0) instead of the BBA-2003 form-factors may lead to an error in Ma 
of 0.050 GeV/c2. 

6.3.2 Measurement of the axial form-factor in MINERi/A 

Current and future high-statistics neutrino experiments at low energies {e.g. K2K, MiniBooNE, J- 
PARCnu and MINERi/A) use an active nuclear target such as scintillator (mostly Carbon) or water 
(mostly Oxygen). The maximum values that can be achieved with incident neutrino energies of 
0.5, 1.0, 1.5 and 2 GeV are 0.5, 1.2, 2.1 and 2.9 (GeV/c)^ respectively. Since K2K, MiniBooNE and 
J-PARCnu energies are in the 0.7-1.0 GeV range, these experiments probe the low < \ (GeV/c)^ 
region where nuclear effects are large (see Figures |^ and l23l and where the free-nucleon axial form- 
factor is known rather well from neutrino data on Deuterium (see Figure fT6b . The low {Q^ < 
1 (GeV/c)^) MiniBooNE and K2K experiments have begun to investigate the various nuclear effects in 
Carbon and Oxygen. 

At higher Q^, as shown by the BBA-2003 fits, the dipole approximation for vector form-factors can 
be in error by a factor of two when > 2 (GeV/c)^. There is clearly no reason to assume the dipole 
form will be any better for the axial form-factor. As shown in FigurefT6l there is very little data for the 
axial form-factor in the high-Q^ region (where nuclear effects are smaller). Both the low-Q^ (Q^ < 
1 (GeV/c)2) and high-Q^ (q2 > 2 (GeV/c)^) regions are accessible in higher-energy experiments like 
MINERi^A, which can span the 2-8 GeV neutrino energy range. MINERi/A's measurement of the axial 
form-factor at high will be essential to a complete understanding of the vector and axial structure of 
the neutron and proton. 

Figure ^\ shows the extracted values and errors of Fa in bins of from a sample of simulated 
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Figure 17: Estimation of Fa from a sample of Monte Carlo neutrino quasi-elastic events recorded in the 
MlNERz/A active Carbon target. Here, a pure dipole form for Fa is assumed, with Ma = 1 GeV/c^. 
The simulated sample and error bars correspond to four years of NuMl running. 

quasi-elastic interactions in the MlNERi/A active Carbon target, for a four-year exposure in the NuMl 
beam. Clearly the high-Q^ regime, which is inaccessible to K2K, MinibooNE and J-PARCnu, will be 
well-resolved in MlNERi/A. Figure^Jand ^]show these results as a ratio of FA/FA(Dipo\e), demon- 
strating MlNERz^A's ability to distinguish between different models of Fa- MlNERi^A will be able to 
measure the axial nucleon form-factor with precision comparable to vector form-factor measurements 
at JLab. 

The plot for FA{q^) is done by writing out the cross section as a quadratic function of FA{q^)- The 
coefficents of FAiq^) are functions of Ey and q^. The constant coefficient is a function of the measured 
(or predicted for MlNERi^A) cross section in a bin, as well as being a function of E^ and q^. The 
coefficients are integrated in energy over the flux and the q^ region of the bin. The q^ of the bin is 
determined by bin centering wrt da/dq^. For the extraction of FA{q^) for MlNERz^A we assume a 
dipole for FA{q^) with a value of Ma given by pion electro-production Ma= 1.014, BBA form factors, 
100% detection efficency and perfect resolution. 

For the extraction of data from Miller, Baker, and Kitagaki, we use the BBA form factors to de- 
termine the coefficients. Their plots of da/dq^ is used to extract FA{q^)- The overall nomalization 
of their data is not given. (In addition they probably do not know their normalization to better than 
10%). Hence, we set the overall normalization of FA{q^) by getting the overall normalization from the 
cross section using a dipole FA{q^) with A/^=1.014. Since MlNERz/Awill measure a normalized cross 
section, MlNERi^Awill not need to determine the overall normalization of FA{q^) by assuming a form 
forFA(g2). 

We recognize that at low Q^, there are large nuclear corrections and the recoil proton cannot be 
well measured. This is why for the measurement of F^Q'^) we plan on integrating the axial form factor 
from Q^=0 to Q^=0.2 to 0.4 GeV. Since the axial form factor is known at Q^=0 from neutron decay to 
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Figure 18: Extracted ratio F4/Fyi(DipoIe) from a sample of Monte Carlo quasi-elastic interactions 
recorded in the MINERz^A active Carbon target, from a four-year exposure in the NuMI beam. The 
MINERz/A points assume this ratio is described by the ratio of G£;(Polarization)/G£;(dipole). Fa is 
extracted from deuterium bubble chamber experiments using the da/dq^ from the papers of FNAL 
1983 Ell BNL 1981 EH, and ANL 1982 ^ 



be 1.26 with high precision, this is not an issue. We show two different models for Fa as a function 
of (which differ by a factor of 5 at high Q^, as indicated by G^/G^j data. At high Q^, the proton 
track is long and ID'd in the scintillator, and very high protons also stop in the side or downstream 
absorber, since we have hermetic detector, the acceptance is 100%. Here again, the background from 
MISId' delta's needs to be considered. Or alternatively, tighter de/dx cuts on the proton track (thus 
trading efficiencies for pion rejection). 

Figure |20l shows a typical quasi-elastic event, as simulated in MINERi/A. 



6.4 Nuclear Effects in Quasi-elastic Scattering 
6.4.1 Fermi gas model 

There are three important nuclear effects in quasi-elastic scattering from nuclear targets: Fermi motion, 
Pauli blocking, and corrections to the nucleon form factors due to distortion of the nucleon's size and 
its pion cloud in the nucleus. Figure|^shows the nuclear suppression versus from aNUANCE|58| 
calculation 156 1 using the Smith and Moniz ll57l Fermi gas model for Carbon. This nuclear model in- 
cludes Pauli blocking and Fermi motion but not final state interactions. The Fermi gas model uses a 
nuclear binding energy e = 25 MeV and Fermi momentum kj = 220 MeV/c. Figure |22lfrom Moniz et. 
al.ll57 l shows how the effective kf and nuclear potential binding energy e (within a Fermi-gas model) 
for various nuclei is determined from electron scattering data. The effective A;^ is extracted from the 
width of the scattered electron energy distribution, and the binding energy e from the shifted location of 
the quasi-elastic peak. 
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Figure 19: Extracted ratio F/i/F/i(Dipole) from a sample of Monte Carlo quasi-elastic interactions 
recorded in the MINERi/A active Carbon target, from a four-year exposure in the NuMI beam. The 
MINERz/A points assume this ratio is described by the ratio of G£;(Cross-Section)/G^(dipole), which 
was the accepted result for G^, before new polarization transfer measurements. The extracted values of 
Fa for the deuterium bubble chamber experiements are the same as the previous figure 



6.4.2 Bound nucleon form-factors 

The predicted distortions of nucleon form-factors due to nuclear binding are shown in Figure |23] as the 
ratios of Fi, F2, and Fa for bound and free nucleons. With a variety of nuclear targets, MINERz^A will 
be able to compare measured form-factors for a range of light to heavy nuclei. 



6.4.3 Intra-nuclear rescattering 

In neutrino experiments, detection of the recoil nucleon helps distinguish quasi-elastic scattering from 
inelastic reactions. Knowledge of the probability for outgoing protons to reinteract with the target 
remnant is therefore highly desirable. Similarly, quasi-elastic scattering with nucleons in the high- 
momentum tail of the nuclear spectral function needs to be understood. More sophisticated treatments 
than the simple Fermi gas model are required. Conversely, inelastic reactions may be misidentified 
as quasi-elastic if a final-state pion is absorbed in the nucleus. Because of its constrained kinematics, 
low-energy neutrino-oscillation experiments use the quasi-elastic channel to measure the (oscillated) 
neutrino energy spectrum at the far detector; the uncertainty in estimation of this background due to 
proton intra-nuclear rescattering is currently an important source of systematic error in the K2K exper- 
iment. 

The best way to study these effects is to analyze electron scattering on nuclear targets (including the 
hadronic final states) and test the effects of the experimental cuts on the final-state nucleons. MINERz^A 
can address proton intra-nuclear rescattering by comparing nuclear binding effects in neutrino scattering 
on Carbon to electron data in similar kinematic regions. Indeed, MINERz^A members will be working 
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MINERVA SIDE VIEW Run Event 20lntTypeQE 

CC/NC 1 Mech. nu-n 
Vertex ( 0.0. 0.0,1336.3) 
PNEU 1 4 ( 0.0000,0.0000.2.2266,2.2266 ) 
PLEP 13 ( -,6252,-.47S6,1,5a04J.7691 ) 
X 0.952 q2 0.817 y 0.205 w2 0.9383 




Figure 20: A simulated charged-current quasi-elastic interaction in MINERz^A. The proton (upper) and 
muon (lower) tracks are well resolved. In this display, hit size is proportional to energy loss within a 
strip. The increased energy loss of the proton as it slows and stops is clear. Note that for clarity the 
outer detector has not been drawn. 



with the CLAS collaboration to study hadronic final states in electron scattering on nuclear targets 
using existing JLab Hall B data. This analysis will allow theoretical models used in both electron and 
neutrino experiments to be tested. Other work in progress, with the Ghent|60| nuclear physics group, 
will develop the theoretical tools needed to extract the axial form-factor of the nucleon using MINERz^A 
quasi-elastic data on Carbon. The ultimate aim is to perform nearly identical analyses on both neutrino 
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Figure 21: Pauli suppression in a Fermi gas model for Carbon with binding energy e = 25 MeV and 
Fermi momentum kf = 220 MeV/c. A similar suppression is expected for quasi-elastic reactions in 
MINERi/A. 

and electron scattering data in the same range of Q^. 
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Figure 22: Extraction of Fermi gas model parameters, the effective Fermi momentum kj and nuclear 
binding energy e, from 500 MeV electron scattering data L57J . Distributions shown correspond to scat- 
tering from (a) Carbon, (b) Nickel, and (c) Lead. 
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Figure 23: The ratio of bound (in Carbon) to free nucleon form-factors for Fi, F2, and Fa from ref 
l59l . Binding effects on the form factors are expected to be small at higher (therefore, this model is 
not valid for > 1 (GeV/c)^). 
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7 Resonance-Mediated Processes 



Inclusive electron scattering cross-sections with hadronic mass W < 2 GeV exhibit peaks correspond- 
ing to the A(1232) and higher resonances at low (see Figure UM. This resonant structure is also 
present in neutrino scattering, although there is little data in this region. In addition to the natural 
interest in probing the nucleon weak current and axial structure via neutrino-induced resonance produc- 
tion, a better understanding of this process is essential for interpreting modern neutrino-oscillation and 
nucleon-decay experiments. This is particularly true for neutrinos in the region around 1 GeV, where 
single-pion production comprises about 30% of the total charged-current cross-section. 



E = 2.445 GeV ^ = 20.00° E = 3.245 GeV tS- = 26.98° 




Figure 24: Inclusive electron scattering showing the A and higher resonances. at the A peak is 
approximately 0.5, 1.5, 2.5 and 3.5 (GeV/c)^ for the four spectra, respectively 
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In this kinematic region, neutrino Monte Carlo programs have reUed on early theoretical predic- 
tions by Rein and Sehgal|70|. Recently Sato, Uno and Lee 1681 have extended a model of A-mediated 
pion electroproduction to neutrino reactions. Also, Paschos and collaborators, using the formalism of 
Schreiner and von Hippel L69J have included the effects of pion rescattering and absorption for resonance 
production in nuclei. 

7.1 Overview of Resonant Electroproduction 

In electron scattering, the behavior of the A(1232) transition form-factor is considered to be a primary 
indicator of the onset of perturbative QCD (pQCD). The behavior expected for a resonant spin-flip 
transition is dramatically different from the helicity conservation characteristic of perturbative descrip- 
tions. Comparison of the measured elastic and resonant form-factors reveals|71 72 T^EDE^ 84| 
that while the nucleon and higher-mass resonant form-factors appear to approach the predicted 
leading-order pQCD behavior around = 2 (GeV/c)^, the A (1232) transition form-factor decreases 
more rapidly with Q^. One possible explanation! 71 1 is that helicity-nonconserving processes are dom- 
inating. The A excitation is primarily a spin-flip transition at low momentum transfer, in which the 
helicity-nonconserving As amplitude is dominant l76ll . If the leading order Ai helicity-conserving 

2 2 

amplitude were also supressed at large momentum transfers, the quantity Q'^F would decrease as a 
function of Q^. 

Electromagnetic helicity matrix-elements correspond to transitions in which the initial state has 
helicity A and the final states have helicity A'. Transitions between a nucleon state | > and a resonant 
state \R > can be expressed in terms of dimensionless helicity matrix-elements ITTI : 



In this equation, the polarization vectors e+'^' correspond to right- and left-circularly polarized pho- 
tons, and longitudinally polarized photons, respectively. Following the formalism used by Stolerf72l 
and others, the differential cross-section may be written in terms of longitudinal and transverse form- 
factors Ge and Gt, as follows: 




(1) 




Ge and Gt are analogous to the Sachs form-factors for elastic scattering. In terms of the dimensionless 
helicity elements above. 



Ge = Go 



(3) 



and 



t*|Gt|' = ^(|G'+|2 + |G_|2) 



(4) 



where 



,2 



(5) 



The recoil factor free is given by 



rec 



EL 



(6) 
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R{W) is the familiar Breit-Wigner expression lTTll for the line-shape as a function of energy: 



"^(^^ - iw^-w^r+w^vi 

The mass and width of the resonance are Wr and Tr. 

Helicity is conserved in vector interactions of free, relativistic fermions. In the limit that a spin-| 
parton is massless and free, its helicity must be conserved in interactions with a vector gluon or pho- 
ton. At sufficient momentum transfer, the constituent quarks within a hadron can indeed be treated 
as massless and free, and the hadron helicity can be replaced by the sum of its constituent quark 
hehcities|78 , 79 1. Therefore, at high Q^, hadron helicity should also be conserved. 

For resonant electroproduction, the scattering can be analyzed in the Breit frame of the A = 3/2 
A resonance. The incoming virtual photon can have positive, zero, or negative helicity. The outgoing 
resonance helicity can be calculated from angular momentum conservation! 80 1 : 

Aa = A^ — Xjy (8) 

Hadron helicity is conserved when the incoming photon helicity is positive, and the A excitation 
emerges with the same helicity (1/2) as the initial nucleon state. This is described by the helicity 
amplitude Ai given by: 



A, = f-^G, (9) 

2 V K 

K is the energy of an equivalent on-mass-shell (real) photon producing a final mass state W: 

K = (1^2 - m2)/2M (10) 
Helicity is not conserved when Aa, given by 



2 



2 



A, = f-^G. (11) 



is the dominant amplitude. 

In terms of helicity amplitudes a dimensionless form-factor F may be defined where: 

= \Gt{Q')\' = -^^{Wl - M%)\Anm' (12) 

Here, 

\Ah{Q'')\^ = \Ar{Q^)\^ + \AdQ^)\^ (13) 

2 2 

At high Q^, the helicity conserving amplitude should dominate the helicity-nonconserving amplitude. 

A3 should be small compared to Ai according to pQCD. 
2 2 
In leading-order pQCD, two gluons are exchanged among the three pointlike quarks. These gluon 

exchanges ensure that the final quarks, like the initial ones, have low relative momenta, so that no 

powers of come from the wave functions. Form-factors calculated in the light-cone frame take the 

form |72l: 

F(g2) = / / dxdy^ixfTHHy) (14) 
JO Jo 
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where x and y are the initial and final longitudinal momentum fractions. $(x) and $(y) are the cor- 
responding quark distribution amplitudes and Th is the transition operator which is evaluated over all 
possible leading-order diagrams. This leads to the dimensional scaling rule llSlI 



G+(xAi(xQ-^, (15) 



2 



or 

F oc Q-^ (16) 

This dependence of the helicity amplitudes may be established up to factors involving lnf(D^) L82J . 
At high Q^, where the quark helicities are conserved, 

G+ oc Q-^ (17) 

777 

Go oc (18) 

and 

2 

777 

G_ oc i-Q^)G+ (19) 

The prediction that F{Q'^) oc l/Q^ if G+ is dominant can be understood by combining the above with 

the definitions of ^3 and ^1 in the dimensionless form-factor expression. 

2 2 

In addition to this dependence of the transition form-factors, pQCD makes definite predictions 
about the relative contributions of the magnetic dipole Mi+, electric quadrupole £'1+, and Coulomb 
quadrupole 5i+ amplitudes. In quark models at low Q^, the N — A transition is primarily due to a 
single quark spin-flip, requiring the Mi_|_ to be the dominant contribution! 83 1. At very low Q^, near 
zero, experiments have confirmed this prediction, evaluating £'1+ and Afi+ at the resonance position. 
However, as noted, only helicity-conserving amplitudes should contribute at high Q^, which leads to the 
prediction that the ratio £'i+/Mi+ = 1. Results from Jefferson Lab 1841 indicate that hadron helicity 
is not yet conserved at = 4 GeV^, finding the transition form-factor F to be decreasing faster 
than and continued Mi_|_ dominance. However, while pQCD apparently does not yet describe 
resonant excitation at these momentum transfers, it is not clear how constituent quark models can be 
appropriate at such high values, and regardless, no single model describes all of the data well. The 
Delta resonance, then, remains an object of intense study at facilities like Jefferson Lab and Mainz, with 
future experiments planned. 



7.2 Weak Resonance Excitation 

Sato and Lee|67| have developed a dynamical model for pion photo- and electroproduction near the 
A resonance which is used to extract N — A transition form-factors. Through this work, the afore- 
mentioned discrepenacy between the A transiton form-factor as calculated from a constituent quark 
model and the measured transition form-factor (a difference of about 35%) has been understood by 
including a dynamical pion cloud effect. Recently this work has been extended by Sato, Uno and Lee 
to weak pion production! 68 1- They show that the renormalized axial N — A form-factor contains large 
dynamical pion cloud effects which are crucial in obtaining agreement with the available data (in this 
case, on Hydrogen and Deuterium). Contrary to previous observations, they conclude that the N — A 
transitions predicted by the constituent quark model are consistent with existing neutrino-induced pion 



52 



production data in the A region. It is interesting to note that the pion cloud effect on the axial N — A 
form-factor is mainly to increase the magnitude. On the other hand, both the magnitude and the slope of 
the Mi+ are significanlty changed by including pion cloud effects. The authors cite the need for more 
extensive and precise data on neutrino-induced pion-production reactions to test their model and to pin 
down the Q^-dependence of the axial-vector N — A transition form-factor - data which MINERz^A can 
certainly provide. 

MINERz^A will measure scattering on nuclei, at least in the first years without a hydrogen target, 
and comparison to improved data on a free proton target will not be possible. Still, as discussed in 
Sectionnol the average dependence of the cross-sections (and, hence, structure functions and form- 
factors) will be magnified by the Fermi smearing of the resonant enhancements. It should be possible 
to map out the Q^-dependence of the axial-vector N — A form-factor. The work of Sato, Uno and Lee 
can be used as Monte Carlo input for MINERi/A, and should be essential to predictions of A excitation 
in nuclei which can be compared directly with MINERz^A data. 

7.3 Nuclear Effects 

Neutrino experiments rely heavily on detailed Monte Carlos to simulate the response of the rather 
complicated target / detector systems involved. The MINERi/A simulation will be greatly enhanced 
by accurate descriptions of the nuclear effects involved. The majority of hadrons produced in inelastic 
scattering are pions, and so the nuclear attenuation of these must be taken into account. In considering 
hadron attentuation results from HERMES, Gaskelll91 1 suggests that a good first step is the one time 
scale parameterization, which goes as (1 — z)ij. The A-dependence could then be taken into account via 
a simple A^^^ scaling in (1 — Ra), where Ra is the ratio of cross-section on nucleus A to deuterium. 

Another relevant nuclear effect, currently being applied in neutrino event generators for protons 
but not pions, is termed color transparency (CT). Color transparency, first conjectured by Mueller and 
Brodsky 1 85 1 refers to the suppression of final (and initial) state interactions of hadrons with the nuclear 
medium in exclusive processes at high momentum transfers. CT is an effect of QCD, related to the 
presence of non-abelian color degrees of freedom underlying strongly interacting matter. The basic 
idea is that, under the right conditions, three quarks (in the case of the proton), each of which would 
normally interact strongly with the nuclear medium, can form an object that passes undisturbed through 
the nuclear medium. This small object would be color neutral outside of a small radius in order not 
to radiate gluons. Unambiguous observation of CT would provide a new means to study the strong 
interaction in nuclei. 

Several measurements of the transparency of the nuclear medium to high energy protons have been 
carried out in the last decade. At Jefferson Lab, CT searches have concentrated on the quasi-elastic 
A{e, e'p) reaction which has several advantages in the search for CT. To date, A(e,e'p) experiments at 
SLAC 1.86J and JLab L87J have found no evidence for the onset of CT at momentum transfers up to 
8.1 (GeV/c)^. However, there is some potential evidence for CT in A{p, 2p) data from Brookhaven 
MM- 

It has been suggested that the onset of CT will be sooner in a qq system than in a three-quark system. 
Thus, the next best reaction in the expectation of CT is the A{e, e'lr) reaction. Current theoretical 
calculations suggest that most of this CT effect should be seen around = 10 (GeV/c)^, well within 
the MINERi/A kinematic range. This effect has not yet been considered in neutrino Monte Carlos, nor 
has it been well studied in other processes. However, it will be well-measured in the Jefferson Lab 
kinematic regime prior to MINERi/A L90J . and should then be incorporable into the Monte Carlo. 
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7.4 Exclusive Channels 



While there is a large body of inclusive (e, e') scattering data in the resonance region on hydrogen, deu- 
terium and nuclei, more exclusive measurements have been rare until recently. With JLab p{e, e'p)TT^ 
spectrometer measurements[72 84 1, the CLAS N* program|65| and CLAS ^^C(e, e'X) data, more 
exclusive reactions are becoming available. This data will help to "calibrate" the vector current part of 
weak resonance/meson production models and to extend Delta resonance models such as that of Sato, 
Uno and Lee to higher resonances. These exlusive measurements are also naturally of interest because 
even to make inclusive measurements with neutrinos, the full final state must be observed and recon- 
structed. With the expected statistics and resolution of MINERz^A, it should thus be possible to extract 
much more information about resonances than what is available in the inclusive channel. 

Figure |25] from the CLAS 1 64 1 is an illustration of the type of just part of the information available 
when one or more reaction fragments are detected in resonance region electron scattering. One item of 
interest in this data is a peak observed near W = 1.72 GeV in the spectrum for the pyr+vr^ final state. 
While an analysis of the angular distribution of this peak gives quantum numbers that agree with the 
PDG ^3/2+ (1720) state, the observed hadronic properties (coupling amplitudes) of this resonance are 
quite different from what is predicted from the PDG state. This illustrates that electro-weak excitiation 
of baryon resonances is an active field and that MINERi^A measurements are timely. 

7.5 Expected Results 

As shown in Tables |4] and |6l over 250,000 resonance production events, with useful statistics to at least 
= 4 (GeV/c)^, are estimated for 4 years of running. Approximately 40% of these events have W < 
1.4 GeV, so a good mix of events mediated by the A (1232) and higher resonances will be obtained. 
The resonance production measurements in MINERz^A can be grouped into several categories: 

1. Measurement of inclusive (f^, /i^) and (z/, v) spectra in the resonance region: As is done in the 
deep-inelastic region, this implies extraction of structure functions which can be compared to 
structure functions and form-factors from electron scattering. The experimental method is the 
same as for DIS events. For each event, we sum up the energies and momenta of the muon and all 
final state hadrons (either pions and nucleons, or all nucleons if the pions from resonance decay 
get absorbed on their way out of the nucleus) to get the total neutrino energy, and calculate 
and W, y etc. This kind of analysis will be done as a function of W for the entire resonance 
region. 

Although these measurements are done on a nucleus, we will be able to compare results to res- 
onance production predictions (such as Rein Seghal[70|) on nucleons, with some guidance from 
electron scattering. Because inclusive measurements are a sum over all final states, nuclear ef- 
fects should be primarily limited to Fermi motion and some Pauli Blocking. Despite the Fermi 
motion and resolution of MINERi/A, the Delta resonance will still be clear so it's form factor as 
a function of can be measured. The higher resonances will be smeared out, but still can be 
compared to "smoothed" behaviour of resonance models and to predictions from duality. One 
practical result of measurements above the Delta resonanance may be to modify the amount of 
non-resonant background in the resonance region models used in neutrino event generators. 

2. Examination of specific final state reaction products (single pion production, inclusive pion spec- 
tra): Specific final states, through the reactions: {v^, ^^p-K^), (f^, ^^nvr^), (f^, /i~p7r°), (f^, v^~mT^), 
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Figure 25: Invariant mass spectra from p{e, e'X) demonstrating the multi-hadron reconstruction capa- 
bility in the JLab CLAS spectrometer. I64II 



{^11,1^11 PT^^), and [u^^u^ n7r+), are useful in selection of a specific final state isospin. 

These final state measurements will rely on an improved understanding of final state interactions 
and will benefit from electron scattering hadron transparency studies and CLAS ^^C(e, e'X) data 
(which includes (e, e'pir^), and (e, e'n7r+), and which are equivalent to two of the above neutrino 
reactions). With these inputs, we will be able to map out the dependence of the Axial vector 
N — A form factor. But even without this better understanding, angular distributions should be 
less affected by final state interactions than overall cross sections. Thus, we will be able to extract 
ratios of weak transition amplitudes to compare to similar electron scattering amplitudes 

As measurement of detailed angular distributions of these final states is possible, the data on 
nuclear targets can also be used to study the feasibility of doing a phase shift analysis of the 
data if a hydrogen target is used in later phases of this experiment. This phase shift analysis, 
recommended by Sato, Uno and Lee[68|, like the JLab CLAS N* program, would be aimed at 
extracting the — A form factor model independently and providing a better understanding of 
low-lying nucleon resonances. 



55 





0.5 




0.45 


;cm 


0.4 




0.35 


b 


0.3 




0.25 


c 


0.2 


=1 n.i'i 


T 
c 


0.1 


> 






0.05 








• Radecky el a]. 
O Barish et al. 
4 Grabosh et aJ. 


















1 2 


3 


4 


5 


6 



EJGeV) 




Figure 26: Total pion production cross-sections. 



3. Using resonance production as a tool to study final state interactions: Having a selection of nu- 
clear targets helps here as the A dependence of the various reactions channels listed above can be 
studied. Another analysis that can be done along these lines is to measure inclusive pion spectra. 
Paschos et.al.i6j| combine resonance production and final state interactions to make predictions of 
pion spectra from neutrino scattering on nuclei. These spectra (Fig. |^ can be easily convoluted 
with neutrino beam energy distributions to produce pion energy distributions that can be directly 
compared with our data. 



7.5.1 Complementary studies at JLab 

The analysis of the above types of measurements will be closely coordinated with complementary ex- 
periments at Jefferson Laboratory (which are led by members of the MINERi/A collaboration). The 
following are the Jefferson Laboratory electron scattering experiments in Hall C that are connected 
with measurments of inclusive scattering in the resonance region at MINERz^A. 

1. JLab hydrogen experiment E94-110 (investigates inclusive F2 and R in the resonance region). 
C.E. Keppel spokesperson (data already taken). 

2. JLab deuterium experiment E02-109, investigates inclusive F2 and R in the resonance region. 
C.E. Keppel, M. E. Christy, spokespersons (approved to run in 2004). 

3. JLab experiment E99-118 investigates nuclear the dependence of F2 and R at low for high 
values of W. A. BruU, C.E. Keppel spokespersons (data already taken). 

4. Jlab experiment E02-103 hydrogen and deuterium resonance F2 data at high approved by Jlab 
PAC24 to run in 2004 (J. Arrington, spokesperon) 

5. Jlab experiment E04-001 to investigate F2 and R in the resonance region with nuclear targets. 
A. Bodek and C. E. Keppel, spokespersons (proposed to run in Hall C together with E02-109 in 
2004) to provide vector resonance form factors and R on the same nuclear targets that are used in 
neutrino experiments (e.g. Carbon, Iron, Lead). 
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Figure 27: Predicted Tr"*" energy distribution for scattering on ^^O of Paschos et. al.l6l. 

The following are collaborative programs between the electron scattering community that are con- 
nected with measurements of final states in the quasielastic region and in the region of the first resonance 
at MINERi/A. 

1. Steve Manly (Rochester) and Will Brooks (Jlab) program to use existing Hall B CLAS data at 
Jefferson Laboratory to study hadronic final states in electron scattering on nuclear targets (e.g. 
Carbon). 

2. Work with the Argonne group of Lee to model first resonance production in the region of the 
first resonance and also Ghent nuclear physics group in Belgium ['60'|, to model both electron 
and neutrino induced final states. In addition, there are other theoretical efforts (e.g. Sakuda and 
Paschos |6 J) on nuclear effects for the hadronic final states in the region of the first resonance. 

3. Comparison of electron scattering data (primarily proton and pion transparency measurements) to 
final state interaction models used in neutrino event generators such as NUANCE and NEUGEN|661. 

7.5.2 Resonant form-factors and structure functions 

The analysis of inclusive data in the resonance region with MINERz^A will be done using the standard 
structure function analysis techniques. The sum of neutrino and antinuetrino differential cross-sections 
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is used to do a Rosenbluth separation and extract F2 and R for a Carbon target. The difference between 
neutrino and anti-neutrino differential cross-sections is used to extract the structure function xF-^. The 
nuclear effects in the resonance region at low values of are not well understood. Electron scattering 
data show that duality works at greater than 1 (GeV/c)^ for hydrogen and deutrium targets. In 
addition, there are indications that the nuclear effects also scale with the Nachtman scaling variable. 
However, these observations have not been tested in neutrino scattering, nor have they been tested 
in neutrino or electron scattering at lower values of Q^. The information from Jefferson Laboratory 
proposal E03-110 will provide this information for nuclear targets for the vector structure functions. 
MINERi^A in turn will be able to extend these duality studies to the axial vector structure functions. 

At present, the axial form-factor for the first resonance is not very well known. MINERz^A will 
have a very high statistics sample in this region, which is equivalent to the sample for quasilelastic 
scattering described earlier. However, since MINERi/A data is on a Carbon target, nuclear effects must 
be understood. The theoretical tools used to model the nuclear effects in Carbon for the final state 
particles in the region of the first resonance in neutrino scattering, will be tested with CLAS Hall B 
electron scattering Jefferson Lab data on Carbon and other electron scattering data. 

7.5.3 Single-pion final states 

Using the angular distribution in the exclusive final states —>■ n^ir^p, we plan to fit for the reso- 
nant and non-resonant amplitudes. The extracted non-resonant amplitude should be consistent with the 
measured value of R in this region (extracted from the inclusive scattering sample). 

By using both neutrino and anti-neutrino data MINERi/A can investigate transitions into isospin 3/2 
states A++ andA^. An analysis of the ratios of various final states. p7r+, nvr^ and piiQ will provide 
additional information. As mentioned ealrier, we plan to do a comparison of resonance production 
with electron scattering on free nucleons to Hall B CLAS data with bound nucleons in Carbon. Within 
MINERi^A itself, we can compare the reactions up i/nvr"*" and vn up'K~ on bound nucleons 
directly, and investigate additional channels in order to better understand the effects of pion and nucleon 
rescattering. 

MINERz^A is expected to have good resolution for single pion events in the resonance region 
{W < 2 GeV). Figure |28] shows and W distributions of single pion events from CH2 in the 
MINERi^A Monte Carlo along with reconstructed distributions that take into account MINERi/A's 
energy resolution for hadrons (Figure l30b . While the Fermi motion in nuclei washes out higher res- 
onances, it is clear that Delta events can be readily identified and separated from higher resonances. 
This expected resolution implies that a differential cross-section da/dQ'^ for Delta production on Car- 
bon equivalent to that for Hydrogen (Figure l3lTl can be obtained with high statistics. Figure l29l shows 
an example of a charged-current neutrino interaction producing a A++ which decays to a pion and 
proton. Distinct muon, proton and pion tracks are all visible showing that the resonance can be well 
reconstructed. 
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Figure 28: W and reconstruction for events with a single 7r+. Top row are "true" W and distri- 
butions from the MINERi/A Monte Carlo. The second row are the reconstructed distributions assuming 
hadron energy resolutions from Figure |30| The invariant mass of the pion and highest energy proton 
give W which along with the muon energy and direction gives sufficient information to reconstruct Q^. 
Bottom row shows the correlation between the "true" and reconstructed quantities. 
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Figure 29: A++ production and decay in a charged-current neutrino interaction in the MINERz^A 
detector. Shown are (top track) the muon and (middle and bottom track) the pion and proton produced 
in the decay. Energy deposition is shown by hit size. For clarity the outer barrel is not shown. 
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Figure 30: Single charged pion resolution derived from MINERi/A Monte Carlo. 
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8 Coherent Neutrino-Nucleus Scattering 



The MINERz^A experiment has the potential to dramatically improve our knowledge of the dynamics of 
coherent neutrino-nucleus scattering. This process, in which the neutrino scatters coherently from the 
entire nucleus with small energy transfer, leaves a relatively clean experimental signature and has been 
studied in both charged-current (v^ + A —>■ +7r+) and neutral-current (u^+A — > u+tt") interactions 
of neutrinos and anti-neutrinos. Although the interaction rates are typically an order of magnitude or 
more lower than other single-pion production mechanisms, the distinct kinematic characteristics of these 
events allow them to be identified. Because the outgoing pion generally follows the incoming neutrino 
direction closely, this reaction is an important background to searches for —>■ v^. oscillation, as these 
events can easily mimic the oscillation signature of a single energetic electron shower. Neutral-current 
coherent production will be discussed in more detail in Section 113.51 here we limit our attention to 
the charged-current channel where the kinematics can be fully measured and the underlying dynamics 
explored. 
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Figure 32: Charged-current neutrino-carbon coherent cross-sections. Results have all been scaled to 
carbon assuming an 

Ai/3 dependence, and a{CC) = 2a{NC) CTOHt . 



8.1 Theory 

It is well known from electron scattering that at low and high v, vector mesons are abundantly 
produced through diffractive mechanisms. These interactions aie interpreted as fluctuation of the virtual 
photon intermediary into a virtual meson with the same quantum numbers, which by the uncertainty 
principle can travel a length 



63 



where m is the mass of the meson in question. For the weak current, similar fluctuations can occur, 
into both vector- and axial-vector mesons. From the Adler relation and "partially-conserved axial cur- 
rent" (PCAC) hypothesis, it is known that the hadronic current at low is proportional to the pion field. 
The hadronic properties of the weak current in these kinematic regions have been investigated through 
the study of nuclear shadowing at low x and the coherent production of vr, p, and ai mesons. Coherent 
scattering therefore allows investigation of the PCAC hypothesis and hadron dominance models of the 
weak current in detail f9T\. 

A number of calculations of coherent scattering, involving substantially different procedures and 
assumptions, have been made over the past thirty years l93 94 95 96 1. These calculations factorize the 
problem in terms of the hadron-like component of the weak current and the scattering of this hadron 
with the nucleus. The calculations assume PCAC as a starting point but quickly diverge when it comes 
to the number of hadronic states required to describe the weak current and how the hadron-nucleus 
scattering should be treated. The Rein-Sehgal model, used by both NUANCE and NEUGEN, describes 
the weak current only in terms of the pion field; the dependence of the cross-section is assumed 
to have a dipole form. Other calculations rely on meson-dominance models [95J which include the 
dominant contributions from the p and ai mesons. Figure|32lshows the coherent charged-current cross- 
section as a function of energy, compared to the model by Rein and Sehgal as implemented in NEUGEN 
and the calculation in L96J . 

8.2 Experimental Signatures 

The kinematics of coherent scattering are quite distinct compared to the more common deep-inelastic 
and resonant interactions. Because the coherence condition requires that the nucleus remain intact, 
low-energy transfers to the nuclear system, are needed. Events are generally defined as coherent by 
making cuts on the number of prongs emerging from the event vertex followed by an examination of 
the t distribution, where t is approximated by: 

-\t\ = -{q- p^f = {^i{Ei - p^j)f - {^^{pi)f (21) 

With its excellent tracking capabilities, the MINERz^A inner detector can measure this kinematic vari- 
able well. 

Figure |33] shows an event display of a coherent charged-current interaction in the MINERz^A inner 
tracking detector. Distinct muon and pion tracks are clearly visible and the vertex location is well 
defined. 

8.3 Expected Results 

To determine the ability of the MINERz^A experiment to measure the charged current coherent cross 
section, a Monte Carlo study was carried out using the GEANT detector simulation described else- 
where in this proposal. Analysis cuts were tuned on a sample of coherent interactions corresponding 
to that expected in a 3 ton fiducial volume for the integrated 4 year run (24630 events). Events were 
generated according to the appropriate mix of low, medium, and high energy running. This study used 
the Rein-Seghal |93 1 model of coherent production, as implemented in NEUGEN3. A 20k low-energy 
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Figure 33: A charged-current coherent event in the inner tracking detector of MINERi^A. For clarity 
the outer barrel detector is not shown. 

beam event sample was used for background determination. This sample included the appropriate mix 
of NC and CC events. Based on published bubble chamber analyses, it is expected that charged current 
reactions are the largest contributor to background processes, in particular quasi-elastic and delta pro- 
duction reactions where the baryon is not observed or is misidentified as a pion. To isolate a sample of 
coherent interactions, a series of cuts are placed on event topology and kinematics. 

Topological Cuts: an initial set of cuts are applied to isolate a sample of events which contain only 
a muon and charged pion. These cuts are based on the hit-level and truth information as provided by 
the GEANT simulation. 

1. 2 Charged Tracks: The event is required to have 2 visible charged tracks emerging from the 
event vertex. A track is assumed to be visible if it produces at least 8 hits which are due to this 
track alone. 

2. Track Identification: The two tracks must be identified as a muon and pion. The muon track is 
taken to be the most energetic track in the event which does not undergo hadronic interactions. 
The pion track is identified by the presence of a hadronic interaction. The pion track is required 
not to have ionization characteristic of a stopping proton (which is assumed can be identified 95% 
of the time). 

3. 7r°/neutron Energy: Because the MINERi/A detector is nearly hermetic we have also assumed 
that neutral particles will produce visible activity which can be associated with the event and 
cause it to be identified as not coherent. Events with more than 500 MeV of neutral energy (vr" 
or neutron) produced in the initial neutrino interaction are rejected. 
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Table 8: Existing measurements on coherent pion production l92ll . 
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7r°/neutron Energy 
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Table 9: Analysis cuts to isolate a sample of coherent interactions. The cuts are described in the text. 

4. Track Separation: In order to make good measurements of the two tracks, it is required that 
the interaction point of the pion be greater than 30 cm from the vertex and that at this interaction 
point at least 4 strips separate the two tracks in at least one view. 

Kinematic Cuts: because of the very different kinematics between coherent and background reac- 
tions, cuts on kinematic variables are very effective at isolating the final sample. In this analysis, the 
true pion and muon 4-momenta were used as the reconstruction values. For the final event rates we 
reduce our overall signal sample by 0.65 to roughly account for this assumption. 

1. x< 0.2: A cut is made requiring that Bjorken-x (as reconstructed from the observed pion and 
muon 4-momenta) be less than 0.2. This cut eliminates a large amount of the background coming 
from quasi-elastic reactions which have x^ 1 . 

2. t< 0.2GeV^: The most powerful variable for the identification of coherent events is the square 
of the 4-momentum transfer to the nucleus. The previous expression relating t to the observed 
particles in the event is used as the estimator of this quantity. 

3. Pt^i 600 MeV: Requiring > 600 MeV effectively eliminates backgrounds from delta produc- 
tion which tend to produce lower energy pions. 
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Applying this set of cuts to our signal sample we find that 7698 signal events pass all cuts, which 
gives an overall efficiency of 31%. Applying the factor 0.65 to account for the fact that we have not 
used fully reconstructed quantities for our kinematic cuts gives us a final event sample of 5004 events. 
Applying these cuts to the background sample we find that 12 events out of 20k pass all cuts. Nor- 
malized to the total event rate this gives an expected background of 4400 events. We note that in this 
analysis other important variables for background rejection, related to associated activity around the 
vertex, were not used. Figure |35] shows the expected precision of the MINERzvA measurement as a 
function of neutrino energy. Here we have only included the statistical error on the signal and assumed 
that the measured value is that predicted by Rein-Seghal. 

Another task for MINERi/A will be comparison of reaction rates for lead and carbon. The expected 
yield from lead will be 1800 charged-current events, assuming the same efficiency. The A dependence 
of the cross-section depends mainly on the model assumed for the hadron-nucleus interaction, and 
serves as a crucial test for that component of the predictions. No experiment to date has been able to 
perform this comparison. For reference, the predicted ratio of carbon to lead NC cross-sections at 10 
GeV in the Rein-Sehgal and Paschos models are 0.223 and 0.259, respectively 11071 . Figure l36lshows 
the predicted A-dependence according to the model of Rein and Sehgal. 
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Figure 34: Topological and kinematic quantities used to define the coherent sample. In all plots the 
solid histogram is the coherent sample and the dashed histogram are background processes. The relative 
normalizations of the two distributions in the initial plot is arbitrary, subsequent plots show the effect of 
the applied cuts. Top Left: Visible charged tracks. Top Right: Distance between the event vertex and 
the location of the pion interaction (in cm). Bottom Left: Bjorken-x as computed from the true pion 
and muon 4-momenta. Bottom Right: Square of the 4-momentum transfer to the nucleus (in GeV^) as 
calculated from the pion and muon 4-momenta. 
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CC Coherent Pion Production Cross Section 
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Figure 35: Coherent cross-sections as measured by MINERi^A compared with existing pubUshed re- 
sults. MINERi/A errors here are statistical only. 
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Figure 36: Coherent cross-sections as a function of atomic number. 
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9 Strangeness and Charm Production 



9.1 Overview 

The MINERz^A experiment in the NuMI near hall will allow high-statistics studies of the rich complex- 
ion of exclusive-channel strange-particle production reactions accessible in the I < < S GeV energy 
regime. We propose precision measurement of cross sections a{Ey) of exclusive associated-production 
reactions (A5 = 0) and Cabbibo-suppressed AS" = 1 reactions. The AS" weak hadronic current will 
be mapped out in detail, including its dependence, resonant structure, and polarizations of produced 
hyperons, to elucidate its coupling strengths and form-factors. A panoramic experimental delineation 
of all near-threshold f^-N strangeness production processes is envisaged which will motivate renewed 
efforts to formulate detailed models of these reactions. The resulting picture will have ramifications in 
other areas of particle physics, for example in estimation of atmospheric neutrino AS" backgrounds for 
nucleon-decay searches at megaton-year sensitivities. A MINERz^A exposure will also enable searches 
for new processes, e.g. unusual baryon resonances such as the recently reported candidate pentaquark 
state in K+n and K^p systems, and neutral-current strangeness-changing reactions. Extended running 
of the NuMI beam will allow exposures that will provide valuable complementary data for many 
neutrino topics. Anti-neutrino exposure will facilitate study of AS* = 1 single-hyperon production 
(A, S, Y*). Study of hyperon reactions will greatly extend the range over which the weak interaction 
form-factors which govern hyperon beta-decay can be examined. Thus a much better determination of 
the form-factors - especially of the three axial form-factors - will be possible. Hyperon polarization will 
provide additional analyzing power here, and the analysis will be free of the 'missing neutrino' problem 
which has hindered examination of the underlying V-A structure using semi-leptonic hyperon decays. 
As a natural extension of strange-particle production studies, we will search for strangeness production 
which accompanies dilepton processes. Such reactions have, in previous neutrino experiments, served 
as gateways to the study of charmed baryon production. 

The NOMAD experiment! 1 14 1 has studied inclusive strange-particle production extensively. MINERi/A 
will not improve significantly on those results, and the physics motivation for attempting to do so is un- 
clear. MINERz^Awill focus instead on exclusive channels; this is relatively unexplored territory, with 
the potential for high impact on future physics. 

9.2 Neutrino Strangeness Production Near Threshold 

In the threshold regime 1 < Ey < ^ GeV, neutrino interactions involving strangeness production yield 
final states containing either one or two strange particles. Exclusive i^-N channels comprise three 
categories, distinguished by reaction type (charged-current (CC) or neutral-current (NC)) and the net 
change in strangeness A5 = Sf — Si of the hadronic system (either AS = or AS = 1). 

The first category comprises charged-current AS = reactions initiated by f^. These are associated- 
production reactions where a strangeness -i-l meson (K+ or K*^) is produced together with a strangeness 
-1 hyperon (A or T,"^'^) or meson (K^ or K*^). Reactions of this category include: 





/i-K+AO 


(22) 




^ /i-vrOK+AO 


(23) 




^ /i-vr+ROAO 


(24) 


v^in - 




(25) 



71 







K+p 






■ROp 




At' 


vr+KOn. 



i^^p ^ /i-K+K%°p (26) 

Among charged-current AS* = reactions, reaction (I22t has the largest cross-section. This reaction, 
and reactions (l23l and (l24li as well, may proceed predominantly via N* production followed by strong 
decay into KA. 

Charged-current AS = 1 reactions make up a second category. For f^reactions, the resulting final 
states contain single strange K-mesons. The reaction cross-sections are Cabbibo-suppressed relative 
to AS = reactions involving similar hadronic masses. Additionally the AS* = AQ selection rule 
applies, and so the produced mesons are necessarily (K+, K*') and not (K^, K*^). Reactions of this 
category include 

(27) 
(28) 
(29) 

Reaction d27l has the largest cross-section among AS" = 1 exclusive f^-N reactions. 

Note that AS = AQ selection restricts AS' = 1 single-hyperon production to 17 rather than u 
reactions, e.g. 

F^N^^+ + (A,S,Y*). (30) 

Strange-particle AS = associated production can also proceed via neutral-current reactions. Ob- 
served channels include 

i/^p ^ i/K+A° (31) 
i/^n ^ i/K°A° (32) 
i/^n i/7r"K+A° (33) 

As with final states of d22l - d24b . it is similarly plausible that the hadronic systems of dSTt through d33l 
are dominated by intermediate N* states. 

9.3 Strangeness Production Measurements at Bubble Chambers 

Cross-sections for many associated-production and AS" = 1 reactions were obtained during the 1970's 
and '80s in experiments using large- volume bubble chambers exposed to accelerator neutrino beams. 
Principal experimental programs were the and exposures of the Gargamelle heavy-liquid (CFsBr) 
bubble chamber at CERN 1109 . 1 101 and the i'f_i,-D2 exposures of the 12-foot diameter bubble chamber 
at Argonne llllll and of the 7-foot diameter bubble chamber at Brookhaven 11121 . Typical samples 
involved less than ten observed events per channel, and cross-sections thereby inferred relate to one or 
a few bins in E,^. Contemporaneous theoretical/phenomenological treatments of reactions (l22l . d28l . 
(Ell, and (Ell can be found in 19 115. 116J . 

Cross-section ratios obtained by the bubble-chamber experiments provide rough characterizations 
of relative rates of occurrence among the strangeness reaction categories. For example, the frequency 
of strange versus non-strange hadronic final states in charged-current reactions is indicated by lllll : 

aiuN u"AK+ + u^p K) 
\ ^ — , = 0.07 ± 0.04 (34) 
a{uN fi-N + pion(s)) 
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The relative contribution of neutral-current versus charged-current reaction to threshold strangeness 
production is indicated by fl4l : 



(7(f^N i^pV^ + anything) 
(T(f/iN /i~V° + anything) 

and 



K+AO) 



0.22 ±0.14 (35) 



0.18 ±0.13 (36) 



cj(/i-K+A) + o-(/i-K+AXO) 

Perhaps the most significant "find" arising from bubble-chamber survey experiments was the first ob- 
servation of CC charmed-baryon production in a AS = 1 final state at BNL I113I : 

z/^p^//- S++ (37) 
S++ ^ A+7r+ ^ A°7r+7r+7r~7r+ (38) 

That excellent spatial resolution is a prerequisite for study of neutrino strangeness-production re- 
actions is illustrated by the bubble-chamber event in Figure |^ The figure shows the tracing of a 
photographic image recorded by one of four separate camera views of this f^-n interaction. The event 
shown was the first example of NC associated strangeness production via reaction ( H^ obtained using 
the deuterium-filled 12-ft diameter bubble chamber at ANL. Within the final state, the flight paths of the 

and and A" from the primary vertex to their respective "vee" decay points are 8.0 cm and 5.5 cm re- 
spectively mill . Fortunately it should be possible, with the lattice of triangular-cell scintillator tracking 
elements of a fine-grained detector, to achieve spatial resolutions near bubble-chamber quality (studies 
currently predict vertex resolutions of less than 1 centimeter in MINERi/A; see Section (15. 5.4t . This 
capability, together with dE/dx ionization imaging and momentum determination by ranging and ex- 
ternal magnetic tracking, will allow MINERzvA to explore exclusive strangeness-production processes. 



9.4 MINERz/A Samples Amenable to Hypothesis Fitting 

As described above, the available data on exclusive channel strange particle production by neutrinos 
is currently limited to samples of few tens of events isolated in bubble chamber experiments of the 
1970's and 80's. With the proposed MINERi/A program this data pool can be boosted by two orders- 
of-magnitude, thereby paving the way for comprehensive phenomenological treatments of neutrino 
strangeness production near threshold. In MINERi/A occurrences of i/^N exclusive strangeness pro- 
duction will comprise only a small fraction of the total event rate in the detector. However these events 
can be readily extracted from the accumulating total sample by exploiting MINERz^A capabilities: 

• The primary charged particle multiplicities are low for all strangeness production channels of 
interest; this feature is readily discernible event-by-event as result of MINERi^As fine granularity. 

• A prompt occurs in more than 50% of exclusive strangeness channels. The subsequent in- 
detector decay of the K+ some tens of nanoseconds later yields a signature in the light yield versus 
time profile of these events. The double-peak signature will be unmistakable in low-multiplicity 
events. 

• All channels of interest which do not have a K+ meson, have a final state K'' meson. In the 
latter reactions, those K'^'s which decay via vt+tt^ can be identified by examination of 
two-particle invariant mass and vertex displacement. 
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Figure 37: Trace of photograph from the ANL 12-ft diameter bubble chamber, of a neutrino neutral- 
current interaction in liquid deuterium yielding NC associated production z^KA. Flight paths to the vee 
decays of the two strange particles in the event are 8.0 cm and 5.5 cm in real space. 

TablefTOlsummarizes the sample populations for exclusive channel reactions obtainable by MINERi/A 
in the initial four-year run with the i/^beam. Note that the listed channel rates are restricted to sub-sets 
of events for which the imaged final states allow kinematic constraints (energy and momentum con- 
servation) to be imposed at the primary reaction vertices. That is, for final states which include V'^ 
particle(s), the only events tallied for Table ITOl are those for which each V'^ particle decays into two 
charged tracks, e.g. A'^ ^pvr^ and —>■ vr+vr^. Additionally, an overall detection efficiency is in- 
cluded which is based upon processing experience with strange particle production reactions in the ANL 
12-ft diameter bubble chamber u^D exposures II117I . Thus the sample populations estimated for Table 
[TOlrepresent events which can be reconstructed and treated using hypothesis fitting. Since the incident 
neutrino direction will be known relatively precisely, conservation of four-momentum will enable three- 
constraint fitting to charged current hypotheses and zero-constraint fitting to neutral current hypotheses 
(with Fermi motion of the target nucleon restricted to an allowed range in the fit). Given the various 
strange particle signatures in the reactions of Table flOl and given the four-momentum constraints which 
can be imposed at decay as well as at primary vertices, backgrounds from non-strange u^N interactions 
can be strongly mitigated. It should be possible to limit their contamination of total sample rates to well 
below 10% for most channels. 



9.5 Expected Results 

The paragraphs below summarize some specific topics involving neutrino strangeness-production pro- 
cesses that can be investigated using MINERz^A. 
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Table 10: Event populations for kinematically constrainable samples of exclusive-channel strangeness 
production reactions, obtainable in a four-year exposure of the three-ton inner fiducial volume of 
MINERi/A. 



Reaction Type 


Exclusive Channel 


No. Events (> constraint) 


AS = CC 


I^AtH - 




10,500 






^-^OK+AO 


9,300 




t'/xll - 


/x-^+K^AO 


6,300 




t'/xll - 


/x-K-K+p 


5,100 






;U~K°K+7r°p 


1,500 


AS = 1 CC 


i^MP - 


p 


15,900 




l^^tH - 


//-ROp 


2,400 






;u-^+kO 


2,100 


AS = NC 




zyK+A'J 


3,600 






z^kOAO 


1,100 






z^kOAO 


2,800 



9.5.1 Backgrounds to nucleon decay 

Current lifetime lower limits for nucleon decay (t//3 > 10^^ years) have not diminished hopes for the 
eventual success of supersymmetric grand unification (SUSY GUTs). Indeed, there is strong motiva- 
tion to proceed with more ambitious experimental searches. For the near future, improved searches will 
be carried out by Super-Kamiokande. Eventually these will be taken up by a next generation of un- 
derground detectors, e.g. by megaton-scale water Cherenkov experiments such as Hyper-Kamiokande 
and/or UNO llTTSl . 

Continued progress, either by improving limits to 10'^^ year lifetimes or discovery of nucleon decay, 
hinges upon improved knowledge of certain neutrino interactions which, when initiated by atmospheric 
neutrinos, can imitate nucleon-decay signals. The most problematic backgrounds to SUSY GUT modes 
arise via neutral-current associated production of strangeness at threshold energies. 

SUSY GUTs predict that nucleon-decay modes proceeding via virtual transitions involving inter- 
generational mixing are favored. Such modes yield final states containing strangeness +1 mesons, e.g. 

p ^ 1/K+ (39) 

n ^ z/K° (40) 

and possibly 

p ^ /i+K° (41) 

p ^ e+K° (42) 

Decays ( l39l and (l40l are thought to hold particular promise for first observation of baryon instabil- 
ity. 

Search for p — > jvK+ The leading nucleon-decay search experiment for the next decade (and per- 
haps longer) will be Super-Kamiokande. Its successor is also likely to be an underground water 



75 



Cherenkov detector with similar resolutions but a fiducial volume approaching megaton scale. In 
Super-Kamiokande, the search for proton decay mode ( B^ is currently carried out using three different 
methods, each motivated by the particulars of the final-state sequence being sought: 



i^0(7p + p + 8n) ^ + 7(6.3 MeV) 



TT^ + vr 







e ' z/zv 



-77 



e ' h'l' 



The three Super-Kamiokande approaches to finding proton decay are: 

i) K+ ^ /x+i/ spectrum search: Looks for an excess of single ^u-like ring events for which the 
reconstructed momentum matches that of two-body K+ decay at rest and the delayed rings 
accompanied by subsequent /i ^ e decay showers. This technique is already background hmited. 

ii) K+ — > /x+i/ gamma search: A candidate event has a signature 6.3 MeV gamma emitted by the 
parent nucleus together with a single /U-like ring having j?^ for a stopped K+ and accompanied by 
/u — > e decay. 

Hi) K+ — > vr+vr'^ search: Candidates have three rings compatible with vr+vr'' with vr*^ 77 from a 
stopped K^" and with a subsequent fi e decay signal. 

Neutrino background for p— > K~^u The combined search sensitivity for p— > uK~^ is dominated 
by the prompt gamma method ii) for which detection of a 6.3 MeV gamma from the nuclear de- 
excitation chain is crucial. Assuming this capability will be retained by next-generation underground 
water Cherenkov detectors, there is but one atmospheric neutrino reaction which may become an irre- 
ducible background in the search for this mode, and that is the neutral-current associated strangeness- 
production reaction dSTT i. That is, in an underground water Cherenkov detector, an atmospheric neutrino 
of or fe flavor may interact with a proton bound in an oxygen nucleus, producing a meson to- 
gether with a A hyperon and an (invisible) outgoing neutrino. Subsequently, the ^^N nucleus which 
is the remnant of the struck ^^O, de-excites producing the 6 MeV signature 7. The final state A is a 
target fragment and will most always have low momentum. When it decays into pvr^ as will happen in 
two-thirds of reaction d3n i occurrences, the daughter tracks will usually be below Cherenkov threshold 
and hence invisible. The final-state K"*" will subsequently decay, usually at rest, to yield a /i+ or vr+Tr*^ 
signature. Consequently the detection sequence in a water Cherenkov experiment indicated above for 
proton decay ( l39t is perfectly mimicked: 

u + ^^0 ^ V + K+ + + + 7(6.3 MeV) 

kP^t^ (below C threshold) 
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7r++7r0 

I 

It is crucial for future, and for ongoing proton decay searches as well, that neutrino background 
posed by (l3lT l and by other neutrino strangeness-production reactions be quantitatively understood. For- 
tunately, the relevant neutrino strangeness-production cross-sections, including their Ey dependence, 
can be precisely measured in MINERz^A. 




9.5.2 Measurement of a{uKK^) 

MINERi/Awill measure the exclusive AS" = neutral-current channel 

da 



dE^ 



v^^^u^K+k), (43) 



from its threshold at 1 GeV through its rise and leveling off to an energy-independent value at E^ 
between 10-15 GeV. For purposes of comparison and as a valuable check on svstematics llllll . the 
IS.S = companion charged-current reaction will also be measured: 

^{u^n^ li-K+K). (44) 

The off-line selections required to isolate reactions d43b and (l44T i are straightforward. Assuming the 
final-state A decays into p vr^ for these reactions, they share the following topological attributes: 



/) The reactions have relatively low charged-particle multiplicities from the primary vertex region. 
Reaction d43b has three charged prongs, including the two daughter tracks from A decay; reaction 
(l44l i has four charged prongs. 

//) The proton track of A decay will appear as a short, heavily-ionizing track from the vertex region 
which stops in the scintillator. 

The final-state K+ mesons will decay at rest or nearly at rest, and consequently a large-angle /i+ 
track wiU result. 



The most distinctive signature, however, arises with the time sequence for light emission in scintil- 
lator elements from these events. For reaction (l43b a "prompt" signal arises from the two-body decay of 
the A into charged tracks; in reaction (l44l the prompt burst is enhanced by the presence of the charged- 
current /x~ in the final state. The prompt signal is followed by a second signal, delayed by some few 
tens of nanoseconds, from two-body decay of the K+. This timing signature, taken in conjunction with 
the three topology attributes above, should yield clean samples of reactions (I43t and (l44li . 
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Figure 38: Time distribution from a neutrino interaction candidate for un fi^AK^, — > //+z^ 
recorded using the 1 kiloton water Cherenkov detector (IKT) at KEK. Two peaks, separated in time by 
a few tens of nanoseconds, signal the occurrence of a K'^ decay subsequent to the primary charged- 
current interaction. 

The feasibility of exploiting the signature afforded by the time profile of these reactions is illus- 
trated in Figure |38] The figure shows the time distribution from Cherenkov light from a candidate 
event for reaction (l44Ti . where the occurrence of two peaks separated by approximately 16 ns is readily 
seen [.119J . At K2K the effective energy reach of the KEK neutrino beam restricts cross-section mea- 
surements io Ey <3 GeV; the atmospheric neutrino flux however extends to higher energies. Thus the 
NuMI beam operated in the "low-energy" configuration will enable a complete picture of cj{Ey) for 
reactions d43l and (l44l to be obtained, providing an observational basis for future proton-decay searches 
to discover or set improved lifetime lower limits on decay modes favored by SUSY grand unification 
models. 

9.5.3 Strangeness-changing neutral currents 

Notably absent from the interaction categories of the previous paragraphs are neutral-current strangeness- 
changing reactions. These have never been observed; their occurrence at rates accessible in NuMI 
would imply new physics beyond the Standard Model. The existing limits on NC IS.S = 1 processes 
are based upon searches for rare K decays. Although there are experimental difficulties with unam- 
biguous identification of such processes in neutrino reactions, there is nevertheless an opportunity for 
strangeness-changing NC search in the neutrino sector. 

Hints that an unrecognized type of neutral-current processes may exist are to be found in discrep- 
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ancies involving hyperon weak radiative decays. These strangeness-changing weak decays have a clear 
disagreements between existing data and a variety of theoretical models - see II122LI123II for recent 
reviews. A long-standing puzzle concerns the large negative asymmetry coefficient observed in S"^ —>■ 
P7 decay, the measured value of which contradicts accepted notions concerning the size of SU(3)- 
breaking. To date, all of the assorted models invoked to describe these decays - including pole models, 
quark models, skyrmion models, vector meson dominance models and chiral models - fail to explain 
either the asymmetries observed or the decay rates of the various hyperons. Very recently, measure- 
ments of asymmetries which are large and negative in the — > and — > A°7 decays by the 
KTeV (Fermilab) and NA48 (CERN) experiments (|124| and 11251 respectively), run counter to the 
theoretical predictions for a sizable positive value [1261. According to the comprehensive analysis of 
Oilman and Wise (TtT], the hypothesis that all weak radiative hyperon decays in the 56 multiplet of 
SU(6) are driven by the single-quark short-distance transition s —>■ d'y, is untenable. 

A search for strangeness-changing neutral-cunent neutrino interactions can usefully clarify the ex- 
tent to which new physics parameters may be missing from the analysis of weak radiative hyperon 
decays. It is plausible that neutrino reactions, in contrast to hyperon weak decays, may provide cleaner 
signals for a new weak current in as much as the multiloop quark-gluon diagrams which complicate 
hyperon decay analysis would be absent. To hope for such a circumstance is perhaps not unreasonable; 
after all, the first clear evidence for existence of the Z° in the guise of neutrino NC reactions preceded 
the direct production of the by ten years. 

Below we list charged-current neutrino interactions which are examples of hyperon production; 
included are two-body final states which represent the inverse of hyperon beta-decay. These CC reac- 
tions require exposure of MINERz^A to an beam. Also listed are "companion" NC neutrino reactions 
which yield single final-state hyperons. The latter include possible strangeness-changing neutral-current 
reactions (labeled SCNC), a subset of which could be the focus of a dedicated search. Note that the 
SCNC reactions are in principle accessible with either or 17^ beams. 



z^"m + P^M+ + A° Uf, + p^iy + J:+ SCNC (45) 

^ ^+ + S° ^ 1/ + ^0 + S+ SCNC (46) 

^^+ + 7r° + S° -^u + K^ + p SCNC (47) 

+ n ^ /i+ + u + n^u + A° SCNC (48) 

^^+ + 7r-+A° + SCNC (49) 



To isolate SCNC interactions, it is of course necessary to distinguish them amongst the predominate 
neutrino CC and NC samples. For certain selected SCNC reactions, this should be feasible. First and 
foremost on our list to identify SCNC process is the absence of a charged lepton since they are only in 
the neutral-current reactions in conjunction with only one strange particle being present. Other methods 
at our disposal is the existence of hyperon resonances without an accompanying meson which requires a 
highly-segmented detector with excellent containment of neutral mesons. Background estimates under 
the hyperon resonances can be accurately determined by off-resonance measurement of pvr^ states that 
would then give the accuracy needed for the resonance region. Events above that expectation would 
yield limits on the SCNC processes. Prerequisite detector requirements are good resonance mass re- 
construction, neutral meson containment and a magnetic field, knowledge of the sideband backgrounds 
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not going through resonances to an accuracy 10 x better than the resonance search region. All of this is 
achievable in the current design of MINERi/A. 



9.5.4 Hyperon beta-decay and inverse neutrino processes 

Hyperon beta-decay A —>■ B e~ Ve provides a window onto weak hadronic current form-factors and 
their underlying structure. In the V-A formulation the transition amplitude is: 

M = -^<B\J^\A> Ue7A(l + l^)uu (50) 

The V-A hadronic current can be written as: 

<B\J^\A>=Ct u{B) { /i(^V + /2('z')^ + /3(g')^ + 

[ gi{q')i^ + g2{q^f-^+gz{q')^]i, }u{A) 

where C is the CKM matrix element and q is the momentum transfer. There are 3 vector form-factors: 
/i (vector), /2 (weak magnetism) and /s (an induced scalar); plus 3 axial-vector form-factors: gi (axial- 
vector), g2 (weak electricity) and (an induced pseudo-scalar). 

Recent high statistics measurements of these form-factors using KTeV hyperon beta-decays have 
been reported! 128 I; the results show that the level of SU(3) breaking is very small compared to expec- 
tations of modern theories! 129 1. These new results have been used to extract the CKM matrix elements 
Vus|130|. Similar physics studies can be done with anti-neutrino interactions that produce hyperons. 
The hyperon decays themselves will have the added feature of a self-analyzing power of the polariza- 
tion vector. Thus the fundamental form-factors and CKM matrix elements will be accessible without 
the hindrance of double solutions due to the missing neutrino energy. On the other hand, in MINERz^A 
there arises the problem of dealing with the nuclear potentials which comprise the environment for tar- 
get nucleons. This consideration might motivate running with liquid Hydrogen and Deuterium targets 
in a future program. 

Although the simplest beta-decays and their corresponding inverse processes provide the predomi- 
nate samples for both hyperon beta-decays and in AS = 1 neutrino interactions, there aie also interest- 
ing albeit more complicated 4-body beta-decay processes Usted below along with some corresponding 
strangeness-producing neutrino interactions: 



A p7r°e"P + P ^ ^+A7r° (51) 

S+ pvr^e+i^ (52) 

— > p7r^/_i^j/ i^/i + P — > /i~S'''vr"^ (53) 

^i^ + V^ fi'pK+ (54) 

E~ pTT^e^P i^/i + P ^ /i+S^7r+ (55) 



Although none of the 4-body hyperon beta-decays have been officially observed, a handful of il^ 
candidate decays may have been isolated; preliminary results were presented at DPF 2003| 131 1. The 
theory behind these decays, LI 32 J and 11331, with their more complicated interaction models, were 
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developed in the 60s. With neutrino interactions these processes should be easier to obtain and hence 
studied for the strength of physics interactions. They would allow for a much more complicated form- 
factor analysis process that involves 16 variables in the V-A style, but in the Standard Model would be 
a way to check the scale of the SU(3) breaking. It should also be noted that there are several types of 
beta-decays in this list that can be studied in neutrino beams before going to anti-neutrinos; these are 
the types shown in Equations |53] and |54] and others not listed. 

There are also forms of 4-body hyperon beta-decays and likewise neutrino interactions that are 
forbidden by A5 = — AQ. These too can be extensively searched for in neutrino-beam running before 
the necessity for anti-neutrino beams. While the existence of the forbidden decays would be exciting in 
hyperon beams, these interactions in neutrino beams would give information about probing the forms 
of the interactions. 

Study of AS = 1 pentaquark states, like those recently announcedl 134|, could be greatly extended 
here. In regard to these pentaquark states (whether 4-quark and a anti-quark bound combination or a 
loosely-bound baryon-meson combination similar to mesonic atoms), with the production of hyperons 
and mesons together a wealth of combinations can be throughly examined for studying the full spectrum 
of the pentaquark familv L135,l of particles as well as other exotic quark combinations such as di-baryons. 

9.5.5 Charm production physics 

Historically, most neutrino scattering experiments found their way into charm-production studies when 
they investigated opposite or same-sign muon pairs generated by neutrino interactions. This signal 
arises because many of the charm particles decay with a muon, giving an extra muon along with the 
one produced by the CC neutrino reaction. The decay muons usually differ substantially from the direct 
CC neutrino muons in both momentum and angular distribution, but in some cases it is not possible 
to discern a difference. In MINERi/A, with its lower neutrino energy beam, the production of charm 
particles will be suppressed compared to previous high-energy physics experiments. Hence the reach of 
MINERi/A will be limited, but its large neutrino flux still allows interesting charm physics to be done. 

An important contribution MINERi/A will provide in charm production is study of the cross-section 
turn-on at or just a few hundred MeV above threshold. This threshold is very sensitive to the bare charm 
mass. With the proposed beam running schedule for MINERi/A we expect ~ 6500 charm events for a 
three-ton detector over the first five years, with an additional ~ 3200 from anti-neutrino beam running 
for Xi? > 0. Most of these charm events (~ 65%) will be produced during the HE beam running 
configuration. As noted, these yields depend strongly on the bare charm mass; varying this parameter 
by 10% results in expected yield changes of 30%. As discussed earlier, neutrino experiments measure 
charm-production parameters by studying opposite-sign dimuon production. From preliminary studies, 
the expected number of dimuons in MINERi/A over five years is 530 ± 50 for a bare charm mass of 
1.3 GeV/c^. For bare charm masses of 1.15 GeV/c^ or 1.45 GeV/c^ the expected yields are 680 it 60 
and 420 ± 40 respectively. The yield assumes charm produced with xp > and a lower momentum 
cut on the decay muon of 1.5 GeV/c. The errors on the yields include the error on the average semi- 
leptonic branching ratio for charm L137 1 and the error on subtracting the background rate from pion 
decay. MINERi/A is at an advantage in being able to determine the sign of the muons via magnetic 
tracking. Background rates can be determined by looking for same-sign dimuons. At MINERi/A beam 
energies the expected number of background events should be approximately equal to the signal values. 
MINERi/A should improve on the charm-quark mass determination currently set by the NuTeV/CCFR 
data at 1.38 ± 0.13 GeV/c^fTTTll. 
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10 Perturbativ^Non-Perturbative Interface 



10.1 Parton Distribution Functions 

One obvious reason for the importance of neutrino data in the extraction of parton distribution functions 
is the neutrino's ability to directly resolve the flavor of the nucleon's constituents: v interacts with d, s, u 
and c while the V interacts with u, c, d and s. This unique abiUty of the neutrino to "taste" only particular 
flavors of quarks assists the study of parton distribution functions. A high-statistics measurement of the 
partonic structure of the nucleon is here proposed, using the neutrino's weak probe, to complement 
on-going study of this subject with electromagnetic probes at other laboratories. 

With the high statistics anticipated in MINERz/A, as well as the special attention to minimiz- 
ing neutrino beam systematics, it should be possible to determine the individual structure functions 
Ff ^(x, Q2)^ Ff^(.T, Q2)^ F2^^(x, Q% if ^(x, Q% xF^^{x, Q2) and xF^^ix, Q^) (where N is an 
isoscalar target) for the first time. 

In leading-order QCD, four of the structure functions are related to the parton distribution functions 
(PDFs) by: 

2F^'^{x,Q'^) = u{x) + d{x) + s{x) + u{x) + 
d{x) + c{x) 

2Ff^(x,Q2) = u{x) + d{x) + c{x) + u{x) + 
d{x) + s{x) 

xF^^{x,Q^) = u{x) + d{x) + s{x) - u{x) - 
d(x) — c(x) 

xFi^{x, Q^) = u{x) + d{x) + c{x) - u{x) - 
d{x) — s{x) 

Taking differences and sums of these structure functions allows extraction of individual parton dis- 
tribution functions in each (x, Q^) bin: 



2Fr-2Fr = [s{x)-s{x)\ + [c{x)-c{x)\ 

2F^^-xF^^ = 2[-u(x)+J(x) + c(x)] 

2Fl^-xFi^ = 2[u{x) + d{x) + s{x)] 

xif^-xif^ = [s(x) + s(x)] - [c(x) + c(x)] 



As the order of QCD increases and gluons are taken into consideration, global fitting techniques 
must be applied to extract of the parton distribution functions. With the manageable systematic errors 
expected for the NuMI beam, the ability to isolate individual parton distribution functions will be dra- 
matically improved by measuring the full set of separate v and i> structure functions with the impressive 
statistics possible in this experiment. 

There are two primary (related) methods for extracting this full set of structure functions. One 
exploits the varying y behavior of the coefficients of the structure functions in the expression for the 
cross section: 



83 



dxdy 



IT 



xy^F^^'\x,Q^) + 



2E, J 

the other uses the "helicity representation" of the cross section: 



^^F^^'\x,Q')± 



y(l-y/2)xF3^(")(x,Q2) 



and 



dxdQ"^ 



d^a^u] 
dxdQ"^ 



]^{f;^{x,q^) + xf^{x,q^)) + 



2ttx 
2 



{F^{x,Q^)-xF;^{x,Q^)) 



2'KX 



{F^{x,Q^)-xFUx,Q^)) + 



{Fi{x,Q^)+xFl{x,Q^)) 
2y^F£(x,Q2)", 



2 

..2 c^P/ 



By plotting the data as a function of (1 — y)^ in a given x — bin, it is possible to extract all six 
structure functions. 

For this sort of parton distribution function study, large anti-neutrino samples are an imperative. 



10.2 Quark Distributions at Large x 

Although a large body of structure function data exists over a wide range of x and Q^, the region 
X > 0.6 is not well explored. For x > 0.4 contributions from the qq sea become negligible, and the 
structure functions are dominated by valence quarks. 

Knowledge of the valence quark distributions of the nucleon at large x is vital for several reasons. 
The simplest SU(6) symmetric quark model predicts that the ratio of d to u quark distributions in 
the proton is 1/2, however, the breaking of this symmetry in nature results in a much smaller ratio. 
Various mechanisms have been invoked to explain why the d{x) distribution is softer than u{x). If 
the interaction between spectator quarks is dominated by single-gluon exchange, for instance, the d 
quark distribution will be suppressed, and the d/u ratio will tend to zero in the limit x —>■ 1 11481 . 
This assumption has been built into most global analyses of parton distribution functions [1491, and 
has never been tested independently. On the other hand, if the dominant reaction mechanism involves 
deep-inelastic scattering from a quark with the same spin orientation as the nucleon, as predicted by 
perturbative QCD counting rules, then d/u tends to w 1/5 as x — > 1 11501 . 

Measurement of structure functions at large x will yield insights into the mechanisms responsi- 
ble for spin-flavor symmetry breaking. In addition, quark distributions at large x are a crucial input 
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for estimating backgrounds in searches for new physics beyond the Standard Model at high energy 
colliders|151 1. 

The QCD evolution of parton distribution functions takes high-a^^j pdf's at low and evolves 
them down to moderate-and-low x at higher . This obviously means that one of the larger contribu- 
tions to background uncertainties at LHC will be the very poorly-known high-x PDF's at the lower 
values accessible to the NuMI beam. The appearance of an anomaly at high x will be discussed below. 
First note that one problem in studying this problem has been accumulation of sufficient data at high 
X, off light targets, to extract the PDF's. The NuMI beam will finally yield the necessary statistics to 
address this important concern. 

Uncertainties in current nucleon parton distribution functions at high x are of two types: the ratio 
of the light quark PDF's, d{x)/u{x), as x — > 1, and the role of leading power corrections (higher twist) 
in extraction of the high x behavior of the quarks. 

Analyses of present leptoproduction data from hydrogen and deuterium targets have been unable to 
pin down the high-x behavior of d{x)/u{x). Part of the problem is due to the still unknown nuclear 
corrections involved in extracting the "neutron" results from deuterium |152|. An analysis by Bodek 
and Yangri531 indicated that the d{x)/u{x) quark ratio approaches 0.2 as x ^ 1. However global 
QCD analyses of experimental results, such as the CTEQ fits|154|, do not indicate the need for this 
higher value of d{x = l)/n(x = 1). Besides the statistical and experimental uncertainties in the 
existing data, a complication with past experimental results was to model nuclear binding effects in 
the deuterium target used. These issues could be avoided with high-statistics exposure of an H2 target, 
which could directly measure the d{x)/u{x) ratio in protons as x — > 1 from the ratio of neutrino-proton 
to antineutrino-proton cross-sections. Such a measurement would require only a small correction for 
the residual sea-quark contributions at high x. 

Measurement of quark densities at high-XBj is closely related to the question of the leading-power 
corrections known as "higher twist effects". The n*^ order higher-twist effects are proportional to 
l/Q^" and reflect the transverse momentum of the quarks within the nucleon and the larger size of 
the probe as decreases, increasing the probability of multi-quark participation in an interaction. As 
for the d/u ratio, different analyses of higher-twist corrections in current data leave unresolved issues 
that new experimental information would clarify. Recent work by Yang and Bodekl l55] seems to indi- 
cate that what has been measured as "higher-twist" in charged-lepton scattering analyses is essentially 
accounted for by increasing the order (NNLO) of the perturbative QCD expansion used in the analysis. 

The only actual measurements of a higher-twist term in neutrino experiments have been two low- 
statistics bubble-chamber experiments: in Gargamelle| 156 1 with freon and in BEBC| 157 1 with NeH2. 
Both bubble-chamber analyses are complicated by nuclear corrections at high-x. However, both found 
a twist-4 contribution smaller in magnitude than the charged leptoproduction analysis and, most signif- 
icantly, negative. 

There are several indications that current parameterizations of the PDFs are not correct at high x. 
Figure |39l shows the ratio of measured Drell-Yan pair production lll58l compared to the latest CTEQ 
global fits, CTE06 I159I . The comparison seems to indicate that the valence distributions are overesti- 
mated at high-XBj ■ This directly contradicts a recent analysis at Jefferson Lab which seems to indicate 
that the valence distributions are underestimated at high x, as shown in Fig urel40l 

Efforts are underway to understand how the d{x)/u{x) ratio enters into the experimental compar- 
ison just discussed, and the large sample of high x events in MINERz^A would certainly help clarify 
these results. 

The principal reason that the d{x)/u{x) ratio is not better known is the difficulty of accessing the 
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structure of the neutron, due to the absence of free neutron targets, and the substantial theoretical uncer- 
tainties associated with extracting information from neutrons bound in nuclei. To overcome this prob- 
lem, the BONUS experiment at Jefferson Lab [1601 has been approved to measure the inclusive electron 
scattering cross section on an almost-free neutron using the CEBAF Large Acceptance Spectrometer 
(CLAS) and a novel recoil detector with low momentum threshold for protons and high rate capability. 
This detector will allow tagging of slow backward-moving spectator protons with momentum as low as 
70 MeV/c in coincidence with the scattered electron in the reaction D{e, e'ps)X. This will ensure that 
the electron scattering took place on an almost free neutron, with its initial four-momentum inferred 
from the observed spectator proton spectrum. These measurements will unambiguously provide neu- 
tron structure measurements, which will thereby also reveal which of the available models best describe 
for instance, on-shell extrapolation for neutrons in nuclei. 

It should be stressed that the BONUS experiment at Jefferson Lab will provide complementary in- 
formation to MINERi^A measurements, overlapping in kinematics, and on a similar time scale. With 
BONUS and MINERz^A combined, most of the questions in large-x nucleon structure, parton distribu- 
tions, and medium modifications, will be solved in the coming decade. BONUS will provide vital input 
regarding the extraction of neutron information from nuclei, while MINERz^A can uniquely provide 
flavor decomposition information. 

10.3 Quark/Hadron Duality 

The description of hadrons in terms of their fundamental quark and gluon constituents is one of the 
major challenges in nuclear physics today. While at present the quark and gluon degrees of freedom in 
QCD cannot describe the structure and interactions of hadrons directly, in principle it should be just a 
matter of convenience whether to describe a process in terms of quark-gluon or hadronic degrees of free- 
dom. This idea is referred to as quark/hadron duality, and means that one can use either set of complete 
basis states to describe physical phenomena. At high energies, where the interactions between quarks 
and gluons become weak and quarks can be considered asymptotically free, an efficient description 
of phenomena is afforded in terms of quarks; at low energies, where the effects of confinement make 
strongly-coupled QCD highly non-perturbative and the final state is guaranteed to consist of hadrons, it 
is more efficient to work in terms of collective degrees of freedom, the physical mesons and baryons. 
The duaUty between quark and hadron descriptions reflects the relationship between confinement and 
asymptotic freedom, and is intimately related to the nature of the transition from non-perturbative to 
perturbative QCD. It has been said that (short of the full solution of QCD) understanding and control- 
ling the accuracy of the quark-hadron duality is one of the most important and challenging problems for 
QCD practitioners todav ll38l . 

Although the duality between quark and hadron descriptions is formally exact in principle, how du- 
ality is manifest, specifically, in different physical processes and under different kinematical conditions 
is a key to understanding the consequences of QCD for hadronic structure. The phenomenon of duality 
is quite general in nature and can be studied in a variety of processes, such as e+e" — > hadrons, or 
semi-leptonic decays of heavy quarks. Duality in lepton-nucleon scattering, historically called Bloom- 
Gilman duality, links the physics of resonance production to the physics of deep-inelastic scaling. This 
duality is illustrated in Figure fTo!3l where the nucleon transverse {2xFi) and longitudinal (F^) struc- 
ture functions, measured in electron-proton scattering, are plotted as a function of the Bjorken scaling 
variable x for the indicated bins. The curves are a fit to the resonance data by Liang, and the parton 
distribution function based parameterization of the MRST Iil39J group at next-to-next-to leading order. 
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corrected for target mass ll40l . The data are in the resonance (from Hall C at Jefferson Lab| 141 1) and 
deep-inelastic (from SLAC I142II ) regimes, as indicated. Duality appears here in the observation that the 
hadronic (resonance) and quark (scaling) strengths are, on average, equivalent. Moreover, this is true 
for all bins observed, and thus the perturbative curve (MRST) apparently describes the average 
dependence of the hadronic, non-perturbative, resonance enhancement region. 
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Figure 41: The nucleon transverse {2xF\) and longitudinal (F^) structure functions, as measured in 
electron-proton scattering, are plotted as a function of the Bjorken scaling variable x for the indicated 
bins. The curves are a fit to the resonance data by Liang (light blue), and the parameterization from 
MRST| 139 1 (dark blue) at next-to-next-to leading order, corrected for target mass| 140 1. The data are in 
the resonance (from Hall C at Jefferson Lab L141J . purple) and deep-inelastic (from SLAC L142J . black) 
regimes, as indicated. 

The proposed MINERzvA experiment is uniquely poised to provide a wealth of data to answer 
where duality works, in what structure functions, in what reactions, and for what kinematics. Duality 
has been well- verified for the proton F2 structure function ll43L observed recently in the separated 
longitudinal and transverse unpolarized structure functions! 141], on nucleons and in nuclei| 162], and 
in polarized structure functions! 144 !. While its fundamental cause remains a mystery, duality appears 
experimentally to be a non-trivial property of nucleon structure. It is, therefore, crucial to test it in 
a variety of reactions - including neutrino-nucleon and -nucleus scattering and the structure function 
XF3. Duality studies of electron-deuteron scattering at low found a resemblance to deep-inelastic 
neutrino-nucleus scattering at much higher Q^, indicating potential sensitivity of duality to the valence 
quarks 1J45J. MINERzvA will allow this observation to be verified and tested for the first time, as data 
from similar kinematic regimes but differing in probe and interaction (from MINERi/A and Jefferson 
Lab) may be compared directly. 

It is important to point out that a revolutionary application of duality, if one understands the work- 
ings of the resonance-deep-inelastic interplay, would be to open-up the region of very high x, which 
has not been possible in any other experiment. As discussed above, the region of « 1 is an important 
testing ground for understanding of the valence quark structure of the nucleon, and it will allow us to 
discriminate between various models for the mechanisms of spin-flavor symmetry breaking in the va- 
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lence quark distributions of the nucleon. A first attempt at such an application is the recent analysis by 
Bodek and YangI 8 1, offering a promising procedure for fitting F2 in the low Q^, high x region. Extrap- 
olating their results through the resonance region yields values of F2 consistent with duality arguments 
and the Jefferson Lab results mentioned above. In addition, with nuclear targets, duality extensions 
to large x would permit measurements of the nuclear-medium modification of the nucleon structure 
function (nuclear EMC effect) at large x, where deviation of the ratio of nuclear to nucleon structure 
functions from unity is largest, and sensitivity to different nuclear structure models greatest. 

Members of the MINERi/A collaboration ai^e currently investigating quark/hadron duality in high- 
statistics electron scattering at Jefferson Lab with the same or similar nuclear targets as those proposed 
for MINERt/A ll61l . This will be followed by a comparison with all existing neutrino data, with the 
aim of continuing these studies with the higher statistics MINER;/A neutrino experiment in the future. 
Note that investigation of quark/hadron duality in the axial structure functions of nucleons and nuclei 
with neutrinos also adds a new dimension to the previous electron studies. Many issues, such as nuclear 
dependencies, should be well understood in advance of the MINERi/A data. 

10.4 QCD Moments 

Figure |42ldepicts the substantial enhancement in the kinematic domain of precision data made possible 
by MINERz^A over a range in x and Q^. This data will serve a variety of purposes, and address long- 
standing questions regarding structure function behavior at low Q^. Perhaps most importantly, the 
range of the data will allow for accurate moments of the structure functions to be obtained. To obtain 
a structure function moment, it is necessary to integrate over the full range in x at a fixed value of Q^. 
The Cornwall-Norton moment of a structure function F, for instance, is expressed as: 



The moments are fundamental quantities, calculable in QCD and recently calculated in lattice QCD 
at = 4 (GeV/c)^ for valence distributions ! 146 1 . If duality is shown to hold, the proposed data may 
provide one of the few available quantities that can be directly compared to lattice QCD calculations. 

Bloom-Oilman duality can be formulated in the language of an operator product expansion (OPE) 
of QCD moments of structure functions, in which contributions are organized according to powers of 
The leading terms are associated with free quark scattering, and are responsible for the scaling of 
the structure function. The 1/Q^ "higher twist" terms involve interactions between quarks and gluons 
and hence reflect elements of confinement dynamics. Duality measurements have been explained in 
terms of a weak dependence of the low moments of the structure functions f l4Tl. This is interpreted 
within the OPE as indicating that non-leading, 1/Q^-suppressed, higher-twist interaction terms do not 
play a major role even at low (« 1 GeV^). It is this interpretation that facilitates comparison to 
lattice calculations, as the latter have no higher twist effects included. 

Large-x (resonance region) data become increasingly important for higher-order moments due to 
the n — 2 weighting of the moment. At n=6, for example, the resonance and large x region above 
X = 0.7 make up 70% of the Cornwall-Norton moment of F2 at = 10 (GeV/c)^. The contribution 
is larger at = 4 (GeV/c)^, where lattice calculations are available. As noted above, there currently 
exist little to no neutrino resonance cross-section data in the resonance region or at larger x, while such 
data will be easily obtainable with MINERz^A. 




(56) 
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Figure 42: Left: Distribution of DIS events with LE beam in MINERi/A Monte Carlo. Right: Events 
where total hadronic energy is contained by MINERi/A. The line is an estimate of the Umit where 50% 
of events do not have containment of hadronic energy. 



It is important to reiterate that, regardless of duality or OPE arguments, the experimental values 
for the moments can be unambiguiously obtained with MINERi/A. For example, it is straightforward 
to note, from Figure |^ that even the low-energy beam provides data covering a large range in x (or 
W"^) for each value up to 10 (GeV/c)^ The higher -energy beams will complement this sensitivity, 
extending the range over which moments can be obtained, and adding statistics to the much of the 
region covered by the low-energy beam. While comparable coverage can be obtained by combining 
electron and muon scattering data from a multitude of laboratories, MINERi/A will uniquely provide, 
for instance, the xFs structure function, valence sensitivity (necessary to current lattice comparisions), 
and flavor decomposition. 

10.5 Expected Results 

The proposed studies of structure function moments and quark/hadron duality are straightforward with 
the proposed MINERi/A experiment. These topics do not have the demanding experimental constraints 
that many of the other proposed topics do. While it is crucial to understand the projected W"^ or x 
resolution, for instance, in studying resonance production behavior, duality studies and moment extrac- 
tions average over these kinematic variables and are therefore virtually insensitive to resolution issues. 
The expected MINERzvA resolutions are more than adequate both to form integrals such as that in 
Equation|56land to study duality using data comparison with perturbative predictions, as in Figure fTOJl 
Additionally, it is not necessary to isolate specific production processes in these studies. It is only 
total cross-section averages that are of interest, making MINERz/A essential to this effort in neutrino 
scattering. 

It has been observed from studies of quark/hadron duality using nuclear targets that data in the 
resonance regime scale even more obviously when smeared by the nuclear Fermi momentum ll62L as 
shown in Figure |^ In Hydrogen, the resonance peaks are prominent, while they are much less so in 
Deuterium, and completely smeared away in Iron. In all cases, however, the resonance region averages 
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to the scaling curves. In Iron, the smearing is such that the resonance data and scaUng curves overlap 
completely; the nucleus performs the duality averaging. 
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Figure 43: Structure function F2 data in the resonance region on Hydrogen (top), Deuterium (center), 
and Iron (bottom) covering a range in 0.8 < Q2 < 3 3 (GeV/c)2, and plotted as a function of the 
Nachtmann scaling variable ^. The elastic (quasi-elastic) peaks have been removed. The curves are the 
MRST and NMC parameterizations of the structure function, with a model of the EMC effect applied 
for Iron. 

With concerns about nuclear effects removed, then, there remain no impediments to studying duality 
for the first time in neutrino scattering with MINERi/A. Similarly, extractions of higher-twist contribu- 
tions and studies of evolution for parton distribution function extraction through the dependence of 
the structure functions will not be rendered ambiguous through the utilization of nuclear targets. 

Most strikingly, it appears that the nuclear effects at large x are the same in the resonance and scaling 
regimes, as evidenced by Figure l44lfrom ll62L where the cross-section ratios of carbon, iron, and gold 
to deuterium obtained in the resonance region (red) are the same as those obtained in the deep-inelastic 
regime (green). Whatever the underlying cause for medium modifications to the structure functions as 
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measured in nuclei, it is the same apparently for both hadronic (resonance) and scaling observables. 
Since the large x region of the EMC effect is ubiquitously attributed to Fermi motion in the nucleus, 
MINERzvA neutrino data should yield similar A-dependent results as the electron data in the figure. 
That is, it is expected (and will be tested) that the proposed data in the larger x and resonance regions 
will have the same EMC effect as data at higher W'^. Therefore, MINERi/A data at large x can be 
used for parton distribution studies, higher twist analyses, and nucleon structure studies with minimal 
nuclear extraction uncertainties. 
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Figure 44: Ratio of electron-nucleus scattering data (from top to bottom. Carbon, Iron, Gold) to that 
obtained from Deuterium scattering, for data in both the resonance (red) and deep inelastic (green) 
regimes. The data are plotted as a function of the Nachtmann scaling variable ^, allowing direct com- 
parison of high W^, DIS data to lower W"^, resonance data. 



The expected numbers of events for the resonance and deep inelastic regimes are tabulated in 
Section 15.41 These will make possible all of the studies here discussed in the perturbative and non- 
preturbative transition region of larger x and lower values. This is an exciting regime, with many 
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unanswered problems both interesting on their own and of import to other high energy applications. 
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11 Generalized Parton Distributions 



One of the main goals of subatomic physics is to understand the structure of hadrons, and in particular 
the structure of the nucleon. The primary approach to this problem has been through measurement of 
the nucleon form-factors, with (quasi-)elastic scattering (for up to a few (GeV/c)^), parton densi- 
ties, through inclusive deep-inelastic scattering (DIS), and distribution amplitudes, through exclusive 
processes. However, the usual parton densities extracted from DIS are only sensitive to the longitudinal 
component of the parton distributions and do not give information on the transverse component, or other 
contributions to the nucleon angular momentum. 

11.1 The Nucleon Spin Puzzle and GPDs 

In the late 1980's, results from polarized DIS showed that a relatively small fraction, about 20%, of 
the nucleon spin is carried by the valence quarks. The obvious candidates for the missing spin were 
the quark and gluon orbital momentum and gluon helicity. However, information on those quantities 
cannot be extracted from DIS. 

In 1997, Jil 163lll64l showed that a new class of nucleon observables, which he called "off-forward 
parton distributions", could be used to determine the spin structure of the nucleon. This work, along 
with developments by others, especially Radyuskinf 165 , 166| and Collins lll67l showed that these dis- 
tributions, now called generalized parton distributions (GPDs), had the potential to give a full three- 
dimensional picture of the nucleon structure. This exciting development has led to an immense amount 
of theoretical work in the last few years. Short reviews can be found in I168II1691 and a comprehensive 
review can be found in 11701 . 

Ji showed that in leading twist there are four GPDs, which he called H, H, E, and E, for each quark 
flavor. H and H are nucleon helicity-conserving amplitudes and E and E are helicity-flipping ampli- 
tudes. The GPDs are functions of x, ^ (a factor determining the "off-forwardness" of the reaction), and 
the total momentum-transfer squared, t. The GPDs can be accessed experimentally through reactions 
proceeding via the "handbag" diagram shown in Figure 145) 

11.2 Deeply-virtual Compton Scattering 

The most promising reaction to measure GPDs identified so far is deeply-virtual Compton scattering 
(DVCS). The DVCS reaction is shown in Figure l46h . An interesting feature of DVCS is that it can 
interfere with the Bethe-Heitler process. Figure l46b . which is completely calculable in terms of the 
nucleon elastic form-factors. This interference causes an asymmetry in the azimuthal distribution of 
the scattered proton allowing some quantities to be determined that would otherwise require a polarized 
target. However, DVCS involves a combination of the four GPD amplitudes, which cannot be separated 
using DVCS alone. Some complementary information can also be obtained from nucleon form-factor 
measurements and deep exclusive meson electroproduction. 

Neutrino scattering provides a very similar reaction to DVCS. In this case, the virtual exchange is 
of a with the production of an energetic photon, a with either a recoiling nucleon or nucleon 
resonance, as shown in Figure]^ This "weak DVCS" reaction is very promising theoretically because 
it provides access to different GPDs than DVCS. It will help resolve the individual flavors, e.g. d in 
neutrino scattering and u in anti-neutrino scattering, and the interference of the V and A currents will 
give access to C-odd combinations of GPDs. 
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Figure 47: Reactions sensitive to GPDs in neutrino scattering. 
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11.3 Measurement of GPDs in MINERi/A 



Studies of the weak DVCS reaction are currently underway by A. Radyuskin, A. Psaker, and W. Mel- 
nitchouk. One very encouraging result so far is that the qq equivalent of the polarized structure function 
gi can be measured without using polarized targets. This would allow separation of the valence and 
sea parts of the spin-dependent GPDs, and help determine the role of the axial anomaly in the proton 
spin puzzle. In addition, although the Bethe-Heitler process is suppressed, it is still present and the 
interference with it would allow measurement of individual GPDs. 

The estimated cross-section for weak DVCS is about 10^^^ cm^. For neutrino energies in the 5- 
10 GeV range, this would yield a few hundred events/year for a 1 GeV-wide bin in neutrino energy. 
Although most of the events will be from nucleons in carbon, any nuclear modifications are expected to 
be small except at very small or very large x. 

At least one other reaction, the hard exclusive production of the Ds has also been proposed ll71l 
as a probe of GPDs. This reaction is sensitive to the gluon structure of the nucleus. Unfortunately, the 
cross-section for this reaction (estimated at 10~^ pb for > 12 (GeV/c)^), is too small to be measured 
with precision in MINERz/A. Nevertheless, over the entire run perhaps a few hundred events would be 
observed over all values of Q^, which would give some information on the gluon GPDs. 
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12 Studying Nuclear Effects with Neutrinos 



In most neutrino scattering experiments, massive nuclear target/detectors are necessary to obtain use- 
ful reaction rates. Neutrino-oscillation experiments, despite the extremely intense beams designed for 
them, must also use very massive Iron, water or other nuclear target/detectors, since they are located 
hundreds of kilometers from the production point. Analysis of neutrino reactions within nuclear media 
requires an understanding of certain processes which are absent in neutrino scattering on free nucleons; 
these processes involve the so-called "nuclear effects". Two general categories of such effects can be 
distinguished. Effects comprising a first category include: 

• The target nucleon is moving within the nucleus and, when incoming neutrino energies are < 
2 GeV, the initial target energy and momentum must be accounted for using simulations which 
include either a target Fermi gas model or, preferably, nucleon spectral functions. 

• Certain final states are excluded as a result of Pauli blocking among identical nucleons. 

• The resulting final state may undergo final state interactions (FSI), including re-scattering and 
absorption; these interactions may significantly alter the observed final-state configuration and 
measured energy. 

The first two effects are either already included in Monte Carlos or are currently being examined 
in collaboration with nuclear theorists and will soon be included. The third effect is perhaps the most 
troublesome for current and future neutrino experiments. There is a dearth of data for which nuclear 
effects on specific hadronic final states (the fragmentation functions) have been isolated, whether for 
neutrino or charged-lepton beams. These effects are likely to be sizable for neutrino energies producing 
a large fraction of elastic and resonant final states l6l. 

A second category of nuclear effects are those by which the neutrino interaction probability on 
nuclei is modified relative to that for free nucleons. These effects occur across a wide range of neutrino 
energies and are normally categorized by the xbj of the quark involved in the scattering, and the 
of the intermediate vector boson exchanged. Nuclear effects of this type have been extensively studied 
in deep-inelastic scattering (DIS) measurements of structure functions using muon and electron beams. 
For example. Figure |48l shows the ratio of the structure function F2 measured on a heavy nuclear target 
to F2 measured for Deuterium. 

With neutrino beams, these nuclear effects have only been studied with low-statistics in bubble- 
chamber experiments. 

12.1 Final-state Interactions 

Distortions which result from FSI depend on the particle type and energy. Of primary concern are 
effects involving final-state nucleons and pions. For nucleons, rescattering is the major effect, resulting 
in i) change of direction and energy loss, ii) production of secondary nucleons, or Hi) neutron or proton 
pickup leading to deuteron emission. For pions, FSI can also lead to scattering with possible nucleon 
emission. The pions can charge exchange or be totally absorbed leading to emission of nucleons only. 
In all of these cases, particles that escape from the nucleus have lower energy than the initial, primary 
particle, and the redistributed energy information is lost due to detector thresholds. 
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Figure 48: The trend of the ratio of F2 measured with a heavy nuclear target to F2 measured using 
deuterium, for charged-lepton scattering, as a function of xbj- 



The most reliable information on FSl for nucleons comes from transparency measurements in 
(e, e'p) reactions on nuclei. Transparency, defined as the probability of escaping the nucleus with- 
out interaction, is measured by detecting the scattered electron and integrating the protons detected 
within the quasi-elastic peak. The most recent results quoted for protons in the energy range 0.5 to 
4 GeV are about 60% for C, to 40% for Pbl l72l . with very little energy-dependence. The composition 
and energy distribution of the final-state particles is not well measured. Two-proton and proton-neutron 
final states should dominate, with each nucleon having about half the total energy of the initial nucleon. 

Pion interactions, especially for pions below a few hundred MeV, are dominated by the A resonance. 
The data on FSl can be inferred from reactions of free pions on nuclei. There is little specific data for 
pions resulting from A or other resonant particle production in the nucleus. 

The significant feature of pion reactions is the strong absorption component - both inelastic scat- 
ters and "true" absorption, when the pion disappears from the final state. The absorption component 
comprises about two-thirds of the total cross-section. True absoprtion ranges from about 25% of the 
total cross-section for C to nearly 40% for Pb, for both positive and negative pions in the 100-300 MeV 
range[173 1. Inelastic cross-sections are comparable, although generally smaller for heavier nuclei. Be- 
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cause of the strong absorption component, pions in this energy range should escape the nucleus only 
about 50% of the time. 

Several experiments have found pion absorption to be a fairly complicated process in complex 
nucleif 1 741 |175.I . Although the first step is believed to be absorption on an isospin-0 np pair (quasi- 
deuteron), even in a nucleus as light as carbon the final state is dominated by three-nucleon emission. 
For heavier nuclei, the final state has a large component of four-nucleon emission. Even 7r+ absorption 
usually includes emission of a neutron, and of course vr" absorption is dominated by neutron emission. 

There is very little information on pion reaction cross-sections for energies above about 500 MeV. 
The total pion-nucleon cross-section drops significantly from the 200 mb resonant peak at 200 MeV 
to around 30 mb for energies above 500 MeV. Since this is not significantly different than the nucleon- 
nucleon cross-section, pion transparency should be comparable to proton transparency at higher ener- 
gies, i.e. approximately half the pions will react through either scattering or absorption. 

Interactions of 1-10 GeV neutrinos will produce pions with a wide range of energies. It should be 
noted that backward decay of the A resonance can produce rather low-energy pions, because the decay 
pions have a velocity in the A rest frame comparable to that of the A in the lab. 

The large absorption cross section (100-200 mb for C, 400-600 mb in Fe) for 100-300 MeV pions 
means that even pions that escape the nucleus may interact again, with absorption rates of a percent/cm 
in scintillator. 

There are other effects which influence the observed transparency of produced secondaries. As 
described in the Nuclear Effects section of Chapter the quantum effect of hadron formation length 
and the QCD effect of color transparency can increase the probability that a secondary escapes the 
nucleus without undergoing FSI. These effects are proportional to the energy and transfer and will 
not influence the transparency of low momentun secondaries. 

12.2 Nuclear Effects and Interaction Probabilities 

MINERz^A will provide the setting for a systematic, precision study of cross-sections and, with suf- 
ficient V, structure functions, on a variety of nuclear targets. Briefly reviewing the nuclear effects on 
measured structure functions (directly proportional to the cross-sections) as a function of xbj reveals: 

12.2.1 Low-x: Nuclear shadowing 

In the shadowing region, x < 0.1, there are several areas where neutrino scattering can provide new 
insights compared to charged-lepton probes. "Shadowing" is a phenomenon which occurs in nuclear 
targets and is characterized by reduction of the cross-section per nucleon for larger-A nuclei, such as 
Fe, compared to smaller- A nuclei such as D2. A recent summary of theoretical interpretation of this 
effect is presented in [ 176J . 

Vector-meson dominance (VMD) is the accepted explanation for shadowing in the scattering of 
charged leptons off nuclei (i.e. fi/e—A)forQ'^ <~ 5GeV^. In essence, the incoming boson dissociates 
into a qq pair which interacts strongly with the nucleus as a meson. Due to the V-A nature of the weak 
interaction, neutrino scattering should involve not only a VMD effect but additional contributions from 
axial- vector mesons such as the oi. Other sources of nuclear shadowing (mainly in larger nuclei) 
involve gluon recombination from neighbors of the struck nucleon, shifting the parton distributions 
toward higher values of x. At higher Q^, shadowing is dominated by Pomeron exchange in diffractive 
scattering. 
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A quantitative analysis of neutrino shadowing effects by Kulagin ll77ll uses a non-perturbative par- 
ton model to predict shadowing effects in v-K scattering. As illustrated in Figure 0^ which predicts 
the ratio of scattering off Fe to scattering off D2, shadowing effects with neutrinos should be dramatic at 
low (the kinematic region of the NuMI neutrino beam) and still significant even at large . Kulagin 
also attempts to determine the quark-flavor dependence of shadowing effects by separately predicting 
the shadowing observed in F2(x, Q^) (sum of all quarks) and xF3(x, Q^) (valance quarks only). These 
predictions should be testable in MINERz^A. 



Nuclear shadowing effect in Fe at = 0.7 GeV^ 
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Figure 49: Expected shadowing effects off an Fe target at = 0.7 (GeV/c)^ with Kulagin's non- 
perturbative parton model emphasizing the difference in shadowing for F2 and XF3. The arrows in the 
vicinity of R = 0.8 indicate the expected shadowing strength at = 15 (GeV/c)^. 



12.2.2 Mid-x: Anti-shadowing and the EMC effect 

Drell-Yan experiments have also measured nuclear effects and their results are quite similar to DIS 
experiments in the shadowing region. However, in the anti-shadowing region where Ra, the ratio of 
scattering off a nucleus A to scattering off Deuterium, makes a statistically-significant excursion above 
1.0 in DIS, Drell-Yan experiments see no effect. This could indicate a difference in nuclear effects 
between valence and sea quarks as also predicted by Kulagin. 
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Eskola et al. lll78l have quantified this difference using a model which predicts that the differences 
between nuclear effects in xF^{x, Q^) and ^2(2;, Q^), identified by Kulagin in the shadowing region, 
should persist through the anti-shadowing region as well. More recent work by Kumano ll79l supports 
these findings using different fitting techniques. 

Based on the various theoretical explanations for the anti-shadowing and EMC effects existing 
today, the measured effects could be considerably different for neutrinos. Neutrino scattering results 
would help clarify the theoretical understanding of this phenomenon. 

12.2.3 High-x: Multi-quark cluster effects 

Analyses from DIS experiments of F2{x, Q^) in the "Fermi-motion" region x > 0.7 have used few- 
nucleon correlation and multi-quark cluster models to fit the data. These models boost the momentum 
of some quarks, which translates into a high-x tail of F2{x, Q^) that should behave as e"'*^. However, 
fits to ;U — C| 180 1 and u — Fe llSlI scattering have obtained two different values for the fitted constant 
a: a = 16.5±0.5 and a = 8.3±0.7ib0.7 (systematic), respectively. This is surprising because any few- 
nucleon-correlation or multi-quark effects should have already saturated by Carbon. A high-statistics 
data sample, off several nuclear targets, could go a long way towards resolving the dependence of the 
value of a on the nucleus and lepton probe. 

12.3 Measuring Nuclear Effects in MINERz/A 

To study nuclear effects in MINERz^A, Fe and Pb nuclear targets will be installed upstream of the pure 
scintillator active detector which, essentially, acts as a carbon target. Two configurations are currently 
being examined. One would have (upstream to downstream) three 2.5 cm Fe plates, each plate followed 
by a module of active scintillator detector. Following this would be six 0.8 cm Pb plates (equal radiation 
thickness to the Fe) again separated by scintillator modules. This would give just over 1 ton of each 
target. The second possible configuration involves a total of six planes only, with each plane divided 
transversely into Fe and Pb segments. As one proceeds upstream to downstream, the Fe and Pb ex- 
change sides on each of the six planes. As always, a scintillator module separates each of the six planes. 
This configuration would also translate to just over 1 ton of each target. For the standard four-year run 
described in Section l53l MINERi^A would collect 940 K events on Fe and Pb and 2.8 M events on the 
C within the fiducial volume of the scintillator. 

MINERz^A's goals in measuring nuclear effects can be summarized as follows: 

• Measure final-state multiphcities, and hence absorption probabilities, as a function of A with 
incoming v; 

• measure the visible hadron energy distribution as a function of target to determine relative energy 
loss due to FSI; 

• measure a{xBj) for each nuclear target to compare x^j -dependent nuclear effects with u and 
charged lepton. 

• With sufficient 17, measure the nuclear effects on F2(x,Q^) and xF3(x,Q^) to determine whether 
sea and valence quarks are affected differently. 
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12.3.1 Multiplicities and visible hadron energy 



The expected average multiplicity of neutrino events as a function of Eh, with no nuclear effects, is 
shown in Figure I^U] As mentioned earlier, FSI will perturb this distribution via pion absorption and 
hadron re-scattering in the nuclear medium. FSI will also distort the initial hadron energy, transfered by 
the intermediate vector boson, yielding less visible energy in the detector. Restricting the study to events 
where all particles stop within the 2 m of active scintillator downstream of the nuclear targets will permit 
measurement of the hadron energy by range to within a few percent. The sample of events meeting these 
criteria is a function of the hadron energy Eh, and is shown in Figure|^ As can be seen, even at higher 
values of the hadron energy u, around 20% of the events have all secondary tracks contained within the 
active scintillator volume. With nearly one million events on each nuclear target in the four-year run, 
there will be sufficient statistics to determine the nuclear dependence of both multiplicities and visible 
hadronic energy. 
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Figure 50: The average multiplicty, excluding neutrons, as a function of the hadron energy of the event. 
The distribution is predicted by the NEUGEN Monte Carlo without turning on FSI. 



12.3.2 X Bj -dependent nuclear effects 

Just over 16% of the total event sample has xbj < 0.10. The (approximate) statistical accuracy for 
measurements of the nuclear effects in the ratios of Fe to C events at small x (shadowing region) are 
summarized in the following table. The columns designated DIS indicate that a cut has been made to 
retain only events with > 2.0 GeV/c^ and > 1.0 (GeV/c)^. For the MINERi/A DIS analysis. 
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Figure 51: The fraction of events which are fully contained within the active scintillator detector down- 
stream of the nuclear targets as a function of the total energy of the hadronic system. The distributions 
are for quasi-elastic, resonant, DIS and all reactions, as noted. 



the first three bins could be combined into two bins to reduce statistical errors. 

Assuming the level of shadowing predicted by Kulagin, the measured ratio of Fe/C and Pb/C, with 
statistical errors corresponding to the data accumulated during the 4-year run, is shown in Figure |52l 
The ratios plotted are for all events. The statistical errors would increase, as indicated in Table fTTl after 
making a DIS cut. 

The baseline 4-year run would be adequate to achieve the physics goals of the nuclear effects study, 
although some would be limited by the kinematic reach of the neutrino beam energies used for MINOS 
running and the minimal V exposure planned for MINOS. 

12.4 Nuclear Effects and Determination of sin^ dw 

There have been many attempts to explain the recent NuTeV I182II measurement of sin^ Qw, which is 
3(T away from the Standard Model prediction. Among the most persuasive are the unknown nuclear 
corrections involving neutrinos ll83l . MINERi/A will be able to directly measure the ratio NC/CC on 
various nuclear targets to explore these nuclear effects experimentally. 



105 



1.0 



— f- 

I 



0.9 



^ 0.8 



Pi 



0.7- 



0.6- 



0.5 



.001 



.01 



Fe/C 



Pb/C 



Figure 52: Kulagin's predicted ratio of shadowing effects off Pb, Fe and C targets with the expected 
errors from all events from the 4-year run. 
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Table 11: Statistical errors on the ratio of fully-contained Iron to Carbon events, assuming the level of 
shadowing predicted by the model of Kulagin, as a function of the xbj bin. 
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13 MINERi/A and Oscillation Measurements 



Over the past decade neutrinos have moved to center stage in the field of particle physics with the 
discovery of neutrino oscillation. Following on the initial discovery of solar and atmospheric neutrino 
oscillation are a new generation of high-precision long-baseline experiments dedicated to mapping out 
the neutrino mixing matrix and mass hierarchy in detail. In this section we address some of the ways 
in which the measurements made by MINERz/A can help these ambitious and expensive experiments 
achieve maximum senstivity. 



13.1 Neutrino Oscillation Landscape 

One accelerator-based experiment to explore the atmospheric oscillation sector has already begun, and 
several more are in the construction phase. The K2K experiment in Japan has seen evidence for dis- 
appearance, and expects to double its sample of about 50 events over the next yearfTH^- The MINOS 
experiment, with a much larger expected event sample, will make the first precision measurement of the 
atmospheric mass splitting, again through z/^ disappearance l204Tl . Finally, the OPERA and ICARUS 
experiments in Europe will attempt to further confirm the z/^ oscillation hypothesis by recon- 

structing actual Vr charged-current interactions in a beam produced as f^. The solar sector is being 
addressed by novel detection techniques of solar neutrinos themselves, and the KamLAND experiment, 
which uses a number of reactors as its anti-neutrino source 1 186 1 . If confirmed by MiniBooNE| 185 1, the 
LSND anomaly would dramatically affect the lines of inquiry for future experiments, demanding precise 
oscillation measurements with both long and short baseUnes, and hence both 1 GeV and several-GeV 
neutrino beams. 

One reason for the flurry of recent activity in neutrino physics is that non-zero neutrino masses and 
mixing have profound implications not only for the origin of flavor in the universe, but possibly also 
the origin of the matter-antimatter asymmetry. Because the lepton mixing matrix seems to have large 
off-diagonal elements, leptonic CP violation could be much larger than observed in the quark sector, 
and may be large enough to explain the matter/anti-matter asymmetry that we see today. 

A three-generation neutrino mixing matrix can be described by three independent mixing angles 
(^'i2i ^23i ^13) and a CP-violating phase (5cp)- The standard notation for this matrix, which transforms 
between the flavor and mass eigenstates is as follows: 
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where if Sij = sin 6*4^ , Cjj = cos Qij, then \J can be expressed as three rotation matrices: 



( I \ / C13 sne'^ \ ( C12 S12 

C23 S23 1 -Si2 C12 

V -S23 C23 / V -sise-^^ ci3 y V 1 



In this parameterization, the Vr oscillation probability, which describes atmospheric neutrino 
disappearance, can be expressed: 



P{y^^ ^v^ = \- cos^ 6*13 sin^ 26*23 sin^ ^ 



Am23L 
4S 



(58) 
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The solar (electron) neutrino disappearance, which has been further confirmed by the KamLand 
reactor (electron anti-)neutrino experiment (with average baseline 100 km), can be expressed as: 

P(i/e ^ z^e) = 1 - Sin' 29u sin' (^^J^) (59) 

The measurements in the solar and atmospheric sectors have shown that the mixing angles Ou and 
023 are large, but there remains one mixing angle which has not been determined, ^13. This angle would 
be manifest by electron neutrino disappearance a few kilometers from a reactor, or electron neutrino 
appearance in a few-GeV i/^ beam a several hundred kilometers from an accelerator. In the latter case 
the oscillation probability in vacuum is 

Pil^^. -^ve) = sin' ^23 sin' 2^13 sin' |^ ^^|'^ ^ + ... (60) 

where the missing terms are suppressed by at least one factor of Amf2/ Am23. 

Although reactors can play an important role in discovering non-zero ^13, this field will rely on 
accelerator experiments to eventually search for CP violation and determine the mass hierarchy. For 
example, the asymmetry in neutrino and anti-neutrino oscillation probabilities, in the absence of matter 
effects, is (to first order): 

P{1^^l Ve) - Pji'^i i>e) _ Aml2L sin 5 ^ ^ ^^^^ 
P{u^ Ue) + P{v^ i>e) E sin 6*13 

When electron neutrinos pass through the earth they can scatter off electrons, which creates an 
additional potential not present for muon or tau neutrinos [187]. This additional potential means the 
effective mixing angle and oscillation length is changed from equation |60l and is changed differently 
for neutrinos and antineutrinos. Moreover, the sign of the asymmetry is determined by whether Am23 
is positive or negative. The asymmetry in v^. oscillation probabilities due to matter effects, in the 
limit of Amf 2 being zero, is (to first order, when E < Er) 



P{v^ Ue) + P{l>^ Ve) 




E Eom 

E 



(62) 



En = ^^^llGeV 



E 



Am'^L 



2tt 

To measure this asymmetry, oscillation experiments will need to search for electron neutrino ap- 
pearance in muon neutrino (and anti-neutrino) beams, and to measure the atmospheric Am^ precisely, 
future (and current) experiments will need to measure the muon neutrino survival probability with cor- 
responding precision. Both kinds of experiments will require extremely long baselines, as well as near 
and far detectors to make the actual probability measurements. Even with an identical near detector, 
oscillation measurements will require reliable neutrino interaction models. For disappearance mea- 
surements, these models will be used to determine the mixing parameters from measured distributions 
in near and far detectors. For U(, appearance measurements, these models will be used to predict the far 
detector backgrounds based on data from a near detector. In both cases, the measurements are compli- 
cated by the fact that the far detector's z/^ charged-current event spectrum is dramatically different from 
the near detector's, due to the large, energy-dependent disappearance probability. 
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13.2 Benefits of MINERi/A to Oscillation Experiments 



With its fine grained, fully-active inner detector, electromagnetic and hadronic calorimetry, and excel- 
lent muon measurement capabilities, MINERi^A will have superb pattern recognition, energy resolu- 
tion, and efficiency. These abilities, coupled with the high flux of the NuMl beam make possible a 
host of improvements that will directly assist oscillation experiments. Indeed, most of the physics areas 
discussed in this proposal, from the precision measurement of quasi-elastic form factors, to exclusive 
channel studies, to coherent production, to nuclear physics, have the potential to improve oscillation 
measurements in one way or another. In particular, improving our knowledge in these areas may help 
avoid problems in the analysis of oscillation data that are difficult to foresee at this point. 

In practice, this benefit will be realized in the development of improved neutrino event generators 
that encapsulate the information learned from MINERi^A and provide a powerful and portable resource 
for all future neutrino experiments. The primary authors of two of the most widely-used, publicly 
available event generators are actively involved in MlNERz^A and see the development of such a 'next 
generation' generator as one of the principal tangible benefits of this experiment. 

The remainder of this chapter will focus on two specific areas where MINERz^A can aid oscillation 
experiments. One is determination of the "neutrino energy calibration", important for Arn^^ measure- 
ments, and the second is measurement of backgrounds to appearance, in a search for ^13. 

While this section focuses on MINOS and the proposed NuMl off-axis experiment, MlNERi/A 
will undoubtedly benefit other future oscillation experiments as well, including the proposed J-PARCnu 
project in Japan. The J-PARCnu beam energy is matched to its shorter baseline (and is therefore lower), 
but neutrino energy calibration and neutral-current vr'^ backgrounds are as essential to J-PARCnu as 
they are at NuMl. Neutrino energy reconstruction in J-PARCnu (as in K2K) is limited by knowledge of 
the non quasi-elastic background induced by inelastic reactions which feed down from higher neutrino 
energies. Similarly, most neutral-current background to appearance in J-PARCnu originates from the 
high-energy tail of the beam. Although MlNERz^A obviously cannot directly measure the J-PARCnu 
beam (as it can for NuMl), its somewhat higher energy reach arguably makes it better-suited for mini- 
mizing J-PARCnu systematic uncertainties from neutrino-interaction physics than the J-PARCnu beam 
itself. 

13.3 Am23 Measurements 

As an example of the importance of neutrino interaction physics to oscillation experiments, consider 
a measurement of the "atmospheric" mass splitting, which is the primary (but not only) goal of the 
MINOS experiment. 

In a long-baseline experiment there are two limiting cases: a far detector without a near detector 
(and maximum systematic uncertainties) and a far detector with a perfect near detector (and negligible 
systematic uncertainties). For a far detector without a near detector, the predicted "no oscillation" 
distributions are determined by integrating the flux and cross section over a smearing function with 
takes into account detector acceptance, reconstruction inefficiencies and measurement resolution. In this 
case the determination of oscillation parameters can only be as good as the understanding of the beam, 
neuttino interaction cross-sections, detector performance, and reconstruction. In the other idealized 
extreme, identical near and far detectors see an identical spectrum. In this case, allsources of systematic 
uncertainty cancel in the near/far comparison. 

For MINOS, the near detector will help reduce many of the important systematic errors, but the 
situation is not quite as perfect as ideal case. The beams (oscillated vs. unoscillated and point source 
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vs. line source), detector shapes (octagon vs. "squashed octagon"), photodetectors (Hammamatsu Ml 6 
vs. M64), electronics (IDE vs. QIE), reconstruction (single neutrino event per readout vs. multiple 
neutrino events per readout), and beam-related backgrounds all differ between far and near. Great effort 
in the areas of calibration, reconstruction, and analysis is invested to understand and correct for these 
small near-far differences. In addition, dedicated measurements of hadron production on the NuMI 
target by the MIPP experiment should reduce the uncertainties on the near to far flux ratio to the few 
percent level. 

13.4 Neutrino Energy Calibration 

Analysis of near detector data aims to predict expected rates and spectra in the far detector in the ab- 
sence of oscillation, and for different values of ^23 and Am^a. Differences between these predictions 
and real data will be used to fit the oscillation parameters. A crucial link in this prodcedure is the 
translation between measured energy and the neutrino energy, as this quantity is directly related to the 
oscillation probability. Even the most elaborate suite of near- and far-detector calibrations can at best 
characterize the response of the detector to known incoming charged particles. There is no equivalent 
test beam for "neutrino calibration"; this final step requires appeal to a model. A reliable understand- 
ing of the spectrum and multiplicity of particles produced in neutrino interactions is indispensible for 
reconstructing the true neutrino energy from the visible energy measured in a calorimeter. There are 
several areas where these models have large uncertainties, which MINERi/A could help to reduce: 

13.4.1 Charged and neutral pion production 

As MINOS responds differently to electromagnetic and hadronic showers, it will be essential to estimate 
the relative abundances of charged and neutral particles. These abundances are determined by isospin 
amplitudes at each point in phase space. As explained in Section 115.11 these amplitudes necessarily 
include a resonant component from the low invariant mass region where specific exclusive channels are 
produced, and the deep-inelastic regime where scattering is described by quark transitions in the frame- 
work of the parton model. The predictions therefore depend on the resonance model adopted, the model 
for fragmentation of inclusive quark final-states (particularly at low invariant mass where standard mod- 
els like the JETSET string model are not applicable), and the treatment of the deep-inelastic/resonant 
overlap region. 

13.4.2 Charged particle multiplicities 

Neutrino energy reconstruction also depends on the charged particle multiplicity, as the rest energy of 
pions disappears via nuclear absorption and the neutrinos produced in vr and decays. Correction for 
these losses is therefore also related to the model(s) of charged pion production. 

13.4.3 Intra-nuclear scattering 

At NuMI energies, intra-nuclear scattering can result in large distortions of the hadronic multiplicities, 
angular distributions and total energies. A feeling for the variation of intranuclear rescattering can 
be gleaned from Figure |53j which shows the 7r+ spectra from 3 GeV neutrino interactions on three 
different target nuclei: Carbon, Iron, and Lead. Little data is available to constrain or validate models 
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of these effects. Measurement of these processes requires a Att detector, hence electron scattering data 
has limited utility. 
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Figure 53: Simulated 7r+ momenta for 100,000 3 GeV i/^ charged-current interactions on Carbon, Iron, 
and Lead. 



13.4.4 Expected results 

Figure IS^l a) shows the "neutrino energy resolution" for 1 GeV and 2 GeV charged-current interactions 
in a perfect calorimeter that measures the kinetic energy of charged pions and baryons but the total 
energy of photons and electrons. The structure in the distribution results from production of one or more 
charged pions, and with corresponding amounts of lost rest energy. Figure l54T b) shows the average ratio 
of visible to true neutrino energy versus neutrino energy for the same detector. For the solid lines there is 
no intranuclear rescattering, while for the dashed lines rescattering in a carbon target is assumed. Both 
models show that a correction as large as 10% is required, which translates into a Am23 uncertainty 
equivalent to the statistical error of MINOS [ 188 1. 

MINERz^A can play an important role in the MINOS Am^ analysis by measuring the charged-pion 
multiplicity as a function of visible neutrino energy. Data from a variety of nuclear targets (includ- 
ing Iron) will strongly constrain and redundantly validate rescattering models. Figure l55T a) shows the 
charged paticle multiplicity for charged-current interactions in the NuMI low-energy beam, for the 
two models described in Figure l54T b). Finally, Figure l55\ b) shows the difference between the oscil- 
lated and unoscillated interaction spectra for the two models of nuclear rescattering. Note that the 
differences below 1 GeV are enormous, but even at 3 or 4 GeV they are sizeable. 

To quantify the effect this would have on a Am^ measurement, a toy monte carlo was used to 
approximate the MINOS Am^ analysis, including neutral current contamination and cuts to reduce that 
contamination. Figure |56lshows the fractional size of the 90% confidence level contour region due to a 
20% uncertainty in the total "neutrino energy loss". Also shown in the figure is the size of the 90% CL 
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Figure 54: (a) Distribution of the ratio of visible to true neutrino energy for 1 GeV and 2 GeV charged- 
current neutrino interactions (b) the average of that ratio versus true neutrino energy, for two different 
models: one withoug intra-nuclear rescattering, and one where rescattering in Carbon is simulated. 



region at sin^ 26 = 1, versus Am^ from a more complete MINOS simulation ( 11881 ). Although this 
systematic error would be lower than the statistical error for the lowest exposure of protons on target, 
it is far from negligible, and dominates for values of Am^ below l.5eV^ for all exposures. Without 
MlNERzvA the error due to nuclear effects may be roughly this size, but with even a small amount of 
data on several targets MlNERi^A would be able to model this effect to much better than 10% of itself, 
making this effect negligible. Future Am^ measurements (such at those proposed at J-PARC or NuMI 
Off Axis), which expect to achieve la statistical errors closer to the 1% level rather than the 3% level, 
will have to understand this effect (among others) even more precisely. 



13.5 Measurements 

A longer-term goal in oscillation physics is to probe leptonic CP violation and the neutrino mass hi- 
erarchy by comparing measurements of — > and — > oscillation. These measurements 
are particularly challenging because of backgrounds. Conventional neutrino beams always include 
contamination from muon and kaon decays in the beamline. In addition, neutral-current interactions 
produce vr'^ which can be mistaken for electron appearance in the far detector. 

In MINOS, neutral-current backgrounds should be several times the intrinsic u^, contamination, 
and 9i3 sensitivity depends strongly on the assumed systematic uncertainties, as shown in Figure l57l 
MINOS is not optimized to separate neutral-current interactions from charged-current interactions, 
and the background for this search is roughly 75% deep-inelastic vr^ production, 15% intrinsic u^, and 
10% charged-current interactions! 189 1. 

Because MINOS can only hope to achieve, at best, a factor two improvement in sin^ 6*13 sensitivity 
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Figure 55: (a) Charged particle multiplicity for the two models of nuclear rescattering discussed in the 
text, (b) Ratio of oscillated and unoscillated event spectra for those two models, assuming Am^23 = 

2.5 X 10--'^ eV^ 




Figure 56: Fractional size of the 90% confidence level contour at sin^ 26 = 1, due to MINOS statistical 
errors for different exposures, and due to a 20% uncertainty on "neutrino energy resolution" which 
would come about due to nuclear rescattering effects. 
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Figure 57: MINOS 3cr sensitivity to non-zero ^13 (assuming Am23 = 2.5 x 10^^ eV^), as a function 
of the systematic uncertainty on the background, for several possible integrated proton luminosities. 



over existing limits, other experiments better suited to this search have been proposed. One proposed 
experiment would use a much finer-grained detector off-axis from the NuMl beamline, where the neu- 
trino spectrum is much narrower and thus the signal to background ratio is higher. In this off-axis 
experiment, the dominant background would be intrinsic Ve, although neutral-current background is 
comparable, and the charged-currents are considerably smaller 1 190|. 

Since both MINOS and the off-axis experiment need to measure the intrinsic z/g and neutral-current 
background rates, they each require near detectors similar to their respective far detectors. Measur- 
ing neutral-current backgrounds directly is difficult at a near detector location, as charged-current 
interactions are far more abundant in the unoscillated beam. 

13.6 Ve Appearance Backgrounds 

MlNERi/A can measure all three types of Vf. appearance backgrounds. Thanks to superior segmentation, 
MlNERi/A can isolate a very clean charged-current sample and directly measure the v^. flux. Simi- 
larly, with its excellent vr*^ identification, energy {oe = Q%/\fE), and angular resolution, MlNERi/A 
can map out all the processes that produce neutral pions. Finally, with a fully active detector and good 
timing resolution, charged-current backgrounds can also be identified in MlNERz^A by exploiting 
the delayed muon-decay signature that is unavailable to oscillation detectors. The remainder of this sec- 
tion sketches two possible analyses illustrating both MlNERz^A's potentially decisive ability to isolate 
appearance backgrounds, and its impressive resolution. 

13.6.1 Beam 

The cleanest signature for charged-current interactions in MINERz/A will be the presence of an 
electromagnetic shower originating near a proton track. Figure |58] shows the distance between the 
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electromagnetic shower origin and the true primary vertex for charged-current interactions and vr 
production. The figure also shows the length of the showers, measured in MINERi/A scintillator planes, 
or 1.75 cm of polystyrene. For neutral pions the length is from the beginning of the first showering 
photon to the end of the second one. 
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Figure 58: (a) The distance in centimeters between the neutrino vertex, which can be determined from a 
proton track, and the start of the most upstream electromagnetic shower, for both electrons and photons 
from neutral pions. (b) The shower length in units of scintillator planes, for electrons and neutral pions. 

The MIPP experiment ! Ill will reduce the uncertainty on the z^e/ flux ratio to roughly 5%, but to 
determine the true intrinsic v^. background, the uncertainty on the quasi-elastic (Section^ and resonant 
(Section cross-sections must also be taken into account. 

With a simple analysis that requires a proton track and an electromagnetic shower depositing over 
0.5 GeV in the detector and starting within 2 planes of the proton track, MINERi/A would collect 
roughly 1500 charged-current events per year in a 3-ton fiducial volume, with a neutral-current 
background about a third the size of the signal. Figure |59l shows the resulting energy spectra for 
the z/g signal and neutral-current background. Further cuts to remove events with an identifiable second 
photon cluster from 7r° decay could reduce this background even further. 



13.6.2 Neutral-current vr^ production 

Neutral-current vr*^ production can occur through a number of mechanisms - resonant production, co- 
herent production, and deep-inelastic scattering. Figure |60| shows a striking example of MINERz^A's 
response to coherent vr*^ production. 

Coherent vr*^ production is a dangerous appearance background, because the neutral pion is pro- 
duced along the direction of the incoming neutrino, and carries away most of the neutrino's energy. See 
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Figure 59: Variables to identify charged current events in MlNERz/A : (a)Distance between vertex 
(determined by proton track) and most upstream converted photon from 7r° decay, and (b)Energy in 
electromagnetic cluster, for charged-current and all neutral-current events, (c) True (hadronic) energy 
for charged-current (all neutral-current) events, after applying the simple cuts described in the text. 
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MINERVA TOP VIEW Run Even! 16lntTypeC0H 

CC/NC 1 Mech. nu-p 
Vertex ( 0.0. 0.0,1336.3) 
PNEU 14 (0.0000,0.0000.2.3S21, 2.3321 ) 
PLEP *** ( 0670, -.065.0,0.04Sa, ****** ) 
X 0.000 q2 0.000 y 0.000 w2 0.0000 




Figure 60: A simulated neutral-current coherent vr production event in MINERi/A. The position of the 
vrO decay vertex can be determined accurately by extrapolating the two photons backward. Notice that 
both photons pass through a number of planes before beginning to shower, distinguishing them from 
electrons. 



Section|8]for a complete discussion of this process. Its cross-section uncertainty is ±50% or worse, and 
it has not been measured accurately at 2 GeV, the relevant energy for the NuMI off-axis experiment. 

Production of A (and other nucleon resonances) is another mechanism for faking a v^. signal, since 
their decay products often include vr*^. Neutral pions from resonance decay are not as energetic or 
collinear as those produced coherently, but their angular distribution mimics that of the signal. Resonant 
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vr are particularly susceptible to final-state nuclear interaction and rescattering, which will be studied 
in detail by MINERi/A using charged-current reactions (see Sections 171 andfT^. 

As a proof-of-concept, a sample of neutral-current single-vr'^ events has been selected using simple 
cuts. For events with two well-separated electromagnetic clusters (Et^ = Ei + E2), each passing 
through at least six planes of the fully-active region, requiring Ej^/Efot > 90% and Efot — E^^ < 
100 MeV efficiently isolates a neutral-current vr^ sample, as shown in Figure |^ After these cuts, the 
contamination of fg and charged-current interactions (combined) is less than 1%. The resulting 
sample contains about 2400 neutral-current tt^ events per 3 ton-yr, of which half are resonant and half 
coherent. 

Coherent and resonant interactions can be cleanly separated by cutting on the tt^ angle to the beam 
direction, as shown in Figure |62j which also highlights MINERi^A's excellent vr'^ angular resolution. 
The overall efficiency for selecting coherent neutral-current is about 40%. 

Finally, some — > backgrounds in oscillation experiments will come from deep-inelastic scat- 
tering, although that sample is easily isolated from the other two processes in MINERi/A because of 
the high multiplicity. Since the mean hadron multiplicity in deep-inelastic scattering is large, and the 

angular distribution rather flat, this channel is less likely to contribute background to a u^. search 
than the other two. On the other hand the cross-section for deep-inelastic scattering is larger, even at 
Ey = 2 GeV than for either resonant or coherent production, and most deep-inelastic interactions in 
MINOS or the NuMI off-axis experiment will fall into the poorly-understood W ^ 2 GeV/c^ transition 
region at the border of resonant production (see Sections ITOt . and fragmentation of low-W hadronic 
systems is not well-modeled by existing simulations like PYTHIA. 
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Figure 61: Selection of neutral-current single-7r" production. The variables plotted are the fraction of 
visible energy carried by the 7r° candidate {ET^/Etot) and the residual energy Efot — E^^. The left- 
hand plots show the backgrounds from f^(top) and fe(bottom). The plot at top right shows the same 
distribution for true neutral-current vr*^ production, and the lower right shows the subset from coherent 
scattering. In the neutral-current plots, notice the dramatic concentration of the coherent vr*^ signal in a 
single bin, in the left-most comer of the graph. All samples shown are normalized to a 3 ton-yr exposure 
of MINERz/A. 
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Figure 62: Angular distribution of neutral-current single-7r° sample. The plot at left shows all events 
passing the cuts on ET^/Etot and Etot—ET^ described in the text, broken down into coherent and resonant 
reactions. The coherent sample is strongly forward-peaked. The plot at right is a close-up of the forward 
region comparing the true and reconstructed 7r° angular distributions from the beam direction. The 
distributions are nearly identical, highlighting the MINERz/A's excellent angular resolution. 
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Part III 

Project Description 
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14 The NuMI Near Experimental Hall 



The MINOS Near Detector Hallfl911 is a fully-outfitted experimental facility that can accomodate the 
MINERi/A experiment with a limited number of additions to the infrastructure. 

The detector hall is 45 m long, 9.5 m wide, 9.6 m high, with its upstream end just over 1 km from 
the NuMI target, at a depth of 106 m below grade. The MINOS detector will be installed towards the 
downstream end of the hall, leaving a free space upstream amounting to, roughly, a cylinder 26 m in 
length and 3 m in radius. The neutrino beam centerline will descend at a slope of 3.3°, and enter the 
MINOS detector at a height of 3 m from the hall floor. The hall has been excavated, and is currently 
being outfitted for the MINOS near detector, with beneficial occupancy expected in early March, 2004. 

Ground water is pumped from the NuMI/MINOS complex at a rate of approximately 320^00 gal- 
lons (1300-1600 1) per minute. The hall floors and walls may be damp in places, and a drip ceiling will 
need to be extended upstream of the MINOS detector to protect MINERi/A. The air will be held at a 
temperature of between 60° F and 70° F (15° C and 21° C), and 60% relative humidity. 

14.1 Utilities 

The MINOS Service Building on the surface houses the access shaft to the Near Detector Hall, and is 
the entry point for electrical, cooling, and data services to the hall. A 15-ton capacity crane, with a 
hook height of 18.5 feet (5.66 m), will be used to lower the 3.47 T MINOS detector planes to the hall. 
MINOS Detector planes will be moved within the hall using an overhead 15-ton crane, with 22 foot 
(6.7 m) hook height and a coverage along the beam axis of approximately 40 m. 

Quiet power to the hall is provided by a 750 KVA transformer at the surface, which branches to 
a 45 KVA transformer for the muon monitoring alcoves, and two 75 KVA transformers for the Near 
Detector hall. The power needs of the MINOS detector account for the capacity of the 4 panelboards 
served by the two 75 KVA transformers, so additional panelboards for MINERz^A will likely be needed. 
The current estimate for MINERz^A electronics and high voltage power is less than 5000 W. It appears 
that overall capacity for the additional load exists within the MINOS hall, but this needs to be verified 
in detail. 

MINERz^A's main non-quiet power need is for the magnet coils, with an estimated ohmic power 
loss of 30 kW. The MINOS magnet coil power supply will be served by a 480 V line with 400 A 
capacity, but will require less than 80 kW of power. This should leave ample capacity for the addition 
of a power supply for the MINERz^A coil on the same line. 

The heat sink for the MINOS LCW cooling circuit is the flux of ground water collected in the 
MINOS sump. The cooling is adequate for MINOS, with an output water temperature of 70°F. This 
should be sufficient to absorb the heat load of the MINERz^A magnet, but would likely be too warm 
to effectively cool the front end electronics. The relatively low heat load of the MINERzvA electronics 
would hkely be absorbed without problems by the MINOS hall air conditioning. 

14.2 Detector Placement 

MINERj/A will be placed with its downstream end 1.75 m upstream of MINOS. This will leave suffi- 
cient work space between the two detectors, and avoid interfering with the MINOS coil, which extends 
approximately 1.5 m upstream of MINOS, to the lower right in the view of Figure |63l To have the 
beam axis intersect the detector axis close to the center of the active plastic target, the lower vertex of 
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the MINERi^A detector would be placed 1. 10 m off the hall floor. The beam centerline would enter the 
detector at an elevation of 3.4 m from the floor (Fi gures IS^ and E5t . 

MINERz/A will impinge slightly on a "stay clear" egress space for the lower MINOS detector elec- 
tronics racks. This could be resolved by either raising MINERi/A by less than 10 cm, or by rearranging 
the layout of the upstream part of the MINOS electronics platform and stairs. 




Figure 63: View of the proposed MINERz/A detector, and the MINOS detector, looking dowstream. 



14.3 Impact on MINOS 

The impact on MINOS of the heat load and power consumption of the MINERi/A detector can be made 
negligible through relatively minor additions to the hall infrastructure. The presence of the detector in 
the neutrino beam will cause an increase in the rate of activity in the MINOS detector, particularly in 
the first 20 planes forming the MINOS veto region. With the current design of the MINERi/A detector, 
the expected event rate in the detector is w 1.4 CC events / 10^'^ POT. For a spill of 2.5 x 10^'^ POT 
this is 3.4 CC events plus an additional 1.0 NC event. Since, in addition, the vectors of all particles 
leaving MINERz^A with a trajectory heading towards MINOS, will be made available to MINOS when 
MINERi^A is taking data, this should be a managable situation. 
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Figure 64: Plan view of the MINERzvA detector (purple outline near top of figure). 
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Figure 65: Front view of the MINERzvA detector. 
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14.4 MARS Simulation of Radiation Flux 



The intense neutrino beam will create a fluence of other particles due to neutrino interactions in the 
rock surrounding the experimental hall. Several physics topics are sensitive to background interactions 
caused, particularly, by neutrons. A MARS 14-based model ll92l has been created to estimate non- 
neutrino background in the detector. The model includes the rock surrounding the experimental hall, 
and the MINERz^A detector located upstream of MINOS. Both the detectors are positioned on the 
NuMI axis. The MINERi/A detector is simulated as described in Section but with the magnetic 
field ignored. 

The muon neutrino energy distributions used in this simulation for the Low, Medium, and High 
Energy beam configurations are shown in Figures l66ti68] There is an admixture of the other types of 
neutrinos from tt", kaon and muon decays shown in (Figurel69l. 
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Figure 66: Muon neutrino energy distributions at 
the near detector hall within a distance of 5 m 
around the axis of NuMI. 



Figure 67: Muon neutrino energy distributions at 
the near detector hall within a distance of 1 m 
around the axis of NuMI. 



The MARS 14 neutrino interaction model tracks the energy and angle of final state neutrinos, hadrons, 
e^, and /i^ from neutrino interactions. These particles, along with the showers they initiate, are trans- 
ported through the user-defined geometry and energy deposition and dose are estimated. MARS dis- 
tinguishes four types of neutrinos: u^, V^, and Vf., all of which are present in the NuMI beam 
(Figure l69t. The interactions included in the model are listed in Table fT2l 

The model is described in detail in 1 196 1. Some notable features of the model include: 

• Recoil of the target is simulated only for elastic interactions. 

• The model does not include inelastic neutrino interactions which produce pions via resonances. The 
cross sections for such processes are relatively small however compared to ones for the deep-inelastic 
(DIS) and coherent elastic scattering (process 8 in Table fT2b. 

• For charged-current DIS (process 1 in Table fT2b. the process of hadronization is simplified. Once 
the momentum of lepton is decided, the total momentum is balanced by a single pion, which is forced 
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Figure 68: Muon neutrino radial distributions at Figure 69: Neutrino beam components, 

the near detector hall. The origin corresponds to 
the NuMI axis. 



to undergo a deep-inelastic interaction in the same nucleus. This coarse "hadronization" is justifiable 
since we are interested in certain gross averages over the showers. 



14.5 Fluxes in MINERz/A 

Particle fluxes in the scintillator part of detector were calculated for three beam configurations and 
various threshold energies (Table [T^. All the fluxes are given for one spill with the beam intensity of 
2 X 10^'^ protons on target/spill. The units are lO^^cm^^. The total integrated path-length of a given 
type of particle is obtained by multiplying the flux by the fiducial volume of the detector. The columns 
"XX tot" and "XX sig" refer to particles coming from "all" sources or from only the central volume of 
scintillator respectively (signal). 
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Table 12: MARS model interactions for i/f and ui, where £ = {e, /x}. 



Particle 


(MeV) 


flux (IQ-^cm-^) 


LEtot 


LE sig 


sMEtot 


sME sig 


sHE tot 


sHE sig 


n 


0.1 


3.10 


0.0465 


6.50 


0.0986 


8.77 


0.1351 




20 


1.26 


0.0271 


2.45 


0.0586 


3.48 


0.0822 




100 


0.49 


0.0140 


0.95 


0.0258 


1.35 


0.0355 


Charged 


0.1 


0.66 


0.0523 


1.39 


0.1059 


2.04 


0.1534 


hadrons 


20 


0.73 


0.0522 


1.34 


0.1012 


1.93 


0.1550 




100 


0.55 


0.0505 


1.10 


0.0877 


1.57 


0.1220 


7 


0.1 


15.94 


0.3150 


33.64 


0.6033 


46.26 


0.9879 




20 


1.08 


0.0583 


1.63 


0.1254 


2.20 


0.1608 




100 


0.26 


0.0238 


0.41 


0.0493 


0.51 


0.0710 


e± 


0.1 


1.28 


0.0614 


2.11 


0.1272 


2.93 


0.1554 




20 


0.44 


0.0436 


0.83 


0.0717 


1.16 


0.1038 




100 


0.16 


0.0163 


0.26 


0.0318 


0.32 


0.0480 


/^^ 


0.1 


1.43 


0.0206 


2.61 


0.0416 


3.08 


0.0493 




20 


1.41 


0.0206 


2.58 


0.0417 


3.09 


0.0491 




100 


1.40 


0.0190 


2.59 


0.0397 


3.07 


0.0472 



Table 13: Particle fluxes averaged over the active taxget.Eth is a threshold kinetic energy used in the 
simulations. 
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15 Monte Carlo Studies and Performance 



This section outlines the event simulation and reconstruction software used to optimize the detector's 
design and quantify its physics capabilities. Much of this software has been borrowed from other ex- 
periments, where it has been thoroughly validated. The detector simulation and reconstruction software 
has been developed specifically for MINERz^A, but is based on widely-used libraries and algorithms. 

15.1 Event Generators 

The MINERi^A simulation software interfaces with two event generators that model neutrino inter- 
actions with matter: NEUGEN[120| and NUANCE(5H1|. NEUGEN was originally designed for the 
Soudan 2 experiment and is now the primary neutrino generator for the MINOS experiment. NU- 
ANCE was developed for the 1MB experiment and is currently used by the Super-Kamiokande, K2K, 
MiniBooNE and SNO collaborations. Both have evolved from "proprietary" programs designed for 
atmospheric neutrino studies into freely-available, general-purpose utilities that aim to model neutrino 
scattering over a wide range of energies and for different nuclear targets. Total charged-current cross- 
sections calculated by NUANCE (FigureEJ and NEUGEN (Figure|51l appear elsewhere in this proposal. 
As the results of the two generators agree with each other (to within the depressingly large range of un- 
certainties in available data) II 081 . they have been used interchangeably for the present studies. 

As in the past, future studies of neutrino oscillation and searches for nucleon decay will rely heavily 
on the best possible description of neutrino interactions with matter. Neutrino event generators are tools 
which encapsulate our understanding of this physics in an easily usable and portable form. Practically, 
they serve two related functions: to allow the rates of different reactions with the experimental target to 
be calculated, by providing total exclusive and inclusive cross-sections, and to simulate the dynamics 
of individual scattering events, by sampling the differential cross-sections. Many comparable packages 
are available to the collider physics community, and have been incrementally improved for decades, 
forming a common basis for discussion of different models and phenomena. One important goal of 
MINERi/A is to improve the quality of neutrino Monte Carlo event generators, and thereby enhance the 
physics reach of many future experiments. 

MINERz^A will attack this problem from both experimental and theoretical directions. Experimen- 
tally, MINERz^A will make definitive measurements of dozens of exclusive and inclusive cross-sections, 
across the range of energies most important for future oscillation and nucleon-decay experiments, with 
a well-controlled flux, and on a variety of nuclear targets. The era of 25% uncertainties and marginally- 
consistent cross-section data for even the simplest neutrino reactions will end with MINERz^A; for the 
first time it will be possible to validate the details, and not merely the gross features, of competing 
models. 

At the same time, MINERz^A will be a natural focus of attention for theorists and phenomenologists 
developing these models. NEUGEN and NUANCE are two of the most sophisticated neutrino-physics 
simulations in the world, but NUANCE models quasi-elastic scattering with the 1972 calculation of 
Smith and Moniz[571, and both programs use the Rein-Sehgal|4| resonant production model which 
dates from 1981. That no other widely-accepted models for these, the most fundamental neutrino- 
nucleon reactions, have emerged in the last quarter century is sobering evidence that an experiment like 
MINERz/A is long overdue. New, high-quality data is the surest way to catalyze theoretical ingenuity, 
and MINERi/A will provide the former in abundance. Through our contacts with these theorists, and 
ability to translate well-tested, state-of-the-art models into universally-available and widely-adopted 
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software, MINERi/A will serve as a conduit for expertise from a diverse collection of disciplines into 
the high-energy neutrino physics community.^ 

15.2 Beam Simulation 

The neutrino fluxes used in the simulation are derived from the GNuMlfT93l program developed for 
the MINOS experiment. GNuMI is a full GEANT| 194|-based simulation of the NuMI beamhne. As 
with all current neutrino beams, neutrinos arise from decay of vr, K and /i mesons that originate in col- 
lisions of a proton beam on a production target. GNuMI simulates all aspects of the neutrino beamline. 
Protons are fired into the target and the interaction products are transported through the focusing and 
filtering elements of the beam. Appropriate care is taken to ensure that the description of the beamline's 
geometry is as complete as possible, and that meson decays proceed with the correct kinematics and 
branching ratios. For this proposal, all fluxes in MINERzvA are taken from the official tables used by 
the MINOS collaboration. 

15.3 Detector Simulation 

The simulation of neutrino interactions in the MINERvA detector is carried out by a GEANT-based 
Monte Carlo program. This program combines a flexible description of the detector geometry, the 
NuMI neutrino beam flux from the beam simulation, neutrino interaction physics from either of the two 
generators and simulation of the scintillator response with the standard tracking and particle interaction 
routines available in GEANT. 

15.3.1 Interface to the GNuMI flux 

The output of the GNuMI simulation of the beamline is a set of files recording the neutrino flux in 0.5 
GeV bins for a nominal number of protons on target. The flux files are in a standard format and hence 
can be interchanged with no additional modifications to the code. In this way different beam configu- 
rations can be easily studied. An option exists to generate interactions with a flat energy spectrum. In 
this case, beam weights are stored in an output ntuple. This is particularly useful if one wishes to study 
the effect of different beam configurations without furthur Monte Carlo running. 

15.3.2 Interface to the event generators 

The Monte Carlo simulation program can be configured to accept neutrino interactions from either 
NEUGEN3 or NUANCE. The results of a neutrino interaction can be passed to the simulation in a 
number of ways. By default, the event generation routines in NEUGEN3 are usually called from within 
the simulation itself. In this mode, the code chooses a neutrino energy from the flux files, samples the 
density of material along the neutrino path; chooses a vertex and nucleus type, calls the kinematics 
generator and inserts the Ust of particles thus obtained into the GEANT data structures. This is not 
the only mode of generation. As a stand-alone generator, NUANCE provides events in either a text or 
ntuple format and so provision is made to read in events from a standard external format. NEUGEN3 

"This trend is already beginning, tlianks to collaborative work sparked by the NUINT series of workshops. The BBA- 
2003 quasi-elastic form-factor fits (see Chapter|6l and Bodek-Yang duality-inspired model of deep-inelastic scattering (Sec- 
tion ll0.3| have recently been implemented in NUANCE, and NEUGEN is exploring Benhar's spectral-function approach ll97l 
to nuclear binding effects. 
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has been modified to write out events in tlie same format, so that the results of both generators may be 
compared in a consistent manner. 

15.3.3 Geometry 

Flexibility drives the design of the detector geometry code. The size, segmentation, material and shape 
of all components of the detector can be set and altered almost entirely from the input datacards. The 
detector is logically divided into longitudinal sections. Each section can have different dimensions, 
strip sizes and absorber widths. In addition the absorbers in each section can be be constructed from 
segments of differing material and widths. The geometry description is sufficiently abstract that minor 
changes in detector design may be accommodated merely by changing the datacard, allowing for fast 
detector reconfiguration and easy bookkeeping. 

15.3.4 Hits and digitizations 

Particles are tracked through the GEANT geometry in the standard manner. When a particle traverses 
a sensitive detector volume the particle type, volume identifier, entrance and exit points and energy 
deposition (including Landau and other fluctuations) are recorded as a hit. When GEANT has finished 
tracking the event, the hits are considered and converted to digitizations. There are as many digitizations 
as there are strips hit. Multiple hits on a single strip are condensed into one digitization, although 
information on which tracks contributed to the digitization is stored. These digitizations are then passed 
to the event reconstruction program. 

15.4 Photon Transport Simulations 

The MINERi^A detector simulation assumes "ideal" light collection, and records the raw energy de- 
posited in each channel. During event reconstruction, the energy deposit is converted to a number 
of detected photo-electrons. The scale factor between energy deposited and expected photo-electrons 
detected is determined by a standalone optical simulation validated for the MINOS experiment; the 
expected number of photo-electrons is smeared by Poisson statistics, and a 10% channel-to-channel 
Gaussian smearing reflecting a conservative estimate of remaining systematics after calibration and 
attenuation corrections. 

In addition to the GEANT Monte Carlo, a photon transport Monte Carlo written by Keith Ruddick| 198 1 
for the MINOS experiment was used to optimize the strip and fiber dimensions. The average light yield 
from a MINOS scintillator module is 4.25 photo-electrons/MIP at a distance of 4 meters, and attenuation 
in the fiber is well described in terms of a double exponential 1 199|: 

N{x) = A(e--'/90 cm + g-./roo cm^ (^3) 

The photon transport Monte Carlo (LITEYLDX) is used to calculate the number of photons trapped in 
the fiber for a MIP entering at a particular position and for a given configuration of strip geometry, fiber 
diameter, and fiber placement. This information is then used to determine a relative light collection 
efficiency for a particular configuration compared to MINOS strips. With the overall normalization and 
attenuation curve from MINOS one can then calculate the amount of light for any particular configu- 
ration. Figure ^2 for instance, shows the relative light output for triangular extrusions when the strip 
thickness, fiber diameter and fiber placement are varied. As expected, light output is nearly proportional 
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Figure 70: Relative light collection efficiency across the 3.35 cm triangular width of the scintillator 
extrusion. 

to the strip thickness, and is greatest when the fiber is placed at the center of gravity of the strip. Having 
a fiber in a groove at the edge only results in a 9% drop in the light level. Figure ^Ol shows the relative 
light collection efficiency for a triangular extrusion where the entry point of the minimum ionizing par- 
ticle is varied across the strip width, and indicates that the collection efficiency varies by ±10% over 
the strip width. 



15.5 Event Reconstruction 

The output of the detector simulation comprises a list of digitizations for each strip. We have developed 
a basic reconstruction program to take this Ust and reconstruct the tracks and vertices in an event. 



15.5.1 Pattern recognition 

Development of a fully-realistic pattern-recognition algorithm to associate hits to track candidates was 
not undertaken, in view of the manpower and time available. We are confident that the three-dimensional 
XUXV modular design of the detector, and its relatively modest occupancy, will allow highly-efficient 
pattern recognition and track identification. Visual inspection of events through the graphical interface 
of the detector simulation program reinforces this conclusion. For our design studies, we have adopted 
"omniscient" pattern recognition based on Monte Carlo truth information. All hits generated by a given 
track (ignoring channels with overlap) are used to reconstruct it. 
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Light Yield for Strip and Fiber Variations 
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Figure 71: Relative light yield for different strip widths and fiber diameters. 

15.5.2 Coordinate reconstruction 

Tracks generating hits in at least six scintillator planes of the inner detector, including three planes of 
the X view, are reconstructed. Coordinates are estimated from the raw, smeared digitizations, using only 
planes which have one or two strips hit. Tracks at high angles to the detector axis may pass through 
more than two strips in a single plane, and it should be possible to recover these higher-multiplicity hits 
with a more sophisticated algorithm. For single hits, the coordinate is taken as the center of the strip. 
For dual hits, the position is interpolated using the charge-sharing between between strips, with a small 
geometrical correction based on the estimated crossing angle. 

The coordinate resolution for a large test sample of single and double hits can be measured directly 
using the residuals obtained when each coordinate is excluded, in turn, from the track's fit. This coor- 
dinate resolution is parameterized as a function of the track's crossing angle, and used to assign errors 
to coordinates in the fitter. 

15.5.3 Track reconstruction 

Reconstructed coordinates are used to fit each track using a Kalman filter algorithm|200|. For this pro- 
posal, tracking performance has only been studied in the non-magnetic region of the detector; the track 
model is perforce a strictly linear one. Neglect of the magnetic field is justified because mission-critical 
resolutions are determined by performance of the fully-active (non-magnetized) volume, and since co- 
ordinate resolution for the strips should not depend on the presence of a magnetic field. The momentum 
resolution for charged tracks in a magnetic field can be reliably estimated from the coordinate reso- 
lution, momentum and field strength li48J . As long tracks may pass through many radiation lengths of 
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scintillator and absorbing material, the Kaknan filter's ability to correctly account for multiple Coulomb 
scattering ("process noise") is essential. The algorithm can optionally be used to exclude outliers from 
the fit. 

Figure EH shows the expected hit residuals, impact parameter and angular resolution for muons 
from a sample of quasi-elastic interactions, assuming triangular strips of 3 cm width and 1.5 cm thick- 
ness (close to the final design values). Hit resolutions of ~ 3 mm and angular resolutions of < 0.5° 
are expected. The coordinate resolution is degraded to approximately 1.5 cm if rectangular strips are 
employed instead of triangular ones, since interpolation based on charge is no longer possible. 

Tracking Performance 





Mean -0.1446E-01 
RMS 0.3406 


\ \ 1 \ ''''II -r-'-rn \ \ \ \ \ \ \ 1 





15000 



10000 



-1.5 



-0.5 



0.5 



1 1.5 2 

Hit residual (cm) 



400 
300 
200 
100 





Mean 
RMS 



0.4100 
0.5080 



0.25 



0.5 



0.75 



1.25 1.5 1.75 2 2.25 2.5 
[I Impact parameter (cm) 




3 3.5 4 

|j, angular resolution (deg) 



Figure 72: Performance of the tracking algorithm on muons from from a sample of simulated charged- 
current quasi-elastic interactions. Shown are (top) the hit residuals, (middle) the impact parameter of 
the muon with the vertex and (bottom) the muon angular resolution. 



15.5.4 Vertex reconstruction 

In this study, reconstructed tracks are associated to vertices using Monte Carlo truth information. The 
vertex positions are then fit using a Kalman filter algorithm. Track directions at the vertex are updated 
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taking account of the constraint. This is equivalent to a least squares fit, but mathematically more 
tractable since it does not involve inversion of large matrices and can be easily extended to a helical 
track model. The primary vertex resolution for a sample of simulated quasi-elastic interactions with 
two visible tracks is shown in Figure ED The vertex postion can be measured to a precision of better 
than a centimeter. 
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Figure 73: Reconstructed vertex resolution for two track charged current quasielastic events. Shown 
are (top) the resolution in the longitudinal position of the vertex (Z) and (bottom) the resolution of the 
transverse position of the vertex (X and Y). 



15.5.5 Particle identification 

Particle identification in MINERi/A will rely on measuring specific energy loss (dE/dx) as well as 
topology (hadron and electromagnetic showers, decay signatures). 

Electromagnetic showers Electromagnetic showers are easily identifiable by their diffuse track and 
characteristic dE/dx profile in the fully-active central detector and energy deposition in the electromag- 
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netic calorimeters. Section fT3.6. ll describes a preliminary technique to separate electrons and photons, 
when the primary vertex is known, using distance to shower onset and shower length (Figure I5HI. 

Muons Energetic muons can be identified by their penetration of material in the calorimeters and/or 
MINOS near detector. Muons with a momentum measurement in the magnetic field, or which stop 
inside the detector can be distinguished from protons and kaons by dE/dx. In addition, the delayed 
fi ^ e decay signature can be detected. 

Hadrons Hadrons can be identified as such by their interactions in the inner detector and/or hadron 
calorimeters. Hadrons which stop without interacting or have their momentum measured by the mag- 
netic field can also be distinguished as it, K or p with good efficiency using dE/dx. 

dE/dx analysis Specific energy loss (dE/dx) will be an important tool for particle identification 
in MlNERz^A. For tracks which stop in the inner detector, the charge deposited near the end of the 
track (corrected for sample length) can be compared with expected curves for, e.g., the vr^, and 
proton hypotheses. This technique does not require an independent momentum measurement, since the 
range (Xstop, in g/cm^) from the stopping point to a given sampling point is closely correlated with the 
momentum at the sampling point. The algorithm is calibrated by fitting the expected dE/dx vs. Xstop, 
and the standard deviation of this quantity, (T^E/dx^ a function of Xstop for the three different particle 
types (see Figure l74t. The measured dE/dx for a track is compared to the expected value at each 
sample, to form estimators reflecting the goodness of fit to each of the three particle identification 
hypotheses: 

^sample 

i=l 

where the sum runs over all measured samples, and a = {vr, K, p}. The hypothesis a with the minimum 
is assigned to the track. The frequency of misidentification can be visualized most easily by plotting 
the difference Ax^ between the correct (for the particle's true type) and the smallest of the two 
(incorrect) others (Figure l75l. With this naive dE/dx analysis, MlNERi^A correctly identifies 85% of 
stopping kaons, 90% of stopping pions, and > 95% of stopping protons. A similar analysis can be 
applied to tracks with momenta measured in the magnetic regions of the detector. 

15.5.6 Energy reconstruction and containment 

The energy of muons from charged-current interactions will be measured using range and/or curvature 
in the magnetized regions of the detector and the MINOS spectrometer. For muons stopping in the 
detector, the momentum resolution will be ^ ~ 5%. If the MINOS detector is used, the momentum 
resolution will be 13%|1911. Preliminary work on hadronic energy reconstruction suggests that the 
energy of hadrons which rangeout in the detector will also be measured to a precision of 5% whereas 
the resolution for isolated showering hadrons will be 35%/-\/i?(GeV). The resolution for hadronic 
showers in deep-inelastic scattering will be approximately 55%/\/^ assuming the "ideal" light collec- 
tion and smearing effects described in section fT54l For electromagnetic showers, the estimated energy 
resolution is 6% / ^/E. 
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Figure 74: The top figure sliows tiie average specific energy loss dE/dx for stopping vr^, kaons and 
protons, vs. range from the stopping point (in g/cm^), for the simulated MINERi^A inner detector. The 
bottom figure shows the estimated standard deviation of the energy loss, which is used to form a 
estimator for particle identification. 
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Figure 75: The three plots show the Ax^ dE/dx estimator for simulated and reconstructed charged 
pions(top), kaons(middle) and protons(bottom) stopping in the inner detector. Tracks with Ax^ < 
are correctly identified. 
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Containment of hadronic energy is a significant design consideration, as it assists in meeting many 
of the experiment's physics goals. Studies show that the visible hadronic component of quasi-elastic 
and resonant events in the fully-active central region of the detector are completely contained, apart 
from secondary neutrinos and low-energy neutrons. Figure shows the fraction of escaping visible 
hadronic energy for deep-inelastic reactions in several hadronic energy ranges, and figure [77\ shows 
the probability that a deep-inelastic event will leak visible energy as a function of the true hadronic 
energy. Only for hadronic energies greater than 8 GeV is there any significant probability of leakage 
and only above 15 GeV is the average fraction of escaping energy greater than 10%. The fraction 
of deep-inelastic interactions with hadronic energies over 15 GeV in the low-energy, semi-medium or 
semi-high energy beams is < 1%, and so visible energy leakage should be insignificant. These estimates 
ignore downstream components beyond the forward hadron calorimeter, such as a muon ranger and/or 
the MINOS detector, and are therefore conservative. 



Energy Leakage for DIS events 



10 



10 



10 



10 



EHad<10GeV 

10 GeV < Ejj^^ < 15 GeV 

15 GeV < Ejj^^ < 20 GeV 

20GeV<Ejj^, 



0.05 



0.1 



0.15 0.2 0.25 0.3 

Fraction of hadronic energy escaping detector 



Figure 76: Fraction of hadronic energy escaping the detector for deep-inelastic scattering in the fully- 
active central region. 



15.6 Event Categorisation 

Particle identification and event classification will play a central role in the analysis of data from 
MINERz/A. One possible method of event classification is use of artificial neural network (ANN) 
techniques. 
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Probability of hadronic energy leal^age for DIS events 




Figure 77: Probability that visible hadronic energy from a deep-inelastic event escapes undetected vs. 
total hadronic energy. 



Event classification will be based on on topological characteristics as well as on particle ID. Separa- 
tion of CC from NC interactions will be based on muon identification. Detection of muon decays for low 
energy muons stopping in the carbon gives the potential for accurate CC identification even at high y Bj ■ 
In each such class further event identification will be based on other particle ID, energy/momentum 
measurements and kinematics. Neural networks are designed for such categorisation and have been 
frequently used in the analysis of data from high energy physics experiments (see, for example, the 
DONUTlSnil experiment). 
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16 Detector Design 



This section describes the basic elements of the MINERi/A detector, including the arrangement of ac- 
tive elements and absorber, photosensors and scintillator strip details, and the electronics. A summary 
of detector parameters along with an estimate of costs and construction schedule are provided in Sec- 
tionfTTl 

Muon Ranger (MR) Outer Detector (OD) Veto 








Figure 78: A schematic side view of the MINERi/A detector with sub-detectors labeled. The neutrino 
beam enters from the right. 

The MINERz^A detector is made up of a number of sub-detectors with distinct functions in re- 
constructing neutrino interactions. The fiducial volume for most analyses is the inner "Active Target" 
shown in Figure ^81 where all the material of the detector is the scintillator strips themselves. In other 
regions of the detector, the strips are intermixed with absorbers. For example, the side, upstream (US) 
and downstream (DS) electromagnetic calorimeters (ECALs) have lead foil absorbers. Surrounding the 
ECALs are the US and DS hadronic calorimeter (HCAL) where the absorbers are steel plates. On the 
side of the detector, it is the outer detector (OD) that plays the role of the HCAL; however, note also 
that the OD is a magnetized toroid which will focus and bend muons, thus allowing a momentum mea- 
surement for muons which exit the detector. Upstream of the detector is a veto of steel and scintillator 
strips to shield MINERi^A from incoming soft particles produced upstream in the hall. Finally, the most 
downstream element, the muon range detector/toroid (MR) gives MINERz^A the capability to fully re- 
construct even high energy muons without the use of the MINOS near detector as an external muon 
spectrometer. The presence or absence of the MR in the final design will depend upon the location 
chosen for MINERi/A. 
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16.1 Overview of MINERz/A Detector Design 

For MINERi/A to meet its physics goals, detector must break new ground in the design of high-rate 
neutrino experiments. With final states as varied as high-multiplicity deep-inelastic reactions, coherent 
single-vr'^ production and quasi-elastic neutrino scattering, the detector is a hybrid of a fully-active fine- 
grained detector and a traditional calorimeter. 

At the core of the MINERz^A design is a solid scintillator-strip detector, similar in principle to 
the recently commissioned K2K SciBarl202J. The plastic inner detector serves as the primary fiducial 
volume, where the precise tracking, low density of material and fine sampling ensures that some of the 
most difficult measurements can be performed. These include multiplicity counting in deep-inelastic 
scattering, tracking of photons, detection of recoil protons in low-Q^ quasi-elastic events, and particle 
identification by dE/dx. 

The scintillator detector cannot contain events due to its low density and low Z, and therefore, 
the MINERi/A design surrounds the scintillator fiducial volume with sampling detectors. At the low 
energies needed to study cross-sections of interest to neutrino-oscillation studies, many of the events 
contain sideways-going and backward-going particles, and therefore these sampling detectors extend 
to the sides, and even to the back of the detector where they also serve as high A targets for studies 
of nuclear dependence in cross-sections. Finally, it is important to contain or measure the final-state 
muon in charged-cunent events, and for this purpose, the outer side detector and downstream muon 
ranger of MINERi/A are magnetized toroids. A side view of the complete MINERz^A design is shown 
in Figure l79l 

The sensitive elements of MINERz/A are extruded triangular scintillator strips, 1 .7 cm height with a 
3.3 cm base, embedded with WLS fibers as detailed in Section (16.41 To improve coordinate resolution 
while maintaining reasonably large strips, these elements are triangular in shape and assembled into 
planes as shown in Figure [80l this allows charge-sharing between neighboring strips in a single plane 
to interpolate the coordinate position. Calorimetric detectors and nuclear targets in the central region of 
the detector are constructed by inserting absorber between adjacent planes as illustrated in Fig ureiU In 
the outer detector (OD), strips of steel absorber and scintillator are assembled in a picture frame around 
the inner detector. In the case of the triangle, the scintillator strips are not the full size, but rather half 
(right) triangles, 1.7 cm in height with a 1.65 cm base and are assembled in doublets between steel 
absorber strips. 

For construction and handling convenience, a single plane of MINERi/A, shown in cross-section 
in Figure |82j incorporates both the inner detector and OD "picture frame" as well as an outer picture 
frame support structure. Groups of four planes (occasionally two planes only in the upstream veto 
and downstream muon ranger components) are ganged together into modules, again as illustrated in 
cross-section in Figure |82l There are three distinct orientations of strips in the inner detector, muon 
ranger and veto, separated by 60°, and labelled X, U, V. A single module of MINERi/A has two X 
layers to seed two-dimensional track reconstruction, and one each of the U and V layers to reconstruct 
three-dimensional tracks. The 60° offset makes the hexagon a natural transverse cross-section for the 
detector, and the size and shape of MINER;/ A are illustrated in Figure 183) 

Except for the upstream veto and downstream muon range (MR) detector, the entire MINERz/A 
detector is segmented transversely into an inner detector with planes of solid strips and an outer picture 
frame magnetized toroid (OD). In Figure^O] the upsteam and downstream most detectors, the veto and 
muon range toroid, respectively, are shown in Figure |84] As shown, the scintillator strips extend the full 
length of the hexagon and range between 205 and 400 cm in length. The toroid steel/absorber is 15 cm 
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Basic element: 1.7x3.3cm triangular strips. 
1.2mm WLS fiber readout in groove at bottom 



Assemble 
into planes 




Figure 80: Assembly of scintillator strips into planes. 




Figure 8 1 : Integration of planes with absorbers in calorimeters or nuclear targets in the inner (above) 
and outer (below) detectors. 
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Figure 82: Plane assembly (left) and module assembly (right) in the active target region for MINERi/A. 
On each drawing, the scale is exaggerated in the horizontal direction to show details. 
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Figure 83: Outline of MINERi/A detector to illustrate shape and scale. Note the locations of the PMT 
readout boxes on top of the detector, coils on the bottom, and the support stands. 

thick in the muon ranger and 10 cm thick in the veto; note also that the final module in the muon ranger 
is constructed without steel to ensure one final three dimensional spatial point free of local multiple 
scattering. The magnetic properties of the OD and the MR detectors are discussed in Section fT6.2l 

Moving towards the center of the detector from each end, the next detectors are the downstream 
and upstream hadronic calorimeters (HCALs), shown in Figure l85l with 2.5 cm absorbers, one per 
plane downstream and one per module upstream. This detector is surrounded by the picture frames of 
absorber and scintillator strips that make up the outer detector (OD). Note that the strips in the OD run 
only in one direction, in the bend plane of the magnetic field. Three-dimensional tracks must therefore 
be matched from the inner detector and extrapolated outwards for an energy measurement or muon 
momentum measurement. A complication of the design is illustrated by the fact that the inner detector 
strips, which range in length from 120 to 240 cm, end inside the OD, and therefore the WLS fibers 
from must be routed out to the detector edge through a grooved plastic guide plate through the region of 
the OD. Note also the holes for the OD muon toroid coil in the lower region of the detector. Magnetic 
flux will be isolated in each region frame of the OD, and will be prevented from leaking into the inner 
detector by a guard ring of stainless steel as part of the HCAL absorber. 

Moving in again from upstream and downstream, the next detector module elements are the electro- 
magnetic calorimeters (ECALs), which have 0.2 cm Pb/Stainless absorbers downstream, one per plane, 
and 0.8 cm Pb absorbers upstream, one per module. Their design is shown in Figure |86l Note that the 
absorber only overlaps the inner detector and not the outer detector where it would represent a negli- 
gible fraction of the absorber material. The fine granularity of the ECAL ensures excellent photon and 
electron energy resolution as well as a direction measurement for each. 

Finally, we reach the center of the detector, the fully-active inner detector (ID), whose plastic core 
represents the fiducial volume for most analyses in MINERz^A. A plane of the active target is shown 
in Figure |^ In the center region, there is no absorber at all; however, 30 cm from the edge of the 
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Figure 84: Muon range/toroid and upstream veto plane design (landscape). 
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Figure 85: The Hadronic Calorimeter Plane Design (landscape). 
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Figure 86: Electromagnetic calorimeter plane design (landscape). 
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ID, there are lead/stainless absorbers identical in thickness to the downstream ECAL, which act as a 
side electromagnetic calorimeter. This part represents the bulk of the detector in length, and the outer 
calorimeter surrounding the fully active planes are the largest part of the detector in mass. 

Note that MINERi^A is, by design, entirely modular along the beam direction. Individual elements 
may be easily lengthened or shortened by omitting modules from the design or adding new modules. 
One configuration that would be attractive is to forgo installation of the muon ranger and perhaps a 
portion of the downstream HCAL in order to move as close as possible to the front face of the MINOS 
near detector, thus allowing MINOS to serve as a calorimeter and muon detector. For the purposes of 
cost and schedule, however, we proceed to make estimates under the assumption that the full stand-alone 
detector will be built. 

16.2 Muon Toroid Performance 

The MINERi/A design calls for toroidal muon spectrometers in the outer detector (OD) and downstream 
of the HCALs in the muon range (MR) detector. (Again, however, it should be noted that if MINERz^A 
were situated immediately upstream of MINOS, the downstream muon toroid may be omitted.) This 
section describes the momentum reconstruction and range capabilities of these detectors for fj, produced 
in the inner plastic fiducial volume. 

The OD has a total of 50 cm of magnetized steel sampled by active planes that are traversed by 
muons in a direction perpendicular to the beam. It is magnetized by a 48 turn coil with 700 Amp 
current. The average magnitude of H in the OD is therefore about 30 Gauss. We plan to use Armco 
specialty steel 12031 for the OD absorber which would give a magnetic field of about 16 kGauss. For 
muons which exit the side of the OD, the fractional momentum resolution measured from the bend 
angle varies from 22% to 30% for muons with an angles of 30° to 90° with respect to the beam. In 
practice, of course, the resolution will be better because of the loss of momentum with dE/dx in the 
OD. The OD will run to focus muons forward with a transverse momentum kick of 0.5 GeV (0.25 GeV) 
30° (90°) angle. Focusing will serve to lengthen the path length through the OD and to direct the muons 
into a downstream muon range detector, be it the MINERz/A MR or the MINOS near detector. 

The downstream MR toroid has a total thickness of 1.2 m of magnetized steel with a 48 turn coil 
and 1200 Amp current, resulting again in an average field of 16 KGauss. This yields a typical pT kick 
of 0.6 GeV and a momentum resolution of 20% from the bend, which is, again, improved by the muon's 
energy loss in passing through the steel. 

In summary, the MINERi/A detector has, on its own, excellent acceptance and momentum resolu- 
tion for muons. This resolution can be improved, especially for forward-going high-energy muons, by 
use of the MINOS near detector as a downstream muon toroid. 

16.3 Photosensors for MINERz/A 

With an inexpensive active detector technology, the dominant equipment costs for MINERz^A are pho- 
tosensors and their associated readout electronics. The path through the parameter space of available 
technologies is determined by the answers to three questions. First, is the light output of the detector 
for a MIP signal sufficient to support a low quantum-efficiency detector such as photomultiplier tubes 
(PMTs) or image intensifier tubes (IITs)? For MINERi/A there is sufficient light to use a 1/6 quantum 
efficiency photocathode with a WLS fiber diameter of at least 1.2 mm as demonstrated in 116.41 Second, 
is timing within the spill important or can a technology that only integrates over a long spill, such as 



155 



nTs be used? We concluded that timing within the spill, both to flag overlapping events and measure 
time of flight and decay times at rest was important for our physics goals. Third, what level of technical 
risk, R&D time and cost is acceptable? We concluded that to allow MlNERz^A to operate as early as 
possible in the NUMI beamline and given the modest size of our collaboration and expected detector 
costs, we should choose low technical risk over lengthy R&D programs designed to reduce those costs 
or improve performance. 

In our design exercise, we considered four technologies for photosensors: multi-anode photomulti- 
plier tubes (MAPMTs), IITs, avalanche photodiodes (APDs) and visible light photon counters (VLPCs). 
Ultimately, we chose to pursue a solution based on MAPMTs which results in a sensor+electronics 
cost (including EDIA and overhead but without contingency) of approximately $40 per channel, which 
breaks down approximately as $15 per channel for the sensor, $15 for the electronics and $10 for EDIA 
and testing. To defend this important decision, we discuss the alternative technologies mentioned above. 

Image intensifying tubes coupled to CCDs as a readout device are an extremely appealing low cost 
solution for reading out bundles of fibers, in part because the CCD itself is the final stage photosensor 
and readout device. This device is well-matched to the pulsed structure of the neutrino beam with 
one readout corresponding to one beam pulse. Costs per channel are largely proportional to the total 
photocathode surface required, which is set by the number of channels and fiber diameter. Cross- 
talk in adjacent channels is a non-trivial issue, but can be addressed because of the high density of 
CCD channels relative to fiber granularity, even with intermediate spatially demagnifying stages. We 
were driven to relatively expensive CCD cameras because of the need maintain reasonable linearity. 
Our candidate system, based on Hamamatsu C8600 2-stage multi-channel plate (MCP) intensifiers 
and C7190 bombardment CCDs, was approximately $15 per channel, including photosensor and CCD 
readout but not including required demagnification optics. Nevertheless, a complete IIT/CCD system 
would still likely be half the cost of the chosen MAPMT solution. Our concerns about the system were 
the smaller effective dynamic range, even with relatively costly IIT/CCD systems, and the relatively low 
mean time to failure per device reported in other large systems (4 years per two stages in the CHORUS 
experiment). However, the missing capability of timing within a single main injector spill was enough 
for us to discard this otherwise promising option. 

Avalanche photodiodes (APDs) were also considered because of their recent successful application 
in the CMS ECAL and their proposed use in the NUMI off-axis far detector. APDs are low gain (~ 100), 
high quantum efficiency (85% for Yll WLS fibers) devices which offer significant cost savings in the 
photodetector. Complications of operation include the need to cool the sensors below room temperature 
to reduce noise, but this is a fairly easily solved problem as cryogenic temperatures are not required. 
The primary problem we identified with APDs for MINER^'A was the need for significant electronics 
R&D to develop a relatively low-cost system capable of controlling noise over the long NUMI spill. 
For MINERt/A we set a requirement of keeping the photosensor and electronics noise well below 10 
delivered photon equivalents to maintain good sensitivity to a MIP (typically 70 photons in a doublet 
of triangular scintillators) and a low rate of detector noise. Over a 12 |Ltsec gate (the NUMI spill plus 
2t^) at -10°C with an operating gain of 100 (optimal), the signal from 10 photons is 850 electrons and 
the noise on the best existing candidate electronics, the MASDA chip, is 900 electrons. To achieve the 
better signal to noise that is the goal of the proposed NuMI off-axis R&D program requires design of a 
new ASIC, which would imply at least a one-year development project. In short, although the APD is 
a potentially promising technology, we were not convinced it could be in production on the timescale 
required for MINER;yA. 

The final option we considered was the VLPC. These have the advantage of successful past deploy- 
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ment and electronics design in the DO fiber tracker and preshower detectors. However, the costs for 
just the VLPCs themselves, even with optimistic assumptions about the outcome of future R&D, are 
expected to exceed $50 per channel, and are thus significantly more expensive than the MAPMT solu- 
tion. Given that the low quantum-efficiency solution gives sufficient resolution, it is difficult therefore 
to justify VLPCs. 

The MAPMT we have tentatively selected as our default photosensor is the Hamamatsu R7600U- 
00-M64. These are an incremental design improvement from the R5900-00-M64 MAPMTs used in the 
MINOS near detector, and we expect much of the experience gained by the MINOS collaboration with 
these detectors to be applicable. In particular, we have confidence in costing the testing, housing for 
and optical connectors to the PMTs because of our ability to scale costs from the MINOS experience. 

Having chosen MAPMTs, a low quantum-efficiency device with good timing, low noise and a large 
dynamic range, the following two sections address the issues of the photoelectron yield from the strips 
married to the MAPMTs and the electronics to readout these MAPMTs, respectively. 

16.4 Scintillator Strips 

The MINOS experiment has successfully demonstrated that co-extruded solid scintillator with embed- 
ded wavelength shifting fibers and PMT readout produces adequate light for MIP tracking and can be 
manufactured with excellent quality control and uniformity in an industrial setting. The performance 
characteristics of the MINOS scintillator modules produced at the three 'module factories' are now well 
known, both through measurements taken with radioactive sources post-fabrication at the factories and 
through measurements of cosmic rays at Soudan. We intend to use the same technology for the active 
elements of MINERz^A. 

The basic active element in MINERi/A is a co-extruded triangular scintillator strip with a wavelength- 
shifting fiber glued into a groove. Like MINOS, the scintillator strips are polystyrene (Dow 663) doped 
with PPO (1% by weight) and POPOP (0.03% by weight), co-extruded with a reflective coating of Ti02 
loaded polystyrene|204j. The strip cross-sections have width 3.35 cm and height 1.7 cm. Strip lengths 
vary throughout the detector and range from 1.4 meters to 2.2 meters in the inner tracking detector 
to 4 meters for the veto and muon ranger sections. The WLS fiber (Kurrary Yll, 175ppm dopant) is 
1.2 mm in diameter, glued into an extruded groove and covered with aluminized mylar tape in the same 
fashion as MINOS. The WLS fibers are brought to optical connectors at the edge of the modules, and 
clear optical cables bring the light to a PMT box. Single-ended readout is used, and the far strip/fiber 
ends are mirrored. 

Physics simulation studies indicate that for a triangular extrusion, average ^ light levels above 3.9 
photo-electrons("PE")/MeV of dE/dx for a minimum-ionizing particle (MIP) are required in the inner 
detector in order to obtain good particle identification as described in Section fT5.5.5l Coordinate reso- 
lution, vertex finding, and track pointing are also affected by light levels, but to a lesser extent. For this 
design we have targeted an average light level of 7.8 PE/MIP on average through the strips. This allows 
for losses expected to be 25% in the clear fiber and connectors and possible effects from the degradation 
of the scintillator over time, the latter of which was measured to be as large as 20% over 10 years for 
MINOS 1 205 ]. 

The overall light levels from 3 lengths of strips, as calculated using the photon transport Monte 

'Note that this is an not only an average over photostatistical fluctuations, but also an average over all locations for normally 
incident tracks to enter the strip. The average light through the full thickness of scintillator in a plane, a doublet of triangles, 
is twice this average. 
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MINERVA Light Yield With Mirrored Strip Ends 
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Figure 88: Light yield vs. distance along strip for MlNERz^A scintillator strips with one-ended readout 
and a mirrored end. Dot-dashed line is light collected from reflections off the mirrored end, dashed 
line is light travelling directly to the readout end; solid line is the sum. The top plots correspond to the 
shortest and longest strips used in the fully-active inner detector, and the bottom plot is for the longest 
strips in the veto and muon ranger. 



Carlo described in Section fT5.4l are shown in Figure[88l Here we have assumed a 90% reflectivity from 
the mirror end of the fiber, and in all cases a 1 meter WLS 'pigtail' from the end of the near end of the 
strip to the PMT face. Clear fiber lengths and connectors are not included. In the MINOS near detector, 
the far strip end was not mirrored; here we assume the strip ends are mirrored with 100% reflectivity. 
Because the light produced in the scintillator is generally collected within a few cm of the MIP crossing 
location, this approximation only affects the calculation of collection efficiency at the very far end of 
the strip. Shown are the light levels predicted for three strip lengths. In each plot, the lowest curve 
corresponds to light collected from reflections off the mirrored end, the middle line corresponds to light 
travelling directly from the MIP to the readout end, and the upper line is the sum. As the figure shows, 
the light level in the inner tracking detector, with a maximum length of 2.2 m, exceeds the design 
requirement of 7.8 PE/MIP over the entire length by about 25%. In the longer strips (only used in the 
downstream muon range detector and upstream veto counters), the light falls slightly short of this target 
at the far ends of the strips; however, because these detectors are not used for particle identification by 
dE/dx, this is still acceptable. 
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16.5 Electronics 



The requirements for the MINERz^A electronics are summarized in Table \T4\ To minimize technical 
risks, we studied a number of existing solutions, including those used for the MINOS design. Major 
components of the electronics system include the front-end boards, the PMT and electronics housing 
and the slow controls and readout systems. 

16.5.1 Front-end boards 

For the front-end digitization and timing, the best performing system with the least required R&D is 
a scheme based on DO TRiP ASIC. The TRiP chip is a redesign of the readout ASIC for the DO fiber 
tracker and preshower originally motivated by the need to run at 132 ns bunch crossings in the TeVatron. 
A production run of one version of TRiP with 7000 produced chips has been completed; however, since 
the TeVatron run plans do not now call for 132 ns operation, these chips will not see their original use. 
These existing chips, however, could be recycled into use in MINERzvA. 

TRiP was designed by Abder Mekkaoui of the Fermilab ASIC group and successfully met the specs 
in its first submission and has undergone extensive testing by DO II206I . Its analog readout is based 
on the SVX4 chip design, and each TRiP chip supports 32 channels for digitization, but only half that 
number of channels for timing. A simplified schematic of the TRiP ASIC is shown in Figure [89l The 
preamp gain is controlled by SW2 and has two settings which differ by a factor of four. The gain of the 
second amplifier stage is controlled by SW3-SW5. We will set the chip to the lowest gain setting for 
the preamp and largest integration capacitor. This gives a linear range with a maximum charge readout 
of 5 pC. The 'ANALOG_OUT" goes into a analog pipeline, which is identical to the one used on the 
SVX4 chip and 48 cells deep. The SVX4 chip can read out four of these 48 buffers, and although the 
TRiP chip was also designed to read out four buffers, it can empirically only read out one buffer. It is 
not known why only one buffer can be read out; however, this is not an issue in MINERi/A as shown 
by the per channel per spill occupancy illustrated in Figure|90l To gain dynamic range, MINERz^A will 
increase the input range of the electronics by using a passive divider to divide charge among two TRiP 
channels with a ratio of a factor of 10. This "high range" channel, then, will give a equivalent total 
readout charge of 50 pC. Each TRiP channel will be digitized by a 12 bit ADC. 

In MINERz^A the integration time for the ADC will be 10-12 fis, much less than the hold time for 



Parameter 


Value 


Comments 


Active Spill Width 


12;usec 


Spill plus 2t^ 


Repetition Time 


> 1.9 sec 




Number of Channels 


37478 




Occupancy per Spill 


0.02 


LE beam, 2.5E13 POT/spill 


Front-end noise RMS 


< 1 PE 




Photodetector gain variation 


4.5 dB 


extremes of pixel-to-pixel variation 


Minimum Saturation 


500 PE 


proton range-out or DIS event 


Maximum Guaranteed Charge/PE 


50 fC 


lowest possible charge at highest gain 


Time Resolution 


3ns 


Identify backwards tracks by TOE 






Identify decay-at-rest 



Table 14: Electronics design requirements and parameters for MINERi/A 
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Q IN 




lANALOG OUT 



DISCR OUT 



SW6 



Figure 89: A simplified sciiematic of tlie front end electronics of the TRIP chip 




ID+OD+Cal Strip Occupancy Per SpUl, 2.5E20 POT, LE+ beam 



Figure 90: The mean number of neutrino interactions in which a strip is hit per 2.5E13 POT spill, NUMI 
LE neutrino beam 
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Histogram of measured signal w/ 5fC injected 
10us gate 
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Figure 91: Response of the TRiP chip to 5 fC injected with a 10 /is gate. 



the charge in the capacitor of 100 /xs. The TRiP chip has been tested explicitly with a 10 fis gate, and 
Figure |^ shows pedestal RMS of 3 fC for the estimated MINERi/A input capacitance of 36 pF The 
MINERz/A design requires no saturation below 500 photoelectron (PE) and RMS noise well below 1 
PE. Matching this to the 5 pC charge limit, the highest gain anodes in a tube would be set at 100 fC/PE 
and therefore the lowest gain anodes would be run at 33 fC/PE. Since the RMS of the noise is about 
3.0 fC, this will put a single photoelectron approximately a factor of 10 above the pedestal RMS, well 
within our design spec. The maximum PMT gain for the lowest gain anode will be 50 fC/PE, safely 
within the desired parameters above. 

Only one of every two input channels to the TRiP chip have a latched discriminator output (latch) 
which can be used for timing information. Hence, only the lower range channels will feed the latch 
whose output will then go into an FPGA. The internal clock of the FPGA can be used to get timing 
with a granularity of 5 ns, and with a delay line scheme this can be improved to below 3 ns. The reset 
for the latch is only 15 ns, so inside the spill the latch will be in the ready state by default. When 
the signal exceeds a threshold of 1.5 PE, the latch will fire. After storing the time, the latch is reset, 
incurring minimal deadtime. Figure |92l shows result of the DO timing test of the TRiP chip using their 
fiber tracker and VLPCs. They get a timing RMS of 3ns for signals with >8 PE. Y-11 (the waveshifter 
in the MINERi^A fibers) has an equivalent decay time to 3HF (the dye used in DO's fiber tracker), and 
hence in a doublet of triangles in the scintillator, we can reasonable expect similar timing resolution. 

Note in the MINERz^A scheme, there is no trigger. Two charges are read from all channels along 
with all latched times at the end of each spill. 

Although individual parts of this system have been tested by DO, the system described above has 
not been tested. Ray Yarema's group at FNAL has begun layout of a board for a vertical slice test of 
this system, and we expect a proof of principle from this test by early summer 2004. 
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Figure 92: Results of the DO test of timing of TRiP using the Test Fiber Tracker with signals > 4pe and 
signals > 8pe 



16.5.2 PMT boxes, readout electronics and controls 

The MAPMTs will be mounted directly on the front-end boards to reduce input capacitance. Therefore, 
the front-end board and PMT need to reside in a single light-tight box with optical cables from the 
detector as input. In our preliminary design, each PMT box will have a single PMT and front-end 
electronics for its 64 channels, along with a Cockcroft-Walton HV generator. 

In addition to the optical input cable, each PMT box has three electrical cables. The slow-control 
cable and the low-voltage cable will travel from box-to-box in a daisy chain. The digitial readout LDVS 
cables will arranged in 16-box 'Token Rings', and will connect to a VME card at the first and last box 
in the ring. In addition to reading out the data after each spill, the token ring will supply the timing 
synchronization signal. The low-voltage cable will likely run at an intermediate DC voltage and step 
down at each box to minimize the role of resisitive losses in the chain. The slow-control cable will 
be a MIL 1553b bus which is in wide use at FNAL. An existing VME card (the 1553 controller) will 
drive the slow control system. Table fTSl shows the number of parts needed for the complete electronics 
system. 

To minimize the length of the clear fiber cables from the WLS fibers to the PMT boxes, we plan to 
mount the PMT boxes directly on the upper parts of the MINERi^A detector, on the two highest sides of 
the hexagon to avoid conflict with the coils or side clearance of the detector. This will require magnetic 
shielding of the MAPMTs. 

The DAQ requirements for this system are trivial as the data rate is expected to be under 100 kByte/second. 
A VME-resident PVIC interface in each of the four VME crates will be readout with a single Linux PC. 
Data will be buffered in a local RAID system and transferred over the network to FCC for permanant 
storage. 
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16.5.3 Whither the TRiP chip? 



DO is likely to redesign the TRiP chip as part of the fiber tracker front end upgrade. As part of that 
upgrade, features desirable for MINERvA, including individual channel discriminator threshold and 
front end buffer gain could be added. A major goal for this submission, independent of MINERz^A, 
would be to get the multi-buffer readout mode, which would be a useful safety valve if rates in portions 
of the detector were higher than we now predict. The MINERz^A readout board would be able to use 
either chip without any modification to the board; hence, development of the front-end board and a 
new TRiP submission could occur in parallel minimizing the scheduling impact. This upgrade is not 
essential for MINERz^A, but could provide enhanced performance. 

We also note that if the new submission fails, DO might need to use the existing chips. If the yield 
is 90%, as was found in a sampling of 100 chips, we would have about 2460 chips after satisfying the 
DO requirements, which is enough for 39000 channels. Hence, in the "worst-case" scenario of both DO 
needing the existing TRiP chips and also a lower yield than the sampling to date would suggest, we 
might have to make more TRiP chips with the existing design and masks. 

16.6 Parameters of the MINERz/A Detector 

MINERz^A combines the fine granularity of an electronic bubble chamber with the final state analyzing 
power of more traditional (but very fine grained) sampling calorimeter and muon magnetic spectrome- 
ters. To maintain the segmentation required to identify each final state particle in a low energy neutrino 
interaction and to accurately track final state photons for vr'^ reconstruction, the number of channels in 
MINERi^A must be large. To contain the produced final state particles, the mass of MINERzvA must be 
large. We attempt to break down the contributions to mass and channels by sub-detector in this section. 

TablefT6llists the total number of channels by sub-detector. Predictably, it is the granuarlity required 
in the plastic, Pb and Fe targets that dominates the channel count, with the downstream calorimeters, 
side calorimeters, the muon and the veto systems contributing 19%, 17%, 7% and 1% of the channels, 
respectively. As shown in Table [nj the situation is very different with the mass apportionment among 
the detectors where the OD and MR dominate the mass. 

The scintillator and optical system system of MINERi/A, though it pales in comparison to long 
baseline neutrino experiments like MINOS, is impressive on the scale of the CDF Plug calorimeter 
or CMS HCAL. MINERz^A will use 19.2 metric tons of extruded polystyrene scintillator, 93 km of 



Component 


Number 


Comments 


Channel 


37478 


WLS Fibers 


PMT boxes 


587 


includes 90 empty M-64 anodes 


Readout Token Rings 


37 


16 PMTs/ring 


VME Readout Cards 


10 


4 rings/card 


VME Slow Control Cards 


20 


30 PMTs/card 


VME Crates 


4 




VME PVIC Interface 


4 


one per crate 


DAQ PCs with PVIC, 


1 


data rate is 120kB/spill 


RAID system 







Table 15: Parts count for MINERi/A Electronics Design 
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Sub-Detector 


Channels in Inner Detector 


Channels in Outer Detector 


Active Target and Side ECAL 


15360 


5760 


US ECAL (Pb Target) 


3072 


1152 


US HCAL (Fe Target) 


1536 


576 


DS ECAL 


2560 


960 


DS HCAL 


2560 


960 


Veto 


426 


n/a 


MR/Toroid 


2556 


n/a 


Totals 


28070 


8408 



Table 16: Channel count by sub-detector 



Sub-Detector 


Mass (metric tons) 


Active Target 


6.1 


Side ECAL 


8.5 


US ECAL (Pb Target) 


3.5 


US HCAL (Fe Target) 


7.0 


OD Framing the Target Regions 


126.5 


DS ECAL 


19.8 


DS HCAL 


26.4 


Veto 


15.1 


MR/Toroid 


90.8 


Total 


302.1 



Table 17: Mass by sub-detector 
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wavelength-shifting fiber and another 46 km of clear fiber in optical cables between the detector and the 
587 M-64 multi-anode photomultipliers. 
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17 Cost and Schedule 



This section describes the cost and schedule associated with the construction of the MINERi/A detector. 

Given the relatively modest cost of the detector, we plan to largely fund the construction of the 
MINERz/A detector from a combination of university program research funds and special program 
funds indepedendent of the Fermilab budget. 

Portions of the project that would be, by necessity, managed and funded by Fermilab would include 
site outfitting and utilities (e.g., magnet and quiet power, cooling), crucial safety items for the NUMI 
hall that must be designed at Fermilab (e.g., low voltage distribution to the electronics, the magnet 
coils), and installation costs associated with bringing modules to the NUMI near hall. At the time of 
the submission of this proposal, we do not have complete evaluations of these costs. As discussed in 
Secti0n fT4.ll these costs have not been estimated. We are encouraged, however, to note that the utilities 
requirements of MINERi/A appear to be within the capacity of the NUMI near hall, and do not appear 
to require major infrastructure upgrades. We expect to update this document with a good estimate of 
these costs by the time of oral presentation to the PAC on December 12, 2003. 

17.1 Description and Summary of Costs 

The cost of MINERi/A is dominated by three major categories of expenses: external materials pur- 
chases, craft durable items and labor to assemble the active elements and absorber into modules. Each 
of these has its own appropriate costing methodology. 

For the large external equipment costs, the MAPMTs, the clear and WLS fiber and the metal plate 
to construct absorbers, we have contacted our preferred vendors directly to obtain quotes. For the pho- 
tosensors, we have shown that the stock specifications of the R7600U-00-M64 PMTs are adequate for 
our application. Similar phototubes are in wide use throughout the lab, and have performed reliably. 
Hamamatsu has provided a quote on our quantity with a three month delivery time, and we are inves- 
tigating cost savings that can be realized through a more efficient custom packaging suitable for our 
Cockcroft-Walton supplies. It is worth noting that another manufacturer, Burle Technologies, manu- 
factures a product which would likely meet our specifications with better channel gain uniformity, the 
Planacon 85011-501, and we plan to pursue this possibility as well. The clear and WLS fiber vendor, 
Kurary, again is a vendor with a long history at Fermilab, and they have provided similar fiber to MI- 
NOS, CMS, CDF, etc. We have also secured a quote on our quantities independent of the concentration 
of WLS dopant should we chose to reoptimze the dopant for our strip lengths. Finally, the costs of the 
absorber were provided by suppliers who have established relationships with the Rutgers Physics De- 
partment machine shop, and variations here would likely result only from movement in the bulk prices 
of the relevant materials. The machining costs have been estimated through the Rutgers shop which is 
ready to perform the work as needed. Because of the relative certainly of these costs, we allow rela- 
tively low contingencies for these items, ranging from 20% (MAPMT and fibers) to 30% (absorbers). 
We have not yet included F&A costs on most of our equipment purchases since we anticipate most 
of these purchases will be made through University fuUy-costed shops which will try to negotiate low 
F&A costs on these large, bulk purchases. 

The second category of costs come from the craft items which must be constructed to assemble the 
detector, including the front-end electronics and associated auxiliary systems, parts for the PMT/front- 
end housing and the extruded scintillator strips. Here the strategy was to identify similar components 
from the construction of the MINOS far or near detectors, or from the CMS HCAL or CDF Plug con- 
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struction, and to attempt to scale the actual project costs. For example, the clear fiber cable costs 
were scaled from the CMS HCAL project which was of similar scope with similar fiber and connector 
scheme. These costs already include actual labor and EDIA costs. Following this approach is very 
useful because one learns interesting and relevant facts about hidden cost drivers. For example, the 
front-end boards require approximately 5 kWatts of power over the entire system, and surprisingly, 
scaling from the MINOS far front-end costs, one calculates a low voltage power supply system cost in 
excess of $135,000, excluding cables. The reason turned out to be the special fire protection require- 
ments imposed by underground operation of power supplies, so at least in this case, we found a cost that 
would otherwise have been missed. For these projects, we assigned contingencies between 40-50% of 
total sub-project cost, based on our scaled estimates. In the case of the electronics, there is an additional 
contingency cost of $70,000 for the case where we have to re-submit the TRIP ASIC due to unexpected 
demand for these chips from DO. 

The final item is technical labor for component assembly and testing for procedures that have not 
yet been prototyped. These have been estimated based on assembly models from the CMS HCAL 
and MINOS far detector projects, and are generally more uncertain than other estimates because of 
differences in construction between CMS HCAL, MINOS and what we proposed in MINERz^A. We 
have assigned contingencies of 50% for these projects. Labor and EDIA costs which dominate here are 
based on FY2005 projected costs for technician and engineering staff on the CMS HCAL project at the 
University of Rochester. 

A summary of the costs is shown in Table (TSl The total project construction cost is estimated to be 
$3.96M, excluding the installation and hall utilities costs. Our calculated contingency, $1.54M, is 39% 
of the total cost. As previously noted, the M&S costs do not include F&A. 

A brief summary of what is included in each sub-project category follows. 

Extruded Scintillator: prototype and production extrusion dies; purchase of plastics; extrusion in the 
Lab 5 facility; Q/C and monitoring. 

Fiber and Glue: WLS and clear fiber (1.2mm, Yll 0-400 ppm, J-type, S-35), BC-600 epoxy. 

WLS Fiber Prep.: Design and construct gluing assembly; cut fibers, mirror one end and glue into 
scintillator, prepare fiber pigtail for connector. 

Optical Cables: Purchase connectors and test equipment; EDIA for fiber termination procedure and 
cable layout; bundle fibers into conduit; insert WLS and clear fibers into two pairs of connectors; 
polish ends; test for transmission and light tight. 

Absorbers: ECAL: purchase and machine Pb and stainless sheets, epoxy, stainless to Pb, ship to as- 
sembly site; OD: order strips pre-cut from vendor and ship to assembly site; HCAL: order partial 
plates pre-cut from vendor, weld and ship to assembly site; Coil: purchase Cu AWG4 wire for coil 
and fabricate bus bars; Plastic fiber router plate for OD: purchase polypropylene sheets, program 
and route groves on CNC router and ship to assembly site. 

Module Assembly: prototype procedures; laminate sub-planes of ID strips; connect inner OD frame; 
connect stainless stop to frame and attach to strongback; construct OD in layers; attach plastic 
routing plates; lay in ID strips and route fiber; add stainless retainer; layer in additional planes; 
join OD at outside layer; attach WLS bundle connector; prepare for delivery to experimental hall 

Photosensors: purchase R7600U-00-M64 PMTs. 
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MAMPT Testing: design and assemble test stand; test for specs. 

PMT Box and Optics: design and prototype PMT box; design and assemble testing station; construct 
box; add connectors; assemble internal optical system; mount front-end board and PMT socket; 
add PMT; add fiber bundle from connector to cookie and attach; attach internal cables to board; 
light tight and Q/C. 

Electronics and DAQ: prototype front-end TRiP design; design front-end board; design VME data 
board; purchase VME crates, controllers, PVIC interfaces and DAQ PC; TRiP checkout; produce, 
assemble and checkout front-end, VME data and slow controls boards; purchase LV system; 
purchase LDVS, slow controls and LV power cables. 

17.2 Schedule 

The MINERi^A collaboration has not yet produced a resource-loaded schedule for the experiment ca- 
pable of reliably predicting the schedule. We plan to present such a schedule at the PAC meeting on 
December 12, and to update this section when it is available. 

This having been said, the schedule driving elements for the experiment to be ready to be installed 
are three: construction of the front-end electronics, assembly of the detector modules, and construction 
of the PMT boxes. We will discuss each of these in turn. With the possible exception of the PMT boxes, 
we have high confidence that the result of the resource-loaded schedule will be to produce detector 
modules ready to install in the NUMI near hall approximately two years from the project start, assuming 
that the bulk of the project funds for M&S items can be expended at the front end of the project. 

Construction of custom electronics with an ASIC would usually be an overriding concern in such 
an aggressive schedule. However, as Section fT6 . 5 l explains . both the ASIC and the bulk of the front-end 
designs are being recycled from the DO fiber tracker upgrade to 132 ns bunch crossings. Design for a 
vertical slice test of a prototype front-end system has already begun, and we are confident that we can 
demonstrate success of this prototype front-end by 2004. The VME data board will require only minor 
modifications from existing designs, and the slow controls board is a stock design which will require 
no modification. Even with the earliest project start date of summer 2004, we would have a completed 
design of all boards by the end of 2004, and be finishing production in the middle of 2005. 

The assembly of the modules is a very complicated task because of the large number of channels 
in the detector, the complicated routing of fibers in the detector and the need to reduce support mass 
in the inner region of the detector. Furthermore, assembly of modules cannot begin until scintillator 
production, WLS fiber installation and absorber production are well underway. We have developed an 
assembly procedure and manpower assessment based on that procedure and the University of Rochester 
CMS HCAL experience that suggests that seven technician-years would be required to assemble the 53 
MINERi/A modules. It is aggressive, but perhaps plausible to attempt to fit those seven technician years 
into twelve or fifteen months, after completing a prototyping Q/C development phase of six months. Of 
the prerequisites for beginning module construction, it is most likely the start of significant scintillator 
production will most severely limit our ability to prototype and construct modules. We would expect to 
be able to begin scintillator production approximately four months after the project start, and therefore 
we conclude module production could be complete 22 - 25 months after project start. Installation of 
the modules in the near hall could proceed in parallel with the completion of the last modules. 

Finally, we are working to develop a complete model of the construction of the PMT box and 
associated optical components. With over 550 boxes to construct and assembly of complicated optical 
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connectors and a fiber bundle in a tight space, the quality control concerns are very non-trivial. We have 
sufficient experience within the collaboration from the design and construction of the MINOS "MUX 
boxes" and the CMS HCAL to address this problem and expect to have a confident assessment of the 
schedule and schedule risks associated with the PMT box by the time of the PAC presentation. 
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Cost (kUSD) 






M&S 


SWF 


EDIA 


Contingency 


Sub-Project 


(no F&A) 


(w/ F&A) 


(w/ F&A) 


(%) 


Extruded Scintillator 


151 


12 


30 


78 (40%) 


Fiber and Glue 


262 


n/a 


n/a 


52 (20%) 


WLS Fiber Prep. 


50 


104 


16 


85 (50%) 


Optical Cables 


77 


162 


11 


100 (40%) 


Absorbers 


310 


67 


32 


122 (30%) 


Module Assembly 


11 


473 


53 


268 (50%) 


Photosensors 


772 


n/a 


25 


159 (20%) 


MAPMT Testing 


6 


45 


n/a 


26 (50%) 


PMT Box and Optics 


278 


95 


51 


212 (50%) 


Electronics and DAQ 


628 


33 


206 


435 (50%) 


Totals 


2545 


990 


423 


1537 (39%) 



Table 18: Summary of MINERz/A detector costs in exclusive sub-project categories 
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A Cryogenic LH2 and LD2 Targets 



Some of the nuclear and fundamental particle physics described in this document may be dramatically 
improved by the inclusion of cryogenic liquid Hydrogen and/or Deuterium targets. The data from 
such targets would allow a detailed comparison with Jefferson lab experiments that are currently using 
Liquid Hydrogen and Deuterium as targets in electron and photon scattering experiments. A comparison 
with these experiments in a similar momentum transfer range with the high precision neutrino cross 
section measurements made possible by MINERi/A would provide for unprecedented studies of nucleon 
structure, particularly at large x where it has heretofore been very difficult. It is, for example, clear that 
the substantial uncertainties on parton distribution functions at large x, which are dominated by nuclear 
corrections and uncertainties involved in the flavor decomposition, would be removed. 

Those measurements described in this proposal which may become Umited by resolution would 
benefit greatly from the inclusion of a cryogenic target system. The interpretation of resonance data, for 
example, would no longer be complicated by uncertainties in nuclear binding and on shell extrapolation. 
This would allow direct comparison with the Sato and Lee pion cloud predictions without additional 
model systematics. The availabihty of a clean nucleon target would remove the complexity of the 
nuclear potential in heavy targets allowing the underlying physics of strange particle production, for 
instance, to be probed without interference and quasielastic studies would be greatly helped by the lack 
of intra-nuclear proton scattering. Furthur, comparison of the data from the liquid Hydrogen/Deuterium 
target with data gathered from interactions in the MINERi^A nuclear targets would be of great use in the 
determination of nuclear effects in neutrino interactions and would help to Umit even furthur systematic 
effects in the oscillation experiments. 

The cryogenic target itself would be small and compact. It could be installed upstream of the 
detector proper and would only require that the veto array be moved to cover it. There exists the 
possibility that the cryogenic target could be converted from a passive to an active target with the 
inclusion of CCD cameras to view the interactions; however, even considered only as a passive target, 
a high-statistics sample of neutrino interactions on liquid Hydrogen or Deuterium would be of a great 
benefit in the understanding of neutrino interactions in this relatively complicated few-GeV region. 
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B Off-Axis Running 



When the MINERz^A detector is located on the beamline axis of the NuMI beamline, it is exposed to 
a broad band of neutrino energies with a peak energy dependent on the momenta of the pions being 
focused by the horns. Because there are high energy mesons that travel through the holes of both 
horns, there is also a long tail of neutrino events with energies well above the peak energy. However, 
as the detector moves off the axis, the peak neutrino beam energy spectrum decreases and becomes 
much more narrowly distributed in energy, and the highest energy mesons are no longer pointing at 
the detector, essentially removing the "high energy tail", as shown in Figure |93] Note that at lOmrad 
away from the beamline axis the event rate is peaked at 2GeV. This is solely a result of the 2-body 
kinematics of the vr —>■ u^fi decay. Although the event rates are highest when the MINERi/A detector 
is on axis, the energy of the incoming neutrino is not known a priori and so by moving the MINERz^A 
detector off axis for a given running period the experiment can make measurements of cross sections 
in a more narrowband beam. This is particularly useful for neutral current measurements, where the 
total incoming neutrino energy cannot be reconstructed because of the loss of the outgoing neutrino. 
Due to the intensity of the NuMI beamline, MINERz^A can collect appreciable statistics for precision 
measurements of low energy neutral current processes when running for short periods of time off the 
beam axis. 
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Figure 93: Distribution of off-axis events that are available at different distances off-axis at the distance 
from the target of the MINOS near detector in the LE beam. In order of the peak from left to right, the 
curves represent event rates 20 m off-axis, 10 m off-axis, 5 m off-axis and on-axis, for comparison. 

The NuMI underground complex itself was excavated primarily by a 21.5 foot diameter tunnel 
boring machine (TBM), and because of this there are large sections of excavated regions underground 
which will be unused once the MINOS near detector is in place. These regions are located anywhere 
from to 20mrad off the NuMI beamline axis, and could house a future off axis near detector. The Off 
Axis experiment is likely to place a near detector in these drifts to be able to predict the appearance 
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backgrounds at a far off axis detector. 

Figure shows three possible locations for the MINERi/A detector to be placed for off axis run- 
ning. We discuss from downstream to upstream, the advantages and disadvantages of each location, 
keeping in mind that none of these locations would require any additional excavation. 




Figure 94: Possible sites for off axis running in the NuMI Underground area 

The most downstream site, site 1 in figure |94j is the easiest site to use, as it is the closest to the 
MINOS near hall. It is just downstream of the shaft, and the floor is flat between the base of the MINOS 
shaft and the MINOS near hall. This location views off axis beams anywhere from 5 to lOmrad off 
the NuMI Axis. The drift itself has an access tunnel on the east side for emergency personnel egress, 
and some cable tray and utilities on the west side, but there is a region in the middle of the drift which 
measures 4.5m wide by 6m tall which is currently vacant, as shown in figure |95] The neutrino energy 
spectrum in this hall would be anywhere from 1.5 to 3 GeV, depending on the location, since at the 
downstream end of the access tunnel the tunnel is nearly on axis, and at the upstream end of the tunnel, 
near the shaft, the tunnel is about 10m off axis. 




Figure 95: Off-axis drift cross-section for site 1. 

Moving upstream, the next site is just upstream of the MINOS shaft (called site 2 in figure 1941. This 
location is also relatively easy to get to since it is near the MINOS shaft and therefore close to utilities, 
but the floor in this region has a 9% slope. The available cross section in this area is also wide, at least 
as wide as the area for site 1 , but some space must be left for access to points upstream, since there is 
no longer an independent egress tunnel as with site 1. Here the mean neutrino energy is about 1.5GeV. 
Finally, the most remote site is located in a drift that was created when the tunnel boring machine had 



178 



to change angles between excavating for the decay volume itself and excavating for the access tunnel to 
the downstream areas. This site is substantially closer to the decay pipe, located about 725m from the 
NuMl target. This site has the widest cross cross section available since it is located in a "dead end", 
and would allow off axis distances from 5 to 15m off axis, which correspond to between 8 and 20mrad 
off axis angles. To run in this location the detector would have to be moved up the 9% sloped floor 
about 200m from the base of the MINOS shaft, and then re-assembled in this new hall, which also has 
a floor with a 9% slope. There are always several meters of earth between the NuMl hadron absorber 
and the detector, but neutron radiation issues would be worse here than in the other two sites. However, 
this site's weakness is also its strength, in that a detector in this location would have the least amount of 
interference with the MINOS experiment during construction. The closer location would also provide a 
significantly higher event rate. The downstream portion of this site has utilities for the NuMl absorber, 
but there are about 10 meters of cleared space upstream of those utilities where the MlNERz/A detector 
could be placed. 
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Party 

Addendum to the MINERi/AProposal 



C Updates to the Physics Case 

(submitted to the Fermilab PAC 30 March 2004) 

C.l Introduction 

As detailed in the MINERz/A proposal, the imminent completion of the NUMI beamline offers the 
particle and nuclear physics communities a new opportunity to study neutrino interactions in an envi- 
ronment unprecedented in granularity of detectors and event rate. The construction of MINERi^A, with 
its fully active core target, would allow a wide variety of measurements in neutrino interaction physics, 
which both support future and current neutrino oscillation efforts, and are interesting for their own sake. 

Since the proposal was submitted to the FNAL FAC at the December 2003 meeting, the MINERz^A 
collaboration has continued to make progress in understanding the ultimate capabilities of the experi- 
ment. This addendum serves to document these continuing studies for the benefit of the PAC as we seek 
approval of the experiment. 

This document is broken down by physics topic, and each section contains a brief summary of the 
physics goals, the status of our understanding at the time of the proposal, and then documents improved 
understanding, since the proposal, in more detail. The topics considered in this addendum are: 



• Quasi-Elastic Cross-Sections and Form Factors 



• Coherent Pion Production 



• Physics Opportunities in the Resonance Production Region 



• Nuclear Effects in Neutrino Scattering 



• The Impact of MINERz^A on Oscillation Experiments 



This document is intended to be read as a supplement to the updated MINERi/A proposal in the main 
body of this document. This document, along with other documentation of the status of the MINERz^A 
experiment, is available from the MINERi/A collaboration web page, h ttp : / / www . pas . rochester 
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C.2 Quasi-Elastic Scattering 



At the lowest neutrino energies relevant for future long baseline efforts, it is the quasi-elastic scattering 
that dominates the charged current interaction rate. As outlined in the MINERi/A proposal"*, there 
are some interesting physics issues to be addressed in quasi-elastic neutrino scattering. The first of 
these is understanding the impact of nuclear effects on the quasi-elastic kinematics at low which 
dominates the signal rate for oscillation experiments. The second is understanding the axial form-factor 
of the proton at high which will contribute to the blossoming body of new measurements on high 

nucleon form-factors where some surprises have already been seen in charged-lepton scattering (see 
main proposal above). Since the proposal, we have made significant progress in simulating our expected 
analysis of quasi-elastics, focusing on the important issues of maintaining high efficiency at low and 
low backgrounds at high Q"^. 

In un — > iJL~p, the outgoing proton carries a kinetic energy that is approximately Q"^ /2Mn- So for 
low Q^, the challenge is identifying events with a very soft recoil proton; for high Q^, this proton is 
high energy and may interact in the detector, making particle identification more challenging. The main 
strategies of the current analysis are: 

• At low Q^, accept quasi-elastic candidates with a single (muon) track, and discriminate from 
background by requiring low activity in the remainder of the detector 

• At high Q^, reconstruct both the proton and the muon, and require kinematic consistency with 
X = 1 and = 

Simple cuts deriving from these ideas allow for reasonable efficiency with good purity, even at high Q^. 
C.2.1 Details of Quasi-elastic reconstruction 

The analysis uses the NEUGEN generation and the hit level MINERz^A simulation and tracking package 
in order to simulate signal selection and background processes. 

The initial event identification proceeds by requiring one or two tracks in the active target. One 
of these tracks must be long range (400 glcrc?) and is the putative muon. If a second track forms a 
vertex with this track, it is assumed to be the proton. No other tracks are allowed to be connected with 
this event vertex. The muon track momentum is reconstructed with a fractional uncertainty of between 
10-20%. 

In the low case, the proton track (if found) would be effectively required to lose energy by range 
since only a limited amount of detector activity not associated with the primary tracks is allowed by the 
event selection. We attempt to recover some of the lost efficiency at higher due to this cut by allowing 
hits on tracks near the proton track to be associated with the proton track itself. Figure |96l shows the 
fraction of hits not associated with the lepton or proton in the quasi-elastic events and in expected 
background processes. For higher events a similar requirement could in principle be applied, but it 
is not particularly effective nor efficient. 

The energy of the proton for the high case (where the proton almost always interacts) is re- 
constructed calorimetrically with an expected fractional energy resolution that is well parameterized by 



Although muons are identified by requiring a single track with a long range in the detector, no 
attempt is made in the analysis to improve particle identification by requiring dE / dx consistent with 

■*see the discussion in Chapter 6 of the main proposal 
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Figure 96: The fraction of hits associated with the muon and the proton tracks in quasi-elastic can- 
didates. The events for the plot may have have one or two vertex tracks, pass additional kinematic 
requirement and are required to have 0.1 GeV^ < < 1 GeV^. 
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the muon or proton tracks. This requirement is expected to be particularly effective for protons of 
0(1) GeV momentum^, and such a requirement can be optimized to improve the analysis in the future. 
In addition, it may be possible to improve the efficiency by allowing a lower range muon with a dE / dx 
requirement without sacrificing purity. 

If a quasi-elastic interaction is assumed, one can reconstruct the event kinematics from only the 
momentum and direction of the final state /x. Neglecting the binding energy of the final state proton, 

^QE ^ MnE, - 4 

Mn - E^+p^ cos Of, ' 

If a proton track is required and its angle and energy are measured, one can additionally require consis- 
tency with the quasi-elastic hypothesis. Two constraints are possible, one on the x of the reconstructed 
interaction and one on the pT of the observed final state. 

If the interaction is truly quasi-elastic, then x = 1, and therefore = 2Mni' where u = E^ad - 
Mnucieon^ and E^ad is the energy of the hadronic final state. In this analysis, we test this by comparing 

reconstructed from the lepton kinematics under the quasi-elastic hypothesis to 2Mni^ and forming 
{Q^ — 2Mni^)/ error where the dominant part of the calculated error for this term comes from the 
smearing of hadronic final state energy. Figure |^ shows this difference significance for two track 
quasi-elastic candidates with observed 1 GeV^ < < 3 GeV^, for quasi-elastic, resonance and DIS 
events. Note that this cut can be applied without identifying a proton track if the visible energy, less the 
muon energy, is assumed to be v. 

The significance {x) cut does not use information on the proton direction, and so we impose 
a second kinematic cut on the of the final state relative to the incoming neutrino direction. This 
selection requires that a proton track is identified and we cut on the significance of the difference from 
Pt = 0. We impose a cut of p^ /error <2 except for > 3 GeV^, for which the cut is 3. Note also 
that if we impose a p^ cut first, the difference cut still improves the result, i.e. both cuts are needed. 

In summary, the selection requirements for quasi-elastic candidates are: 

• One or two tracks for < 1 GeV'^ and two tracks for > 1 GeV"^. 

• One track must have 400 g/cm^ range (muon). 

• {Ql - 2Mu)/ (error) < 2. 

• Pt/ (error) < 2 for < 3 GeV"^ and px/ (error) < 3 for > 3 GeV"^. 

• Hit fraction associated with muon and proton > 0.9, for < 0.5 GeV'^, or > 0.85, for 

0.5 Gey2 < q2 ^ 1.0 GeV^. 

C.2.2 Results 

TablefTOlshows the efficiency and purity of the quasi-elastic sample for different bins after each cut. 
Using the calculated efficiency and purity, we have updated the uncertainties on Fa derived in the main 
proposal which did not include efficiency or background effects. 

Figure |98] shows the extracted values and errors on Fa in bins of from a sample of simulated 
quasi-elastic interactions in the MINERz^A active carbon target, for a four-year exposure in the NuMI 

'see Section 15.5.5 of the main proposal 
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Figure 97: The significance of the difference between from the quasi-elastic hypothesis and from 
the final state energy 



Q2 bin 


Effic Purity 


{Q'f, - 2Mu)/en: 
Effic Purity 


Pt /err 
Effic Purity 


Hits 
Effic Purity 


0.1-0.5 
0.5-1 

1- 2 

2- 3 
3 - 10 


0.926 0.246 
0.775 0.199 
0.600 0.199 
0.456 0.146 
0.689 0.123 


0.918 0.442 
0.765 0.410 
0.541 0.416 
0.400 0.375 
0.600 0.310 


0.866 0.559 
0.624 0.486 
0.397 0.555 
0.344 0.554 
0.467 0.420 


0.775 0.842 
0.528 0.685 
0.338 0.598 
0.278 0.676 
0.311 0.700 



Table 19: Efficiency and purity in bins for quasi-elastic candidates 
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QE scattering, u^, Fa(Q^) 




Figure 98: Estimation of Fa from a sample of Monte Carlo neutrino quasi-elastic events recorded 
in the MINERz^A active carbon target. Here, a pure dipole form for Fa is assumed, with Ma = 
1.014 GeV/c^. The simulated sample and error bars correspond to four years of NuMI running. Also 
shown is Fa extracted from deuterium bubble chamber experiments using the da / dq^ from the papers 
of FNAL 1983 E?! BNL 1981 El, and ANL 1982 El 
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Figure 99: Extracted ratio Fa/ FAi^ipole) from MINERi/A active target (4 NuMI-years) under two 
scenarios for the true Fa- (Top) ^^/^^(Dipole) is G^/G^(dipole) from polarization transfer mea- 
surements. (Bottom) F^/F^CDipole) is G^/G^(dipole) from Rosenbluth separation technique. Also 
shown is Fa extracted from deuterium bubble chamber experiments using the da / dq^ from the papers 
of FNAL 1983 |13 BNL 1981 01, and ANL 1982 I25j 
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beam. Figure l99l shows these results as a ratio of F^/F^(Dipole), demonstrating MINERz^A's ability 
to distinguish between different models of Fa- Note that resolution effects are still not included in this 
extraction of Fa', however, the typical resolution for quasi-elastic events at high is ~ 0.2 GeV^ 
which is smaller than the bin size. 
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C.3 Coherent Neutrino-Nucleus Scattering 

With its high statistics, fine segmentation, excellent tracking, good particle ID, and range of nuclear 
targets, the MINERi/A experiment will be able to obtain data samples of coherent interactions several 
orders of magnitude larger than those published from previous experiments. The study of neutrino co- 
herent scattering is an excellent tool for exploring the 'hadronic' nature of the weak current, particularly 
the axial component. 

The capabilities of MINERi^A for coherent scattering are compelling for a number of reasons: 

1. Coherent vr production is largely a neutrino-specific process, as the coupling is dominated by axial 
current. In the Rein-Seghal model ll93l . for instance, the coupling is entirely axial. Therefore 
there is no electron-scattering analogue to constrain models for this process. 

2. Existing measurements in the energy range of the NuMI beam are quite poor. There is only one 
experiment (SKAT) with data on charged current coherent production below of 10 GeV, and 
it only had 59 events. 

3. The charged current is an extremely clean measurement in MINERz^A because distinct muon and 
pion tracks will be clearly visible and will allow precise measurement of the interaction point and 
kinematic quantities. 

4. With nuclear targets spanning A=12 to A=207, MINERi/A will be the first experiment to measure 
the A-dependence of the coherent cross section. 

5. An improved knowledge of coherent cross sections will be crucial for future oscillation experi- 
ments, as the neutral current reaction producing a single vr*^ is one of the main backgrounds for 
subdominant — > mixing searches. 

6. The tt'^ reconstruction capabilities of MINERi/A make it possible to statistically separate coherent 
from non-coherent vr'^ production to directly determine this background and to check the expected 
relationship to the charged-current process. 

MINERi/A selects coherent events in both the charged and neutral current by relying on the distinct 
kinematics of the events. Because the coherence condition requires that the nucleus remain intact, the 
process is tagged by low-energy transfer (|t|) to the nuclear system which is reconstructed by 

-\t\=-{q- P.? = {^i{Ei - ph? - (Mpi)?- (64) 

Candidate events are generally selected as coherent by requiring a low final state multiplicity and by 
requiring kinematics consistent with low \t\. 

The MINERz^A event selection is described in detail in the proposal^. The efficiency and purity of 
the selection in the charged current case is summarized in Table IC.31 and the expected uncertainties in 
the final cross-section as a function of energy is shown in Figure 135) 

C.3.1 A-Dependence of the Cross Section 

As noted above, there are no existing measurements of the coherent cross-sections on light nuclei (e.g., 
H, C, O). These measurements are important for planned appearance experiments. Current predic- 
tions must rely on extrapolation in A in addition to extrapolation to lower energies. 

^Chapter 8 of the main proposal 
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Cut 


Signal Sample 


Background Sample 




5000 


10000 


2 Charged Tracks 


3856 


3693 


Track Identification 


3124 


3360 


TT^/neutron Energy 


3124 


1744 


Track Separation 


2420 


500 


x<0.2 


2223 


100 


t<0.2 


2223 


19 


Ptt < 600 MeV 


1721 


12 



Table 20: Analysis cuts to isolate a sample of coherent interactions. The cuts are described in the text. 
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Figure 100: Coherent cross-sections for 5 GeV neutrinos as a function of atomic number. Solid curve 
is the prediction of Rein-Seghal, sohd circles are the predictions of Paschos-Kartavtsev for carbon, iron 
and lead. Crosses indicate the expected measurement errors for minerva assuming the measured cross 
section is that of Rein-Seghal. The shaded band indicates the region explored in previous experiments, 
primarily from measurements on aluminum (A=27), neon (A=20), and freon (<A>=30), all performed 
at higher energies than are relevant for MINERi/A or future oscillation experiments. 
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Figure 1 1001 shows the predicted A-dependence of the charged current coherent cross section from 
Rein-Seghal and Paschos-Kartavtsev I93II96I models, and the expected measurement errors from MINERi/A. 
As this analysis indicates, the high statistics and large dynamic range in A of the MINERz^A experi- 
ment will make possible detailed examinations of the coherently produced meson-nucleus interaction. 
Although K2K and MiniBooNE are expected to significantly improve the knowledge of coherent vr pro- 
duction on CHn at very low energies, only MINERz^A will be able to reach this level of precision and 
to perform this systematic examination of the A dependence of the cross-section. 
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C.4 Measuring Nuclear Effects with the MINERi^A Detector 



As indicated in our proposal, to study nuclear effects in MINERz^A, carbon, iron and lead targets will be 
installed upstream of the pure scintillator active detector. The currently preferred configuration involves 
a total of 6 planes, with each plane divided transversely into C, Fe and Pb wedges. As one proceeds from 
upstream to downstream, the C, Fe and Pb targets exchange (rotate) positions. As always, a scintillator 
module of four views (x,u,x,v) separates each of the planes. The total mass is over 1 ton of Fe and 
Pb and somewhat over 0.5 ton of C. Since the pure scintillator active detector essentially acts as an 
additional 3-ton carbon target (CH), the pure graphite (C) target is mainly to check for consistency. For 
the standard four-year run described in the proposal, MINERzvA would collect 1 M events on Fe and 
Pb, 600 k events on C as well as 2.8 M events on the scintillator within the fiducial volume. In this 
section we give more experimental background than was given in the proposal', as well as describe 
more specifically the analysis technique that will be used to measure nuclear effects. 
MINERz/A's goals in measuring nuclear effects can be summarized as follows: 

• measure final-state multiplicities, and hence absorption probabilities, as a function of A with 
initial u; 

• measure the visible hadron energy distribution as a function of target nucleus to determine the 
relative energy loss due to final state interactions (FSI); 

• investigate if the correction factors for observed multiplicity and hadron energy are a function of 
the muon kinematics for a more directed application of nuclear effect corrections. 

• measure a{xBj) for each nuclear target to compare x^j-dependent nuclear effects measured with 
both u and charged lepton beams. 

• If sufficient V running is available, measure the nuclear effects on F2(x,Q^) and xF3(x,Q^) to 
determine whether sea and valence quarks are affected differently by the nuclear environment. 

C.4.1 Pion Absorption Effects in Neutrino Interactions 

Interactions of few GeV neutrinos in nuclei easily produce resonances which decay to pions. Any at- 
tempt to reconstruct the incident neutrino energy based on the total observed energy must take into 
account the interactions of the pions in both the interaction nucleus and the detector. Current neutrino 
interaction Monte Carlos (such as INTRANUKE |209 1) handle pion interactions crudely and have gen- 
erally not yet incorporated the latest knowledge of pion interactions. We would like to summarize here 
the cuiTcnt knowledge of pion interactions and discuss plans for using MINERz^A to better account for 
pion interactions for neutiino-nucleus interactions. 

Our concern here is mainly with pions in the energy range of 100 to 500 MeV, where the interaction 
cross sections are the highest. In this range the pion-nucleon cross section is dominated by the very 
strong A(1232) resonance. The A is a fairly broad (about 100 MeV) resonance, and the pion-nucleon 
cross section reflects this, with a peak near 200 MeV pion energy which drops quickly above and 
below this. The pion nucleus cross section exhibits a similar behavior, with a less pronounced drop- 
off at higher energy. The charged pion nucleus cross section has four important components in the 
intermediate energy range: elastic scattering (nucleus left in ground state), inelastic scattering (nucleus 

^see Chapter 12 of the main proposal 
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left in excited state or nucleon knocked out), true absorption (no pion in the final state), and single 
charge exchange (neutral pion in the final state). 

Neutrino detectors are mainly iron (absorber), oxygen (water) and carbon (scintillator). The total 
pi-carbon cross section is 600 mb, with elastic and inelastic cross sections about 200 mb each, and 
absorption about 160 mb. The total pi-iron cross section is about 1700 mb, with elastic and absorption 
about 600 mb each, and inelastic about 400 mb. Cross sections for positive and negative pions are 
nearly the same because nuclei contain about the same number of protons and neutrons. These very 
large cross sections mean that essentially all pions will undergo some nuclear reaction, many within 
the interaction nucleus, and in most circumstances nearly all will be absorbed rather than stopping or 
exiting the detector. The absorption probability within the interaction nucleus is order 30% while the 
absorption probabilities in the detector are about 1%/cm in scintillator and 4%/cm in iron. Figure (ToTI 
1173 1 shows absorption cross sections for various nuclei as a function of pion energy. 

In elastic and most inelastic reactions the scattered pion will not, because of its light mass, lose 
much energy. However, absorbed pions will lose all of their kinetic and mass energy. Most of that 
energy will go into nucleons. We want to discuss here what happens to that energy. 

Pion absorption cannot occur on a single nucleon due to energy and momentum conservation. The 
simplest absorption mechanism is on two nucleons. Because absorption appears to proceed mainly 
through the N — A intermediate states, an isospin zero (np) pair is the primary candidate. Such an 
absorption for a positive pion would give two energetic protons whose kinetic energy nearly equaled 
the total pion energy. However, early studies of pion absorption found that was not the most probable 
mechanism. 

In the 1990's two large solid angle detectors, the LAMPF BGO Ball and the PSI LADS detector, 
were built to study pion absorption. Both detectors had large solid angles (both more than 90% of the 
full solid angle) and low proton thresholds (about 20 MeV for each). The LADS detector also had 
reasonable neutron detection efficiency and energy measurement. The somewhat surprising result from 
both detectors was that pion absorption was dominated by three body absorption II207I . For positive 
pions, the absorption on a pnn triplet (leading to a ppn final state) was the most common. This was 
observed even in ^He. The absorption in heavier nuclei also appears to proceed mainly through a 
three-body mechanism, although increased initial state interactions (pion re-scattering) and final state 
interactions (nucleon re-scattering) result in four to five nucleons being emitted. Typically the final 
state contains more neutrons than protons. The absorption process, which is still not well understood 
theoretically, largely fills the available phase space thus giving a wide range of nucleon energies with 
little angular dependence. 

Because much of the energy is in neutrons, the observed energy is well below the total pion energy. 
Fi sure fT02l and Figure fT03] 1 2081 show missing energy (total pion energy minus the total proton kinetic 
energy) for absorption of 250-500 MeV positive pions on ^^C and ^^Ni. As can be seen, even in carbon 
more than half the energy is lost to unobserved particles, a fraction which increases with pion energy 
and with A. 

The situation is of course worse for negative pions. Charge symmetry would indicate that the 
primary absorption should be on a ppn triplet leading to a pnn final state. In this case, most of the pion 
energy would be in neutrons, and hence not directly observed. However, if the interaction vertex and 
one proton energy is known, and the angles of the outgoing neutrons are known, the total energy of the 
three nucleons can be estimated. Monte Carlo studies with realistic absorption models will be needed 
to determine the accuracies of such estimates. 

Although neutral pions escaping the nucleus will decay, usually to two photons, the mean distance 
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Figure 102: The missing energy (total pion energy minus the total proton kinetic energy) for the ab- 
sorption of 250-500 MeV pions on carbon. 
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Figure 103: The missing energy (total pion energy minus the total proton kinetic energy) for the ab- 
sorption of 250-500 MeV pions on nickel. 
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Figure 104: Probability for the outgoing proton to escape the nucleus as a function of Q^. 

traveled before decay is a few nanometers, much greater than the size of the nucleus. Thus the absorp- 
tion of neutral pions in the interaction nucleus must also be taken into account in any study of resonance 
production. 

We have begun studies with INTRANUKE to determine the sensitivity to the probability of pion 
absorption in the interaction nucleus. We are currently in the process of modifying Monte Carlo routines 
to treat pion absorption more realistically. Because MINERi^A will also have good neutron detection 
capability, we expect to be able to substantially improve the determination of the incident neutrino 
energy. 

We also wish to note that there are essentially no measurements of pion absorption above 500 MeV. 
The fine spatial resolution and full solid angle detection capability of MINERi/A will allow us to study 
these interactions, especially in carbon. 

C.4.2 Nuclear Transparency in Neutrino Interactions 

A second nuclear interaction process which affects the observed final state energy is the final state 
interaction of a nucleon in the struck nucleus. An outgoing nucleon has a substantial probability of 
interacting in the nucleus. These probabilities have been measured, most recently at JLab, with some 
precision. The experiments used (e, e'p) coincidence reactions. The cross section for finding the scat- 
tered electron in the quasi-elastic peak was compared to the cross section for finding the coincident 
proton. A summary of the results are shown in Fig. 11041 

In contrast to pion absorption, there is little available information on what happens to the scattered 
nucleon. Of course, most either scatter from a single nucleon quasi-elastically or produce a pion (for 
protons above 600 MeV). Improving Monte Carlo routines to model this interaction should allow us 
to better estimate the total final state energy. As with pion absorption, the good resolution, neutron 
detection capability, and full solid angle coverage of MINERzvA should allow us to experimentally 
determine the actual final states and constrain the Monte Carlo routines. 
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C.4.3 Proposed Experimental Analysis 

The NEUGEN Monte Carlo has been used to study the sensitivity of the MlNERi/A experiment to nu- 
clear effects. The nuclear effects in the NEUGEN Monte Carlo are controlled by the INTRANUKE 
processor. This processor incorporates a probability for pion absorption based on earlier eletropro- 
duction absorption studies and lower-statistics Ne/if2 neutrino bubble chamber data. The observed 
phenomena of hadron formation length, which increases the transparency, is incorporated as well. The 
particular model used for pion absorption, which is currently being improved and updated, assumes that 
the absorption process eliminates a pion and the resulting nucleons are themselves either absorbed in 
the nucleus or are too low in energy to be observed in the detector. 

To determine the sensitivity of MlNERz^A measurements to the predictions of this model, the as- 
sumed probability for pion absorption in INTRANUKE has been increased by three standard deviations 
and then decreased by the same amount. The multiplicity and a very crude estimate of the visible 
hadron energy have been examined under these extreme conditions. In the next series of figures the 
predicted asymmetry in the multiplicity and visible hadron energy are shown. Currently we have only 
a very small sample of 2500 generated events available for this analysis. The asymmetry is defined as 
the percentage change under these extreme assumptions. That is, the bin contents at plus three stan- 
dard deviations minus the bin contents at minus three standard divided by bin contents at minus three 
standard deviations. Figure fTOSl shows the predicted change in the multiplicity distributions for carbon, 
while Fi gure IT06l shows the same distribution for iron. Both figures are consistent with the model that 
the decrease in observed multiplicity is caused by an increase in absorption probability and therefore 
the effect should be stronger in Fe than in C. The final determination of the visible hadron energy will 
be an involved process for this experiment. For now, we are using the most primitive estimate of this 
quantity, namely an uncorrected version derived from the total light output of the hadron shower. In 
the data analysis this will be refined for example, through measurements of stopping/decaying particles. 
With this crude estimate, the change in hadron energy for iron is shown in Fi gure (ToTl and for lead and 
Figure 11081 There is a dramatic change, even with these low simulation statistics, in the lowest energy 
bin. MlNERz^A will collect several times these statistics and should be capable of measuring this effect 
at even higher hadron energy. 

Since the incoming neutrino energy is not a priori known, the measured muon kinematics will be 
tested as a basis to compare characteristics of the visible hadron shower across nuclear targets and to 
determine whether the nuclear-effects correction-factor is a function of the observed muon kinemat- 
ics. The muon is relatively free from nuclear dependent effects and serves well as an A-independent 
normalization. For example, the quantity: 

Q' = 4.E^sir?{e/2) (65) 

is representative of the square of the 4-momentum transfer to the nucleon or quark, weighted (inversely) 
by Ey. This quantity then reflects the energy-momentum transferred to the hadronic vertex. The dis- 
tribution of events in this quantity are peaked toward low Q' with half the events below Q' = \ .0. The 
Monte Carlo statistics are still too low to make any conclusions on this study at this point, however the 
trend is encouraging and the study will continue with increased statistics. 
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Figure 105: The fractional change in multiplicity distributions between the two values assumed for pion 
absorption on carbon described in the text. 
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Figure 106: The fractional change in multiplicity distributions between the two values assumed for pion 
absorption on iron discussed in the text. 
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Figure 107: The fractional change in the visible hadron energy distributions between the two values of 
pion absorption on iron discussed in the text. 
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Figure 108: The fractional change in the visible hadron energy distributions between the two extremes 
in pion absorption on lead discussed in the text. 
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C.5 Resonance Production at MINERz/A 



Resonance production is an important issue both for future oscillation experiments, where resonant final 
states comprise much of both the signal and backgrounds, and for MINERi/A itself where the compar- 
ison of resonance production in electron and neutrino scattering provides an important motivation for 
our physics program in the low energy beam at NuMl. For the latter studies, a joint effort has been 
launched between the neutrino and electron communities to better utilize the available electron scatter- 
ing data for improved neutrino cross section modeling, considering nuclear effects as well as production 
mechanisms. This collaboration was recently approved by the Jefferson Lab Program Advisory Com- 
mittee to measure low , separated structure functions in the resonance region on a variety of targets of 
interest to the neutrino community (JLab Hall C Experiment E04-001). Many details of the MlNERz/A 
program in resonance production are given in the proposal^. 

Broadly, the resonance production measurements will focus on two areas. Fundamentally, MlNERi^A 
will measure pion production cross sections on a variety of nuclear targets nuclear targets, carbon and 
heavier. Secondly, the understanding of nuclear effects gained from these measurements will allow 
us to "extrapolate" the measurements to the free nucleon. Comparing these measurements to electron 
scattering will allow a better understanding of the axial component of neutrino resonance production. 
As nuclear effects will be important in all the measurements, we review below approaches to under- 
standing the nuclear effects. This is followed by a discussion of ongoing studies of event and particle 
identification techniques relevant to single pion production. 

C.5.1 Nuclear Corrections 

The interaction of neutrinos in nuclei produce secondary particles which propagate though the nucleus. 
Among them are protons produced in quasi-elastic scattering and pions produced in the resonance 
region. These effects are usually accounted for in Monte Carlo programs, such as the MINERi/A 
simulation, and in an analytic method described here. 
The propagation is viewed as a two step process: 

• A proton (QE) or a pion is produced on a bound proton or neutron, corrected for Pauli blocking 
and Fermi motion. 

• The produced particle propagates through the nucleus performing a random walk and it may 
exchange its charge in each interaction. 

The assumptions allow a general description of the problem based on charge symmetry and isospin 
conservation. The main property is a factorization of production and subsequent propagation. For the 
produced pions for instance we can write. 



where the subcript i denotes the number of pions produced initially at the neutrino interaction and, 
similarly, the subscript / denotes the number of pions emerging from the nucleus. The matrix M has a 

*see chapters 7, 10 and 13 of the main proposal 




/ A^i(vr+) 
M Nii-K^) 



) 



(66) 
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special form, which follows from isospin symmetry. 



M 



f l-c-d 
d 
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l-2d 
d 
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l-c-d 



c 



] 



(67) 



The parameter A contains absorption of the pion and the Pauli blocking in the pion interactions. The 
parameters c and d involve the effects of charge exchange in the re-scattering. The formahsm can be 
used in two ways: 

• To observe the yields in a specific nucleus and compare them to production of pions on free pro- 
tons and neutrons. Simple algebra relates the parameters to the observations. Thus determining 
the parameters in electroproduction allows one to use them in charge and neutral current neutrino 
reactions. 

• To calculate A, c and d theoretically by solving the random waUc problem. The numbers produced 
by the two methods are compatible with each other. 

More results are expected from the present experiments and MINERi/A will contribute on this topic. 
The method is needed to interpret the data and to determine oscillation parameters. 

As an illustration, we show in Figure 11091 the yields of pions in electroproduction and neutrino 
experiments on an ^^O target. The curves in the upper graph show the yields for vr'^s. The dotted curve 
shows the production on free isoscalar target,i.e. l/2(n + p). The solid curve shows the reduction 
produced on an oxygen target, where the change is substantial. The lower graph shows yields of vr+s 
where the nuclear corrections are much smaller. 

The MINERz^A monte carlo simulation, which uses the NEUGEN neutrino event generator also 
includes a type of random-walk simulation of pion re-interactions. These pion re-scatterings are im- 
plemented using the INTRANUKE 1 209 1 pion cascade model. With this model, the correlation of pion 
momenta before and after final state interactions can be studied. For ^^C, (Figure fTTOt . while some 
pions do scatter and loose momentum, many retain their original momentum. This correlation sug- 
gests that pion production data on carbon will be able to constrain free nucleon resonance production 
cross sections. While these correlations are not observable, MINERi/A will be able to test the pion 
re-scattering models by studying the A dependence of pion yields and distributions. 

C.5.2 Event Identification and Particle ID 

In electron scattering, resonance excitation spectra can be measured by detecting just the electron. 
However, in neutrino scattering, W and must be reconstructed using energies of the final state 
hadrons. Since the proposal, we have begun a program of studying resonance reconstruction in the 
MINERi/A detector using using our hit-level Monte Carlo simulation and the NEUGEN neutrino event 
generator. The techniques initially being studied are topological cuts to identify overall reaction type, 
and particle identification of individual tracks. 

We have performed a first analysis of single pion production detection efficiencies using simple 
topological cuts. The event selection here is very simple; tracks in the fully active target pointing to the 
vertex are counted and are identified as charged or neutral by the distance between the vertex and the 
first hit of a track. Tracks with a first hit less than 5 cm from the vertex are labeled as charged, while 
tracks further from the vertex are labeled as neutral (photons from vr'^ decay or neutrons). 
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Figure 109: Dotted lines show electroproduction cross sections for tt and tt"*" on an °0 target, with no 
final state interactions. The soUd curves show these cross sections after the nuclear effects of pion re- 
scattering and absorption, according to the model of Paschos et.al.[] The final state interactions greatly 
attenuate the vr'^ yield. However, charge exchange from the larger vr*^ production cross section serves to 
roughly balance the losses for the tt"*" channel, keeping the the vr^ yield about the same before and after 
pion reinteractions. Except for the overall scale of the cross section, these same results are predicted for 
neutrino neutral current reactions. 
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Figure 110: Correlation of pion energies before and after intranuclear re-scattering as implement in 
NEUGEN for the MlNERi/A simulation. Events are generated within the central plastic region (carbon 
and hydrogen scatterers) of the detector. Each plot shows on the vertical axis, the pion energy after it 
exits the nucleus versus the energy at the point that it is produced inside the nucleus. The first row of 
figures if for CC interactions while the second if for NC scattering. The left column is neutral pions, 
while the middle and right columns are for positive and negative pions. Pions that do not re-scatter in 
the nucleus fall along a diagonal line, while re-scattered pions fall below the line. These correlations 
indicate that while re-scattering is important, pions exiting the nucleus will still carry information about 
their initial state. The top right figure shows the correlation for negative pions arising from CC scatter- 
ing. Production of single negative pions is forbidden due charge conservation, so negative pions will 
only be created by charge exchange from neutral pions. 
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Figure 111: Event display example of charged current event with a muon, proton and neutral pion. This 
event passes a simple topological selection of vr'^ events. The proton has hits close to the vertex, while 
the vr*^ photons convert well away from the vertex. 
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With these simple track multiplicity and vertex cuts, for example, fi^p-K events are identified 
with over a 90% efficiency nearly independent of pion energy. Background processes comprise less 
than 2% of events passing this two charged track cut even without particle identification cuts. Events 
with neutral pions are also identified with good efficiency with these cuts. For example, a topological 
cut designed to isolate a sample of neutral pions accompanied by zero or one nucleons and zero or one 
muons, such as the fi^pn^ event illustrated in Figure fTTTI results in identification of a sample containing 
75% of single neutral pion resonance events with only 10% background. This is in the absence of 
particle identification requirements on the vr^, such as invariant mass reconstruction techniques that are 
described in the MINERz^A proposal^ . 

For charged pions, this topological analysis depends on detecting a nucleon with the pion, which 
biases against pions with low momentum nucleons. Since the full tracking will likely be able to dis- 
tinguish pions by dE/dx, particularly for low energy pions, it should be possible to produce inclusive 
energy and angle spectra of single pions that can be compared directly against pion re-scattering models 
such as Paschos, et.al., described above, or the pion nuclear cascade models in neutrino event gener- 
ators. The shapes of initial and re-scattered pion spectra predicted by NEUGEN are shown in Figure 
11121 Additional capabilities in particle identification, such as identifying production of pions in "for- 
bidden charge" states, e.g., /i~p7r— , will provide a check of pion reinteraction models as these can not 
arise from single pion production on a nucleon through the charged current reaction. Understanding 
these reinteractions is ultimately very important for oscillation analysis as it is these reinteractions that 
limit the ability to transfer knowledge of resonance production from charged lepton data to neutrino 
scattering. 

In addition to identification of resonance production events by simple vertex topology, particle iden- 
tification will be required to differentiate the proton and charged pion in a p-K^ final state, and possibly 
identify pion charge. As described in the proposal, particle identification in MINERz^A will rely on 
measuring discriminating particle characteristics, such as specific energy loss (dE/dx) as well as topol- 
ogy (hadron and electromagnetic showers, decay signatures). As noted above, photons will appear as 
showers disconnected from the vertex, and vr^ will appear as two photons pointing back to the vertex, 
with invariant mass consistent with the vr'^ mass. 

Of particular consideration are low energy particles that stop in the target area without entering 
the magnetic field. We have already demonstrated in the proposal our ability to separate vr^ from 
and protons. However, we would like to be able to distinguish 7r+ from 7r~ in the inner tracker 
alone. Indeed, there are physics processes visible in the highly segmented low density target can allow 
such separation. Shown in Figure fTT3] is the display of a neutrino interaction in the target area before 
GEANT digitization. A muon (red), proton (blue), 7r+ (green) that decays in a vr — > ^ e sequence 
(green-magenta, the muon is not seen). However, a vr^ (black) scatters before stopping without visible 
interaction or may charge exchange with nuclei in the target to produce a neutral pion. 

In order to study the vr^/vr^ differences we generated pion interactions at 0.1, 0.4, 0.8, and 1 GeV/c 
in our target using the MINERz^A simulation. As a case study, the two interactions in Figure fTT4l show 
a 1 GeV/c vr^ in MINERz^A. The final 7r+ decays (vr ^ /i ^ e) with the muon not seen, but with a 
clear positron (magenta). A proton (blue) is also seen coming out of the vertex as weU as two gammas 
(magenta) pointing to the vertex. The it —>■ jj. —>■ e sequence is characteristic to all stopping vr"*". In the 
second interaction, we see a vr^ interaction with charge exchange vr^p —>■ n^n as seen in the central 
detector (after GEANT digitization). The characteristic (two showers magenta) accompanies low 

'see Sections 8.3 and 13.6.2 of the main proposal 
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Figure 112: The top row shows pion energy distributions for events that contain a single tt^ and zero or 
one nucleons, while the bottom row is of events with a single charged pion (positive or negative) and 
a single proton. The left column shows the distributions using known/truth particle identification from 
the monte carlo. The second and third column are events that have been identified using the simple 
topological cuts discussed in the text. The middle column is events that have been correctly identified, 
while the right column is events of other types that have been misidentified as the respective single pion 
event types. The simple cuts correctly identify about 75% of the vr^ events, admitting about a 10% back 
from other types of events. The charged pion efficiency is about 95% with less than 2% contamination 
of background events. Within the statistics of this sample, the efficiencies of event type identification 
are independent of pion energy. 
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Figure 1 13: A 7r+ Decay Chain in MINERi/A 



208 



1 Minerva: x-z | 


















120 








































100 






















80 






















60 






















40 






















20 
°( 


^ , , 1 , , , 1 


1 , 


1 


1 


1 


1 


1 


, 1 






) 20 40 


60 


80 


100 


120 


140 


160 


180 


200 


1 Minerva: x-z | 


















120 








































100 






















80 






















60 






















40 






















20 


^ , , 1 , , , 1 


1 , 


1 


1 


1 


1 


1 


, 1 






) 20 40 


60 


80 


100 


120 


140 


160 


180 


200 



Figure 114: tt"*" and vr Interactions in MINERz^A 
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;ure 115: Electron shower multiplicity and energies for 7r+ and tt Interactions in MINERz^A 
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energy tt~ . These differences in final state can be observed by studying the electron shower multiplicity 
and energy associated with charged vr tracks as shown in Figure 11151 We conclude that although these 
issues need still more study that the MINERz^A inner tracker is capable of statistical vr^ separation even 
without a magnetic field. 
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C.6 MINERz/A and Oscillation Measurements 



Examples of oscillation measurements that will benefit from MINERz^A are the Am^ determination 
by MINOS, and the Vf. oscillation probability measurement by the NuMI Off-axis v^. Appear- 

ance (NOi^A) and the T2K experiments. In chapter 13 of the proposal we presented an estimate for 
how much uncertainties in nuclear effects would contribute to uncertainty in Am^ in MINOS; here 
we update that estimate with an improved treatment of nuclear effects, and consider two consequences 
of these effects separately. After discussing the Am,^ measurement we then describe how NOz^A can 
utilize MINERi^A measurements of cross sections to minimize the systematic error on its oscillation 
probability measurement. In the proposal we described the challenge of measuring the — > os- 
cillation probability due to the uncertainty in backgrounds that must be rejected at a high level, and 
improvements that MINERz^A can provide to neutrino event generators that will used to predict those 
backgrounds. We present here a new analysis (which is also part of the NOz/A proposal, chapter 9) 
which shows quantitatively how uncertainties in cross sections translate to uncertainties in background 
predictions. We also show that information MINERz^A provides will be important regardless of the size 
of oscillations that NOz^A observes. If NOz^A sees no evidence for Vf. appearance then MINERz^A's 
most important contribution will be to improve the determination of the background at the far detector. 
However, if NOz^A does see a signal, then the reduction of the cross section uncertainties plays an even 
larger role in helping NOi^A achieve the best precision on its measurement of the —>■ Vf. oscillation 
probabihty. 

C.6.1 Nuclear Effects and a Am^ measurement 

The key to a precise measurement of Am^ is the ability to measure the oscillation probability as a 
function of neutrino energy. Although MINOS has undergone an extensive program to determine the 
response of its near and far detectors to specific charged particles, it cannot measure the likelihood with 
which those particles are produced in a neutrino interaction. At these low neutrino energies, there are 
two effects that become important, and therefore contribute significantly to the uncertainty in a mea- 
surement of Am^. One effect, which is independent of the target nucleus, results from the rest masses 
of the secondary particles which contribute an important fraction to the reconstructed neutrino energy. 
A measurement of final state particle multiplicities and species as a function of hadron energy, which 
cannot be measured in MINOS, is therefore important for accurate reconstruction of the neutrino energy 
spectra at both MINOS detectors. Secondly, as shown in the proposal, the multiplicity distribution is a 
function of target nucleus, since secondary particles can either scatter in the nucleus, or be completely 
absorbed. Either process results in a change in the visible hadron energy for an event, again contributing 
to the uncertainty in the measurement of the neutrino energy spectrum. 

Figure II 161 shows the changes in the ratio of visible to total neutrino energy for uncertainties in 
absorption (top) and scattering (bottom) separately. In the proposal these two effects were not shown 
separately. Furthermore, nuclear effects were assumed to be the same for carbon as for steel, and a flat 
20% uncertainty on the nuclear plus rest mass effect was assumed. The visible energy is defined as the 
sum of the kinetic energies of all the charged final state particles, plus the total energy for the neutral 
pions, (since it is assumed they deposit all their energy in the form of electromagnetic showers). In the 
top plot, for a steel target, the parameter in the event generator that describes pion absorption is varied 
by three standard deviations. In the bottom plot pion absorption is turned off, and the visible energy 
to total energy ratio is compared for steel, carbon, and lead. Figure [iTt] shows the effect of changing 
the pion absorption on the visible neutrino energy spectra for both near and far detectors (with Am^ 
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Figure 116: Ratio of visible to tine neutrino energy for several different models of nuclear effects. The 
top plot shows the ratio for steel (solid) with the nominal pion absorption, as well as the same ratio for 
the pion absorption cross section varied by plus and minus three standard deviations. The bottom plot 
shows the differences in the average of this ratio for three different target nuclei, where the absorption 
effects are turned off to see more clearly the effects of pion rescattering. 
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Figure 117: (Top left) visible neutrino energy spectrum for the near detector for j/^ CC events, for 
nominal and three sigma high pion absorption in iron. (Top right) visible neutrino energy spectrum for 
the far detector for CC events, again for nominal and three sigma high pion absorption. Am^ in this 
case is assumed to be 2.5 x 10~^eV^. (Bottom) Ratio of changed divided by nominal pion absorption 
model, for the far (solid) and near (dashed) energy spectra, as well as the ratio of far over near. Note 
that the effect of pion absorption cancels somewhat in the ratio of far over near event spectra, but not 
completely. 
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is 2.5 X lO^^eV'^). The effect is also shown (Figure fTTTl bottom) as a ratio for both the near and far 
spectra as well as the effect of the change on the far/near ratio. It is clear that although the visible energy 
difference produces the same kind of effect in near and far spectra, it only partially cancels because of 
the different underlying neutrino spectra. 

Evaluating the appropriate uncertainty in the size of nuclear effects in neutrino scattering is non- 
trivial, because the only precise data on differences between nuclei come from charged lepton scat- 
tering, and one has to use theoretical models to translate the effects from the charged leptons to the 
neutral leptons. The "three standard deviations" in pion scattering were measured in a bubble chamber 
experiment, not on steel, therefore, a theoretical extrapolation must be made for the target in this case. 
If we naively take the differences described above as "one standard deviation" for pion absorption, and 
the differences between steel and carbon or lead as "one standard deviation" for pion rescattering, we 
can determine how this systematic error would compare to the MINOS statistical error. The top plot 
in Figure fTTSl shows how raising and lowering pion absorption would lower and raise, respectively, the 
measured Am^. The central plot shows how assuming nuclear effects for lead or carbon, compared 
to steel, would again lower or raise, respectively, the measured Am^. The bottom plot in Figure fTTSl 
shows the two sets of errors added in quadrature. 

If these errors are appropriate, then they are comparable with the statistical error expected by MI- 
NOS for 7.6 X lO^'' protons on target, as shown in figure [TT9l Also shown on this plot is the previous 
uncertainty that was shown in the proposal, which was evaluated assuming the nuclear effects were 
known to 20%, and which assumed the size of nuclear effects on carbon were the same as those on 
steel. As described in the section on nuclear effects in this addendum, MINERz^A will be able to mea- 
sure these effects in neutrino scattering directly and precisely, for several different targets. 



C.6.2 Measurements of i/^ Ve oscillation Probability 

In order to understand the full importance of MINERi/A cross section measurements for the NOz^A 
experiment, it is helpful to revisit how experiments will determine the i/^ — > oscillation probabilities. 
The number of events in the far detector can be described as 



^far = ^t^P^^i^ Ue)aeeeMf^^ + Bj (68) 

Where (f)^ is the muon neutrino flux at the far detector, P is the oscillation probability, cTg and e^, are 
the electron neutrino cross section and efficiency, respectively, and Mjar is the far detector mass. The 
backgrounds at the far detector, Bjar, can be expressed as 
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Where the notation is the same as equation |68l but ej is the efficiency for a neutrino of type i to be 
misreconstructed as an electron neutrino. Backgrounds come from both muon and electron neutrinos, 
and from several different neutrino interaction channels. Both equation|68land|69lmust be summed over 
neutrino channels, as well as integrals over neutrino energy. 

The error on the oscillation probability, in this simplified notation, can then be expressed as 
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Figure 118: (Top) Increase (decrease) in the measured Am^ due to a decrease (increase) in the pion ab- 
sorption cross section on Iron. Middle plot: Increase (decrease) in the measured Am^ due to assuming 
the target had the nuclear effects of Carbon (Lead) compared to Steel, when pion absorption is turned 
off. (Bottom) The errors due to increases (decreases) in Am^, added in quadrature, as a function of 
Am^. 
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Figure 119: Fractional size of the 90% confidence level region at sin^ 2^23 = 1 for the changes in 
nuclear effects described earlier, assuming the uncertainties for nuclear effects in neutrinos are three 
times the uncertainties coming from the charged lepton measurements. 
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Process 


Statistics 


QE 


RES 


COH 


DIS 


8a jo 




20% 


40% 


100% 


20% 


Signal i^e 


175 (sin2 26*13 = 0.1) 


55% 


35% 


n/I 


10% 


NC 


15.4 





50% 


20% 


30% 




3.6 





65% 


n/I 


35% 


Beamz/g 


19.1 


50% 


40% 


n/I 


10% 



Table 21: List of the signal and background processes than can contribute events in the NOi^A far 
detector, for a 50kton detector located 12km from the NuMI axis, 820km from Fermilab, assuming 
a Am^ of 2.5 x V^^^eV^ . Also given are the cross section uncertainties on those processes before 
MINERi/A runs. 



The two terms in equation ^Ol suggest two regimes: in the case where the number of events in the 
far detector is comparable to the background prediction, the error on the probability is dominated by a 
combination of statistics and the uncertainty on the background prediction. The background prediction 
uncertainty is dominated by the uncertainty in the background process cross sections and efficiencies. 
In the other extreme, where the number of events is dominated by the signal events, the uncertainty on 
the probability comes from the statistics, and the uncertainties on the signal channel cross sections and 
efficiencies. 

For the NOz^A detector simulation running at an off axis location 12km from the NuMI beamline 
and 820km from the source at Fermilab, the signal and background statistics for the nominal 5 year 
run are given in table |^ As has been described in more detail in both the NOi^A and MINERz^A 
proposals, the three categories of backgrounds are V(.\ originally produced in the neutrino beam, the 
neutral current events with energetic neutral pions which can fake electrons, and z/^ charged current 
events, where the final state muon is low energy and the event contains a high energy neutral pion. Also 
given in Tablel^Tlare the fractions that each neutrino interaction process contributes to the events of that 
type that pass all cuts, as well as the cross section uncertainty on that process. 

Without a near detector, the errors from cross sections, for the case that there are no oscillations, 
are 16%, which is equivalent to the statistical error for that case. For the case of mixing at the level 
indicated in the table, the statistical error on the probability would be 8%, while the errors from cross 
section uncertainties alone would be 31%. 

In the case where there is a near detector that is identical to the far detector, one can try to cancel out 
these uncertainties. Consider first the prediction of the background events. The events in an identical 
near detector that pass the same analysis cuts as those made at the far detector can be described as 

^near — ^i=e,fi4'i^i^i-^near (71) 

And then one can use the simulation to predict the number of backgrounds at the far detector by the 
following equation: 

Mfar 

^far = ^nearj7 R (72) 



mear 



where. 



R = ""--^^J-r-^'^ (73) 
'^i=e,ii4'i,nearO'i^i 
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For different near detector off axis angles, there are different fractions of background events that 
pass all cuts, but at no value of the angle is the mix of backgrounds the same as that in the far detector. 
Figure fT20l shows the fractional change in the variable R defined above due to the cross section errors 
listed in table |^ This translates directly into an uncertainty on the far detector background prediction 
due to current cross section uncertainties. 

Note that for low off axis angles the systematic error due to the charged current uncertainties 
is minimum, yet the error due to the neutral current and electron neutrinos is maximum there. Figure 
11211 shows the above errors added in quadrature, as a function of near detector off axis angle. Note 
that at best the cross section uncertainties can be reduced from 16% with no near detector to about 
10% with an identical near detector, but given that this is only one of several systematic errors in the 
NOi^A experiment, it is an unacceptably large fraction of the total error. The MINERz^A experiment 
can significantly reduce the errors due to the charged current processes, because although it is on the 
NuMI axis where the neutrino flux is different from the off-axis fluxes, for the charged current processes 
the final state energy is close to the neutrino energy, and the flux prediction from the hadron production 
combined with NuMI horn B-field measurements means that charged current cross section channels 
will be measurable at the 5% level overall. So with the presence of MINERz/A the error due to cross 
section uncertainties on the background at the far detector can be reduced from 10% to better than 
5%. Furthermore, by measurements of charged current coherent pion production on axis one can infer 
the neutral current coherent pion process off axis using theoretical models constrained by the charged 
current process. 

For the case where the number of signal events is well above the number of background events, 
the challenge to keep the uncertainties due to cross section errors low is even harder, since in that 
case the composition of events near to far is even more different than it is in the case of no signal 
events. Furthermore, because the total number of events is higher, the improved statistical precision 
will required an even more precise far detector prediction. 

As a final note, it should be pointed out that the only cross section errors considered here are 
uncertainties in the overall levels of each of the processes. There is an additional uncertainty in the 
energy dependence of these processes, which will again contribute uncertainties in the far detector 
predictions because of the differing spectra. 

C.6.3 Summary 

This section has quantified how MINERz^A's cross section measurements will have important conse- 
quences on both the current and future generations of neutrino experiments with the NUMI beamline. 
These improvements to MINOS and NOz^A measurements are important regardless of the ultimate value 
of Am^ that MINOS measures, and regardless of the size of the — > oscillation probability that 
NOi^A ultimately measures. 
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Figure 120: The fractional change in the far detector background prediction coming from an identical 
near detector, as a function of near detector off axis angle, for a far detector located 12km off the NuMI 
axis and 820km from Fermilab. The top plot shows the fractional change when the neutral current 
cross sections are varied by their uncertainties, the middle plot shows the fractional change when the 
charged current cross sections are varied by their igg^rtainties, and the bottom plot shows the fractional 
change when the v^. charged current cross sections are varied by their uncertainties. 
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Figure 121: The fractional error in the background prediction at the far detector from uncertainties 
in each process (Quasi-elastic, resonance, deep inelastic scattering, and neutral current coherent vr^ 
production), added in quadrature for each source (NC, v^CC, beam Vf.), plotted as a function of near 
detector off axis angle. 
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